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STOCHASTIC HOMOGENIZATION OF NONCONVEX UNBOUNDED 
INTEGRAL FUNCTIONALS WITH CONVEX GROWTH 

MITIA DUERINCKX AND ANTOINE GLORIA 


We consider the well-travelled problem of homogenization of random integral functionals. When 
the integrand has standard growth conditions, the qualitative theory is well-understood. When it 
comes to unbounded functionals, that is, when the domain of the integrand is not the whole space 
and may depend on the space-variable, there is no satisfactory theory. In this contribution we develop 
a complete qualitative stochastic homogenization theory for nonconvex unbounded functionals with 
convex growth. We first prove that if the integrand is convex and has p-growth from below (with 
p > d, the dimension), then it admits homogenization regardless of growth conditions from above. 
This result, that crucially relies on the existence and sublinearity at infinity of correctors, is also new 
in the periodic case. In the case of nonconvex integrands, we prove that a similar homogenization 
result holds provided the nonconvex integrand admits a two-sided estimate by a convex integrand (the 
domain of which may depend on the space variable) that itself admits homogenization. This result 
is of interest to the rigorous derivation of rubber elasticity from polymer physics, which involves the 
stochastic homogenization of such unbounded functionals. 
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1. Introduction 

Let O be a bounded Lipschitz domain of d,m > 1. We consider the well-travelled problem of 
homogenization of random integral functionals : VL^’P(0;K™') —)• [0, oo] given by 

leiu) = j W , Vu(x)^ dx, 

where PL is a random Borel function stationary in its first variable that satisfies for almost every 
?/ G and all A G the two-sided estimate 

^\K\P-C<V{y,K) < W(y,A) < C{l + V{y,k)), (1.1) 

for some C > 0, p > 1, and a random convex Borel function 1/ : x [Q,oo]. The originality 

of the growth condition we consider here is that V{y, ■) may take infinite values and that its domain 
may depend on y, so that the domain of the homogenized integrand W (if it exists) is unknown a 
priori. The motivation for considering such a problem comes from the derivation of nonlinear elasticity 
from the statistical physics of polymer-chain networks, cf. [4, 26, 20]. Indeed, the free energy of the 
polymer-chain network is given by two contributions: a steric effect (for which proving homogenization 
is one of the most important open problems of the held), and the sum of free energies of the deformed 
chains. The free energy of a single chain is a convex increasing function of the square of the length 
of the deformed polymer-chain, which blows up at hnite deformation. The corresponding problem in 
a continuum setting would be the homogenization of the nonconvex integrand 

W(y, A) = Viy, A) + g(det A) < C(1 + V{y, A)), (1.2) 

where V is an infinite-valued convex stationary ergodic integrand whose domain depends on the space 
variable, and ^ is a convex function (in this paper we assume that g is controlled by V, which 
unfortunately rules out the hnite compressibility of matter). 

Homogenization of multiple integrals has a long history, and we start with the state of the art when 
V and W are periodic in the hrst variable: 

(i) The hrst contribution (beyond the linear case) is due to Marcellini [34], who addressed the 
homogenization of convex periodic integrands satisfying a polynomial standard growth condition, 
that is (1.1) for V{y,k) = |A|P. 

(ii) Marcellini’s result was then generalized to nonconvex periodic integrands satisfying a polynomial 
standard growth condition, by Braides [11] (which covers in addition almost-periodic coefficients) 
and Muller [36, Theorem 1.3], independently. 

(hi) In [36, Theorem 1.5], Muller also addressed the case of a convex periodic integrand satisfying 
a convex standard growth condition (1.1) for V{y,K) = 14(A) with V : —)• M+ a convex 

hnite-valued map, and p > d. 

(iv) In [12, Chapter 21], Braides and Defranceschi treated the case of nonconvex periodic integrands 
(see also [18] in the convex case) satisfying (1.1) where V is convex periodic and satishes the 
polynomial non-standard growth condition 

< V{y,A) < C(l + |Ar) 

for some q < p* (with p* the Sobolev-conjugate of p > 1), and the doubling property 

V{y,2A) < C{l + V{y,A)). 
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(v) In [13], given a collection of well-separated periodic inclusions, Braides and Garroni treated the 
case of nonconvex periodic integrands satisfying a polynomial standard growth condition as well 
as the (strong) doubling property 

W{y,2A) < CW{y,A). 

outside the inclusions, but only satisfying inside the inclusions a convex standard growth condi¬ 
tion (1.1) for some possibly unbounded map V{y,A) that is convex in the A-variable. 

(vi) More recently Anza Hafsa and Mandallena studied in [7] the homogenization of quasiconvex 

periodic integrands satisfying a standard (unbounded) convex growth condition, that is, (1.1) for 
V{y,A) = 1^(A) with V : [0,oo] a convex infinite-valued map such that t^(A) > |A|^, 

and with p > d. Note that in this case the domain is fixed. 

When W is random, the results are more sparse: 

(vii) The first contribution beyond the linear case is due to Dal Maso and Modica, who addressed 
the homogenization of convex random stationary integrands satisfying a polynomial standard 
growth condition [19], generalizing Marcellini’s result (ii) to the random setting. 

(viii) Messaoudi and Michaille later treated the homogenization of quasiconvex stationary ergodic 
integrands satisfying a polynomial standard growth condition [35], following Dal Maso and 
Modica’s approach. 

(ix) In their monograph on homogenization, Jikov, Kozlov and Oleinik treated the case of convex 
stationary ergodic integrands satisfying a polynomial non-standard growth condition 

1|A|P-C < W{y,A) < C(1 + |A|'?) 

for some q < p* (with p* the Sobolev-conjugate of p > 1), see [31, Chapter 15]. 

For scalar functionals, that is, when m = 1, results are much more precise, and we refer the reader 
to the monograph [17] by Carbone and De Arcangelis (which is however only concerned with the 
periodic setting) and [15, 16]. When the domain of V{y, •) is not the whole of the authors call 
an unbounded functional. The main technical tool for scalar unbounded functionals is truncation of 
test-functions (see also [31, Section 15.2]), which cannot be used for systems in general (see however 
the end of this introduction). In particular, such truncation arguments replace the Sobolev embedding 
we shall use for systems and allow one to relax the assumption p > d for scalar problems. 

In this contribution we give a far-reaching generalization of (i)-(iv) and (vi)-(ix) for systems in the 
random setting, by relaxing the assumption that the domain of K(y,-) in (1.1) is independent of y. 
Our contribution also generalizes (v) by relaxing all geometric assumptions. We argue in two steps. 
For convex integrands, our result shows that homogenization holds without any growth condition 
from above (cf. Theorem 2.2 for Neumann boundary conditions, and Corollary 2.4 for the more 
subtle case of Dirichlet boundary conditions), so that we may homogenize the bound V itself. We 
proceed by truncation of the energy density (following the approach by Muller in [36]), and first prove 
in Proposition 3.8 that homogenization and truncation commute at the level of the definition of the 
homogenized energy density. The proof relies on the existence of correctors with stationary gradients 
for convex problems and exploits quantitatively their sublinearity at infinity, see Lemma 3.4 (which 
is a substantial difference between the periodic and random cases, and makes the latter more subtle). 
The second main technical achievement is the construction of recovery sequences in Proposition 3.10 (a 
gluing argument based on affine boundary data trivially fails since the domain of V (y, •) depends on y 
— this difficulty is already present in the periodic setting, and prevents us from using the standard 
homogenization formula with Dirichlet boundary conditions). In the case of nonconvex integrands 
with a two-sided convex estimate (1.1), we show in Theorem 2.6 that homogenization reduces to the 
homogenization of the convex bound V. The first obstacle in this program is the definition of the 
homogenized energy density itself. Indeed, in the absence of correctors (which is a consequence of 
nonconvexity), one usually defines the homogenized energy density through an asymptotic limit with 
linear boundary data on increasing cubes. As pointed out above, such an approach fails in general 
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when the domain of the integrand is not fixed. Instead, in Lemma 4.1, we use the (well-defined) 
corrector of the associated convex problem as boundary data for the nonconvex problem on these 
increasing cubes. Next we argue by blow-up in Proposition 4.3 for the P-liminf inequality (following 
the approach by Fonseca and Muller in [25]), and make a systematic use of the corrector of the 
convex problem to control the nonconvex energy from above. Then, for the F-limsup we argue in 
Proposition 4.6 by relaxation (following the approach introduced by Fonseca in [24] for relaxation 
and hrst used in homogenization by Anza Hafsa and Mandallena in [7]), making a similar use of the 
corrector of the convex problem in the estimates. 

To conclude this introduction, let us go back to our motivation, that is, the homogenization of (1.2). 
On the one hand, we have reduced the homogenization for such integrands to the homogenization for 
the convex integrand V. On the other hand, we have proved homogenization for convex integrands 
without growth condition from above, and therefore proved homogenization for (1.2). In the specihc 
setting of (1.2), we can sharpen the general results described above, by simplifying the definition of the 
homogenized energy density W, cf. Corollary 2.8. We believe a similar approach can be successfully 
implemented in the discrete setting considered in [4] for the derivation of nonlinear elasticity from 
polymer physics. 

In the scalar case m = 1, combining our approach with truncation arguments, we may further 
refine our general results, in particular relaxing the condition p > d. Our approach then improves 
(and extends to the stochastic setting) some scalar results of [15, 17]. 

The main results are given in Section 2. The proofs of the results for convex integrands are displayed 
in Section 3, whereas Section 4 is dedicated to the proofs for nonconvex integrands. In Section 5 
we turn to various possible improvements of our general results under additional assumptions. In 
the appendix we prove several standard and less standard results on approximation of functions, 
measurability of integral functionals, and Weyl decompositions in probability, that are abundantly 
used in this article. 


2. Main results 

Let (f4,T', P) be a complete probability space and let r := {Ty)y^-^d be a measurable ergodic action 
of (M'^,-|-) on (n,T', P), that are hxed once and for all throughout the paper. We denote by E the 
expectation on with respect to P. 

Consider a map VF : x x fl —)■ [0, oo] that is r-stationary in the sense that, for all A £ 

all w £ fl, and all y,z € 

W{y,A,T-^u}) = W{y + z,A,u}), (2.1) 

and assume that W{y, -, 00 ) is lower semicontinuous on for almost all y,a;. We also assume in 

the rest of this paper that, for almost all a;, the map y 1 —)• W{y,A + u{y),uj) is measurable for all 
u £ Mes(M'^; M™^'^), and that, for almost all y £ M'^, the map u W{y^A + v{uj),lo) is measurable 
for all V £ Mes(f4; M™’^'^). Continuity in the second variable and joint measurability (in which case 
W is called a Caratheodory integrand) would ensure these properties; weaker sufficient conditions for 
these properties are given in Appendix A.l. Such integrands W will be called r-stationary normal 
random integrands. 

We further make the following additional measurability assumption on W: 

Hypothesis 2.1. For any jointly measurable function / : x 12 —)• M and any bounded domain 

O C 


W{yfe, f{y, u) + Vu{y),uj)dy 


is T-measurable on fl. 


oj eA inf / 
«GVKo’^(0;K"‘)iO 


□ 
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As discussed in Appendix A.4.2, this last hypothesis is always satisfied if W is convex in the second 
variable, and more generally if it is sup-quasiconvex (in the sense of Definition A. 11), which includes 
e.g. the case of a sum of a convex integrand and of a “nice” nonconvex part. 

For any bounded domain O C we define the following family of random integral functionals, 
parametrized by e > 0, 

4(-,-;0) : X [0,oo] : {u,oj) ^ h{u,uj-0) := [ W{y/e,Vu{y),oj)dy. (2.2) 

Jo 

The aim of this paper is to prove homogenization for Ig as e 0 under mild growth conditions on W, 
which we formulate in terms of F-convergence for the weak convergence of iy^’^(0;M™') (for some 
p > 1). When A i—)• W{y, A,uj) is convex for almost all y,uj, we say that VF is a r-stationary convex 
normal random integrand, and shall use the notation V and instead of W and 1^, that is, for every 
bounded domain O C and e > 0, 

j£(-,-;0) : fy^’^(0,M™) X [0,oo] : {u,uj) J,{u,uj]0) ■.= f V{y/s,Vu{y),uj)dy. (2.3) 

Jo 

The notation W and Je will be used for nonconvex integrands. We start our analysis with the 
case of convex integrands, then turn to nonconvex integrands, present an application to nonlinear 
elasticity, and conclude with a discussion of several possible improvements of these general results 
under additional assumptions. 

2.1. Convex integrands. In this subsection we state homogenization results for with (essentially) 
no growth condition from above. We start with Neumann boundary conditions, and then address the 
more subtle case of Dirichlet boundary conditions. 

2.1.1. Homogenization with Neumann boundary conditions. Our first result is as follows. 

Theorem 2.2 (Convex integrands with Neumann boundary data). Let C : X x D —>■ 

[0, oo] be a r-stationary convex normal random integrand that satisfies the following uniform coercivity 
condition: there exist C > 0 and p > d such that for almost all oj and y, we have for all A, 

^\A\P-C < V{y,A,u:). (2.4) 

Assume that the convex function M := supess^ ^^ yiU) has 0 in the interior of its domain. Then, 
for almost all oj ^ Q and all bounded Lipschitz domains O C the integral funetionals 
T-converge to the integral functional J{-]0) : IF^’^(0;M™’) —[0, oo] defined by 

J{u-0)= [ V{Vu{y))dy, 

Jo 

for some lower semicontinuous convex function V : —)• [0,oo] characterized by the following 

three equivalent formulas: 

(i) Formula in probability: for all A G 

F(A)= inf E[C(0,A + ff,.)], (2.5) 

where the spaee of mean-zero potential random variables TpQ^(D) is recalled in Seetion 3.1.2. 

(a) Dirichlet formula with truncation: for any increasing sequence t ^ of r-stationary convex 
random integrands that satisfy the standard p-growth condition 

i|Ar - C < V\y,A,u:) < Ck{l + |A|?>) (2.6) 

for all y,u},A, and some sequenee < oo, we have for almost all uj, all A, and all bounded 
Lipschitz domains O C 

C(A) = limlim inf V’^{y,A-\-V(j){y),uj)dy. 

fctoo eto Jo/e 


(2.7) 
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(Hi) Convexification formula: for all A and all bounded Lipschitz domains O C we have, for 
almost all oj, 

y (A) = lim lim inf -r V{y,tA + V4>{y),u})dy. (2.8) 

t\l e4,0 </)ewi’P(0/£;K’") Jo/e 
io/e 

As a consequenee of convexity, the limit f f 1 can be omitted when A ^ 5domy. □ 

Comments are in order: 

- The limit t f 1 cannot be omitted in (2.8) in general for A S ddomV. Indeed, let V coincide with 

a convex map V : —)■ [0, oo] with a closed domain, and which is not lower semicontinuons 

at the bonndary of its domain. In the interior of domi/, V coincides with V. However, since V 
is necessarily lower semicontinnous on its domain, it cannot coincide with V on cldomH. 

- In the proof we take (2.7) as the dehning formnla for V, following the approach by Miiller in [36]. 
Formula (2.5) is interesting in two respects: hrst, it is intrinsinc (no approximation is required), 
and second it is an exact formula (there is no asymptotic limit involved). The equivalence of 
both formulas, which can be interpreted as the commutation of trunction and homogenization, 
is the key to the proof of the T-convergence result. 

We may extend Theorem 2.2 in two directions: 

- The extension of Theorem 2.2 to the case of domains with holes (or more generally to soft 
inclusions, for which the coercivity assumption (2.4) does not hold everywhere) is straighforward 
provided we have a suitable extension result. When holes are well-separated, such extension 
results are standard (see e.g. [31, Sections 3.1]). For general situations however, this can become 
a subtle issue (see in particular [31, Sections 3.1 and 3.5]). In the particular case of the periodic 
setting, there is a very general extension result [1], which is used e.g. in [13]. 

- The assumption p > d is crucial in the generality of Theorem 2.2, which is used in the form of the 

Sobolev embedding of VF^’^(0;M™') in L°°(0;]R™'). In the case 1 < p < d, the conclusions of the 
theorem still hold true provided that V{y, A, cj) < M (A) for some convex function M : M 

that satisfies the growth condition lim|yY|_j,oo M(A)/|A|'? < oo for q = dp/{d — p) if p < d or for 
some q < oo A p = d. In the scalar case m = 1, the use of the Sobolev embedding can be avoided 
by a truncation argument, as explained in Corollary 2.9 below, see also [17, 15[. 

2.1.2. Diriehlet boundary conditions. We now discuss the homogenization result in the case of Dirich- 
let boundary conditions (the case of mixed boundary data can then be dealt with in a straightforward 
way, and we leave the details to the reader). A first remark is that Diriehlet data have to be well- 
prepared, as the following elementary example shows. 

Example 2.3. Consider random unit spherical inclusions centered at the points of a Poisson point 
process, choose the integrand V to be equal to |A|P outside the inclusions and to have a bounded 
domain V C inside the inclusions. Given a (nonempty) bounded open set O C for almost 

all CO, the realization of the inclusions corresponding to lo intersects d{0/e) for infinitely many e > 0, 
and hence limsup^, Js{u + A ■ x, 0 J-,O) = oo for all u G Hq’^(O; M"*), due to the Diriehlet boundary 
condition. In contrast, if the intensity of the underlying Poisson process is not too big, it is easily 
seen that the homogenized integral functional J dehned in Theorem 2.2 is finite-valued. This proves 
that, for all A ^ 2? and almost all co, Je{- + A - x, (o; O) cannot P-converge to J(-; O) on ITg^’^(0; M™), 
due to the intersection of some rigid inclusions with the boundary of the domain where the Diriehlet 
condition is imposed. □ 

We propose two ways to prepare Diriehlet data: 

- by relaxing the boundary data so that the energy remains finite for all e > 0 while ensuring that 
the boundary data are recovered at the limit e 0 — we call “lifting” this procedure; 

- by replacing the integrand H by a softer integrand on a neighborhood of the boundary where 
the boundary condition is imposed — we call this a “soft buffer zone”. 
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Corollary 2.4 (Convex integrands with Dirichlet boundary data). Let V, M, J^, and J be as in 
Theorem 2.2 for some p > d. Then, for almost all uj ^ Ll and all hounded Lipschitz domains O C 
we have 


(i) For all boundary data u G VC^’P(0; M™') such that J{au;0) < oo for some a > 1, there exists a 
lifted sequence (rte)e with ^ u in such that we have on VCq’^(0; 


J{- + u\0) = F-lim r-liminf J^{-+ tu£,(jj]0) 

t\l etO 

= r-lim r-limsup Je{-+ tUe,0J]O). 

*ti £4,0 

In particular, 


inf J{v + u]0) 


lim lim inf inf (v + tue ; O) 
m efo veWg’^iO) 

limlimsup inf J^{v + tUs]0). 
Cl £4.0 v^Wq’P{0) 


If u satisfies the additional condition jQM{Vu{y))dy < oo, then we may take = u, and if 
this condition is strengthened to M{a'Vu{y))dy < oo for some a > 1 then the limit tfl can 
be omitted. 

(a) For all boundary data u G M”*) such that J{u;0) < oo, we have on VbQ^’^(0; M”^); 

J{- + u]0) = F- lim F-liminf J'^{-+ tu,(w,0) 

ttl,r)4-0 ei'O 

= F- lim F-lim sup J^{-+ tu,uj',0), 
tfl.rjto £4,0 


where is the following modification of on an rj-neighborhood of dO: 






In particular. 


Jo 

f V{yfe, A, w), if dist(y, dO) > p; 
[|A|P, z/dist(2/,(90) < r/. 


inf J{v + u;0) 

veWg-PiO) 


lim lim inf inf J^{v + tu-,0) 

ttl,r;4.0 £4,0 veW^’^iO) 

lim lim sup inf Jf(v + tu:0). 
ttl,r;4-0 £ 4,0 v£Wq'^{0) 


(2.9) 


If u satisfies the additional condition J{au-, O) < oo for some a > 1, then the limit t f 1 can be 
omitted. □ 


Comments are in order: 

- The results of Corollary 2.4 are not completely satisfactory. Indeed, if one makes a diagonal 
extraction of t and p wrt e to obtain a F-convergence in e only, then the extraction for the 
F-limsup depends on the target function v-\-u and not only on the boundary data u as we would 
hope for. This dependence is however restricted to a dilation parameter only in the case of the 
lifting. In the case of the buffer zone, the result can be (optimally) improved to % = 9e for 
any fixed 6 > 0, under some additional structural assumption in the form of the existence of 
stationary quasi-correctors (see Proposition 5.1). For specific examples for which this assumption 
holds, see Corollary 2.10 below. 

- In the specific situation when domV = domM (this is trivially the case when the domain is fixed, 
i.e. domV{y,-,oj) = domM for almost all y,uj), then the strong assumptions M{'S/u{y))dy < 
oo and Jq M[aS/u{y))dy < oo (for some a > 1) in part (i) of the statement reduce to the simpler 
assumptions J(u;0) < oo and J(au]0) < oo (for some a > 1), respectively. In particular, in 
that situation, the lifting can always be chosen to be trivial, := u for all £. 
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- In [13] (see also [12, Chapter 20]), Braides and Garroni prepare the boundary data in a different 
way in the specific case of stiff inclusions. In particular they introduce an operator which 
acts on functions u as follows: on each stiff inclusion R^{u) has value the average of u on the 
considered stiff inclusion, away from all inclusions R^{u) coincides with u, and in between R‘^{u) 
is an interpolation between u and the average of u on the inclusion. Such a construction can 
be used here as well, but seem to admit no natural generalization in other settings than stiff 
inclusions. 


2.2. Nonconvex integrands with convex growth. In the case when W is nonconvex and admits 
a two-sided estimate by a convex function (which may depend on the space variable), we show that 
a F-convergence result similar to the convex case holds. Before we precisely state this result, let us 
recall the notion of radial uniform upper semicontinuity (which is trivially satisfied by convex maps). 


Definition 2.5. A map Z : [0, oo] is said to be ru-usc (i.e. 

semicontinuous) if there is some a > 0 such that the function 


A|(t) 


Z(tA) - Z(A) 

AsdomZ ® “1“ ^\A) 


radially uniformly upper 


satisfies limsup^^^ A^(t) < 0. A r-stationary normal random integrand W is said to be ru-usc if there 
exists a r-stationary integrable random field a : x Id —)• [0, oo] such that the function 


A^(t) :=supess 


sup ess sup 

AEdomVF(j/,-,a;) 


W{y,tA,uj) -W{y,A,uj) 
a{y,uj) + W{y,A,uj) 


satisfies limsup^^i A^(t) < 0. 


□ 


The following result is a far-reaching generalization of [36, Theorem 1.5] to a wide class of random 
and nonconvex integrands; it is also a substantial extension of [7, Corollary 2.2]. 


Theorem 2.6 (Nonconvex integrands with convex growth). Let IT : x x —)• [0,oo] be 

a (nonconvex) ru-usc r-stationary normal random integrand satisfying Hypothesis 2.1. Assume that, 
for almost all oj, y, and for all A, 


V{y,A,co) < W{y,A,Lo) < C{1+ V{y,A,io)), 


( 2 . 10 ) 


for some C > 0 and some r-stationary convex random integrand 17 : x X 11 —)• [0, oo] that 

satisfies the assumptions of Theorem 2.2 for some p > d. Then, for almost all co & LI and all hounded 
Lipschitz domains O C the integral functionals T-converge to the integral functional 

[0,oo] defined by 

= [ W{Vu{y))dy, 

Jo 

for some ru-usc lower semicontinuous quasiconvex function W : —>■ [0, oo] that satisfies V < 

W < C{1-LV), where V is the homogenized integrand associated with V by Theorem 2.2. In addition, 
the results stated in Corollary 2.4 for also hold for fi. 

For all A G let g\ be the potential field in probability minimizing E[17(0, A -|- •)] (cf. (2.5)j, 

and note that x i-A is a gradient field on for almost all td G H, which we denote by 

'\Apa{x,uj). The homogenized integrand W is characterized for all A G by 


1T(A) = liminf liminf lim inf 

ttl A'^A elO ^eW'o’P(o/£;R’") 


O/e 


W{y,tA' + tVq:)A'{y,uj) + Vv{y),uj)dy. 


( 2 . 11 ) 


for any bounded Lipschitz domain O C and almost every a; G H. 


□ 


In general, we do not expect that the limits 11 1 and A' —)• A can be dropped in (2.11), see however 
Corollary 2.8 below under Hypothesis 2.7. 
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2.3. Application to nonlinear elasticity. In the example from the statistical physics of polymer- 
chain networks, the integrand has the specific decomposition (1.2). Moreover, the nonconvex part of 
the integrand satisfies the following assumption, which in particular implies that W satisfies Hypoth¬ 
esis 2.1 (see indeed Lemma A. 12), as well as the ru-usc property. 

Hypothesis 2.7. There exists a r-stationary convex map H X x 11 —)• [0,oo] and some 

p > 1 such that 

W{y, A, a;) = V{y, A, a;) + W^'^{y, A, cu), 

where —)• [0, oo] is a (nonconvex) r-stationary normal random integrand satisfying 

the p-th order upper bound and the local Lipschitz conditions, for some C > 0, 

(i) for almost all y,uj, for all A, 

W^\y,A,u:) < CdA^ + 1); 

(ii) for almost all y,uj, for all A, A', 

- W^'^{y,A',u;)\ < C{1 + jAr^ + |A'r')|A - A'|. 

□ 

Under this assumption, the following variant of Theorem 2.6 holds and yields in particular a simpler 
formula for the homogenized energy density. 

Corollary 2.8. Let W : x x H ^ [OjCo] be a (nonconvex) r-stationary normal random 

integrand satisfying Hypothesis 2. 7. Assume that, for almost all lj, y, and for all A, 

V{y,A,oj) < W{y,A,vj) < C{1+ V{y,A,u)), (2.12) 

for some C > 0, where U : x x U —>• [0, oo] satisfies the assumptions of Theorem 2.2 for 

some p > d. Then the conclusions of Theorem 2.6 hold true. In addition we have for all A S 
for almost all tv 

1^(A) = lim inf / W{y,A + VipA{y,u;)+Vv{y),uj)dy, (2.13) 

where ipA is ns in Theorem 2.6. □ 

The behavior of W close to the boundary of its domain is of particular interest. From a mechanical 
point of view, in the case of (1.2) with V satisfying U(y,A„) —>■ oo as dist(A„, (9doml/(y, •)) —t 0, 
we expect that W{An) —>■ oo as dist(A„, (9domlU) —)• 0. Except in a few explicit examples, we do 
not know under which condition and how to prove such a property. In the case of the geometry that 
saturates the Hashin-Shtrikman bound (obtained by a stationary and statistically isotropic Vitali 
covering of by balls, cf. [31, Section 6.2]), where each ball is the homothetic image of a unit 
ball with a given spherical inclusion, the corrector gradient field is explicit when A is a dilation, and 
t e^W{t Id) indeed blows up when t approaches the boundary of the domain. 

2.4. Some improved results. The general results above naturally call for various questions con¬ 
cerning possible improvements: 

- What about the subcritical case 1 < p < d? 

- What is the minimal size % of the soft buffer zone needed in the presence of Dirichlet boundary 
conditions (see Corollary 2.4(ii))? Under which conditions can we take = 6e for some constant 
0 > 0? _ 

- What about periodic boundary data? In particular, can one approximate V by periodization (in 
law) in the spirit of [22]? 

These three questions are partially addressed below under various additional assumptions. 
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2.4.1. Subcritical case 1 < p < d. The first improvement concerns the growth condition from below 
in Theorem 2.2. It is relaxed here to any p > 1 in the convex scalar case m = 1 under the addi¬ 
tional assumption that V has fixed domain. The idea is to avoid the use of Sobolev embedding by 
using suitable truncations (in the spirit of e.g. [17, proof of Lemma 13.1.5]), which are indeed only 
available in the scalar case with fixed domain. We recover in particular in this way the results of [17, 
Section 13.4] in Sobolev spaces. 

Corollary 2.9 (Subcritical case). Let V and M satisfy the assumptions of Theorem 2.2 in the scalar 
case m = 1, for some p > 1. Also assume that d.ovaV{y, -^oj) = domM for almost all y,uj. Then the 
conclusions of Theorem 2.2 and Corollary 2.4 hold true for this p > 1. □ 

2.4.2. Minimal soft buffer zone. The second improvement concerns the size of the buffer zone for 
Dirichlet boundary data, at least for affine target functions. The minimal size % = 6e, for any 
constant 0 > 0, is achieved under the technical structural assumption that stationary quasi-correctors 
exist, cf. Proposition 5.1. Understanding the validity of this assumption in general seems to be a 
difficult question of functional analytic nature. It is trivially satisfied in the periodic case. It is also 
valid provided that truncations are available, which holds in the scalar case with fixed domain. 

Corollary 2.10 (Minimal soft buffer zone). Let V, and W,Is,I be as in Theorems 2.2 and 2.6 

for some p> 1. Also assume that one of the following holds: 

(1) p > d, and, for all A,u}, U(-,A,a;) and W{-,A,uj) are Q-periodic; 

(2) m = 1, and domU(y, -, 0 ;) = domM is open for almost all y,U!. 

Then, for all A, for almost all a; G 11, 

U(A)=lim inf -I V^’^^{y,A + Vv{y),uj)dy, 

eio vew^’^iO) Jo 

where is defined as in (2.9) with rj = 9e. The same result also holds for V,Js,J replaced by 

W,Ie,I (forp>d). □ 

2.4.3. Approximation by periodization. The last improvement concerns the approximation of V by 
periodization. As for Dirichlet boundary conditions, periodic boundary conditions require “well- 
prepared data”. In this case, the well-preparedness is formalized in terms of periodization of the law 
of the energy density (as opposed to the more naive periodization in space, cf. Figure 1): 

Definition 2.11. A collection (U'^)/j>o of random maps x x D —>• [0, oo] is said to be 

an admissible periodization in law for V if 

(i) Periodicity in law: for all R > 0, is Qr = [—■§, Y)^'P6i'iodic and is stationary with respect 
to translations on the torus Tji = (M/i?Z)'^ (more precisely there exists a measurable action 

of Tjj on Q such that, for all y = yi + y 2 with yi G and y 2 G Qr, we have V^{y, -, 0 :) = 

V^{y2,;u)=V^{0,;r^y^u;)); 

(ii) Stabilization property: for all 9 G (0,1) and for almost all lo, there exists Re{oj) > 0 such that 

for all R > Reiuj) we have U'^(-, •,a;)|Qg^xR-x‘i = ^)lQeflxR-xd. □ 

Remark 2.12. It may seem quite unnatural to make in this definition the action depend on 
R. Another definition would consist in keeping the same integrand V and action r but letting the 
probability measure vary according to Fr := (rR)*P, which amounts to defining V^{y, A,uj) = 
V{y, A,Tr(uj)), for some projection Tr : D —)• D with suitable properties such as the r-invariance of 
the set Tr{LI). As shown in the examples below, Tr will typically be the (natural) Q^-periodization 
of an underlying Poisson process, and Tr(LI) will be the set of Q/j-periodic realizations (that is of 
P-measure 0). Note that this is a particular case of the formalism of Definition 2.11 above, for the 
conjugated action := Ty oTr. 
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Figure 1. For Poisson random inclusions: periodization in space (left) versus peri¬ 
odization in law (right), four periods. 


The main goal here is to prove an approximation formula for V and W by periodization in law, 
that is an asymptotic formula of the following form: for all A S for almost all w. 


F(A) ~ lim inf 

^too u^W^giQu) 


/ 

J Qr 


V^{y,K + Vu{y),uj)dy, 


(2.14) 


and similarly for W. Such a formula would be of particular interest for numerical approximations 
of V and W. The notion of periodization in law was introduced in [22]. It crucially differs from 
“periodization in space”, which would consist in replacing the argument of the limit in (2.14) by 

inf f V(y,A + Vu(y},uj)dy. 
u&WpfriQn) JQr 

The difference between both types of periodization is illustrated on Figure 1 for Poisson random 
inclusions. Empirical results in [22] tend to show that periodization in law is more precise than pe¬ 
riodization in space for the approximation of homogenized coefficients in the case of discrete linear 
elliptic equations, whereas [27, 29] provide with a complete numerical analysis (convergence rate, 
quantitative central limit theorem) of this method when the coefficients are independent and identi¬ 
cally distributed. In the case of the unbounded integral functionals considered here, it is not clear 
whether the approximation by periodization in space converges in general. Note that for problems 
with discrete stationarity (like the random chessboard or periodic structures) such periodizations are 
not needed but the sequences e, R have to be discrete: = l/k,R = k with A: € N. 

Before we state our results, let us shortly discuss periodization in law on the example of random 
inclusions. Let 


V{y,A,u}) := a(A)liRd\£;(^)(?/)-h6(A)lE(^)(y), (2.15) 

where a, b : —)■ [0, oo] are convex functions and E{uj) is a stationary random set (to be thought 

of as a collection of inclusions). A typical choice for E is E = UnLi-S(9n)) where p = {qn)n is 
a stationary point process in For R > 0, set = {qn,R)n, := B{qn,R), and define 
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V^{y, A,uj) := a{A)l^d\Eii{uj)iy) + K^)l_E«(aj)(y)- For to be an admissible periodization in law 
of V, it is sufficient that {p^)r is an admissible periodization of p in the following sense: 

(i’) Periodicity in law: for all i? > 0, is Q/j-periodic and r^-stationary with respect to transla¬ 
tions on the torus T/j = 

(ii’) Stabilization: for all 9 £ (0,1) and for almost all oj there is Rg{u}) > 0 so that for all R > Rq{uj), 
p^{u)) n Qor = p{oj) n Qbr. 

Let us now describe three typical examples, for which the periodization in law can be explicitly 
constructed: 

(a) Poisson point process. If p is a Poisson point process on a suitable periodization in law is 

given by p^ = + pH Qr). Property (i’) is satisfied by complete independence of the 

Poisson process, while property (ii’) holds for 9=1. 

(b) Random parking point process. We now let p be the random parking measure on M'’* defined in [38]. 

Penrose’s graphical construction is as follows: p is obtained as a transformation p := T'(po) of a 
Poisson point process po on x R+. Following [30, Remark 5], we may periodize the underlying 
Poisson process on Qr as above by setting p^ := + Po H {Qr x define 

a periodization of p as p^ := T{pq). Property (i’) holds by construction. Property (ii’) is more 
subtle and relies on the stabilization properties of the random parking measure. By a union 
bound argument we may rephrase the exponential stabilization of [39, Lemma 3.5] as follows: 
There exist K,k > 0 and for all R > 0 there exists a random variable tR such that 

^ P Ll Qr does not depend on po H ((M'’* \ Qt^) X R'*'), and in particular, for any locally finite 
subset F C (R"^ \ Qt^) x R"*", we have 

Qr n T{F U (po n (Qi^ X R+))) = Qij n r(po n {Qt^ x R+)); 

- for all L > 0, F[tR > L] < 

Hence, for all 9 G (0,1), 

CX) 

P[p n Qbr does depend on po H ((R'’* \ Qr) x R'*’)] 

R=1 

OO CO 

< P[t0ij >R]< < OO. 

R=1 R=1 

By the Borel-Cantelli lemma, pCiQbr does eventually not depend on pon((R'’*\(5i{) xR+) as ii f oo 
almost surely. In particular, for all 9 G (0,1) and for almost all iv, there exists some R^^ > 0 with 
the following property: for all R > R^, p{uj) n Qbr does not depend on po(w) n ((R'^ \ Qr) x R^), 
so that p^{oj) n Qbr = p{oj) n Qbr, that is, (ii’). 

(c) Hardcore Poisson point process. As noticed in [21, Step 1 of Section 4.2] (in a more general 
context), a hardcore Poisson point process p on R'^ can be obtained by applying the transformation 
T of (b) to a Poisson point process po on R'’^ x (0,1), that is, p = T{pq). The same argument as 
in (b) (with however easier bounds) then yields the desired periodization in law of p. 

In terms of applications, example (b) is of particular interest since it is used in [30, 26, 20] as a 
suitable random point set for the derivation of nonlinear elasticity from polymer physics started in [4]. 

We now turn to our T-convergence results with periodic boundary data under periodization in law 
(yielding in particular an approximation of the homogenized energy density of the form (2.14)). The 
main difficulty is the following: in order to carry out the analysis as before, we would somehow need 
the uniform sublinearity of the correctors associated to the family of periodized integrands (R'^)h>o- 
It is not clear to us whether this uniformity holds in general. In what follows, we give an alternative 
argument in the scalar case with fix domain, as well as in the very particular example of well-separated 
random stiff inclusions in a soft matrix. Note that because of the separation assumption this does 
not include the example (a) of Poisson inclusions; however we believe that up to some technicalities 
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the proof could be adapted to that case. For simplicity, we consider spherical inclusions, bnt random 
shapes could of course be considered too. 


Corollary 2.13 (Approximation by periodization). Let V, and he as in Theorems 2.2 

and 2.6 for some p > 1. Also assume that one of the following holds: 

(1) m = 1, and domC(y, -, 0 :) = domM for almost all y, u; 

(2) p > d, and V{y,A,uj) < C{\ + |A|p) for all oj and all y ^ , for some random stationary set 

‘T satisfying almost surely, for all n, and some constant C > Q, 


ini 

ivz m^m^n U 


Rn<C. 


Let (C^)_R>o he an admissible periodization in law for V in the sense of Definition 2.11, and for all 
e > 0 and all A £ denote by J|’®^(-,A, •) the following random integral functional on the unit 

cube Q = 

jr(u,A,a;):= [ vR%y/e,A + Vu{y),u;)dy, uGW^i{Q;Rn, 

Jq 

where Wj^er(Q;K™) denotes the closure in M™') of the set of smooth periodic functions on Q. 

Then, for almost all lo and for all A, the integral functionals ®^(-, A, w) T-converge to J(- + A-x,a;) 
on the space lFper((5;K"*). In particular, we have the following approximation ofV by periodization: 
for all A G for almost all lo, 

y(A)=limlim inf -r V^{y,tA Vu{y),Lo)dy 

tn 

= lim lim E inf -r V^{y,tA+ Vu{y),-)dy . 
tn Rtoo JQr 

By convexity, if A £ intdomC, then the limits t f 1 can be dropped. Note that domC = holds 

in case (2). The same results also hold for V, J^, J replaced by W, I, provided that p > d. □ 


3. Proof of the results for convex integrands 

This section is dedicated to the proofs of Theorem 2.2, Corollary 2.4, and Corollaries 2.10 and 2.13. 
Let P be a convex r-stationary normal random integrand. Up to the addition of a constant, we may 
restrict to the following stronger version of (2.4): for almost all lo, y, we have, for all A, 

^\A\P < V[y,A,uo), (3.1) 

for some C > 0 and 1 < p < oo. We assume that 0 belongs to the interior of the domain of the convex 
fnnction M := snpess^^ ^^ Viv, 

Following the strategy of [36, Theorem 1.5], we proceed by truncation of V. We let {V^)k be an 
increasing sequence of r-stationary convex normal random integrands x x 11 —)• M'*' 

such that, for almost all lo, y, and for all A, we have 

lim U^(y,A,a;) = U(y,A,a;), and hA\P < V^{y,A,uj) < C'^{\A\p + 1), (3.2) 

fctoo O 

for some C > 0 and some sequence f oo (see [36, Lemma 3.4] for such a construction). Let 
Ho U ^ be a subset of maximal probability on which all these assnmptions (about V and the U^’s) 
are simultaneously pointwise satisfied. 

We shall prove the existence of a subset 11' C Hq of maximal probability such that, for all a; G H' 
and all bounded Lipschitz domains O C the functionals J^[-,lo',0) F-converge to the fnnctional 
J{-',0) on VF^’^(0;M™), where we recall the definitions 

Je{u,u}-,0) := [ V{y/e,Vu{y),u})dy, J{u;0)-.= [ V{Vu{y))dy. 

Jo Jo 




14 


M. DUERINCKX AND A. GLORIA 


As usual, the proof of F-convergence splits into two parts: the proof of a lower bound (F-liminf 
inequality) and the explicit construction of a recovery sequence which achieves the lower bound (F- 
limsup inequality). 

3.1. Preliminaries. We first need to briefly recall the standard stationary differential calculus in 
probability (first introduced by [37, Section 2]), as well as some results on ergodic Weyl decomposi¬ 
tions. 

3.1.1. Stationary differential calculus in probability. Let 1 < p < oo. For all 1 < f < d, consider 

the partial action (T^)/igR of (M, -|-) on L^(fl), defined by {T^f){oj) = f{T-heiOj), for h G M. The 
actions {TffjheM. (for 1 < i < d) commute with each other and are unitary and strongly continuous by 
Lemma A.5. For all i, we may then consider the infinitesimal generator Di of defined by 

A/=lim^^, /eL^>(L!), 

whenever the limit exists in the strong sense of L^(n). By classical semigroup theory, the generators Di 
are closed linear operators with dense domains Pj C L^(n), and the intersection VF^’^(fl) := nf=i 
is also dense in L^(n). Moreover, is endowed with a natural Banach space structure. 

For / G LF^’^(fl), its stationary gradient is then defined by Df := (Pi/,..., P^^/) G L^(n;M'^). 
Through the usual correspondence between random variables and r-stationary random fields as re¬ 
called in Appendix A.2.1 (for all g G L^(n), write g{x,uj) := g{T-xUj), defining g G L[(^^(M'^;L^(n))), 
we may define Df[x,u) := Df{T-xOj) for all x. By unitarity of the action T, the operator P is 
skew-symmetric, so that the following “integration by parts formula” holds, for all / G and 

5 G Wi’^’'(L!),p'=p/(p-l), 

E[P/] = 0 and nfDg] = -E[gP/]. 

As explained in Section A.2.2 (see in particular Lemma A. 7), for almost all w, the function P/(-,a;) 
is nothing but the distributional derivative of /(•,tu) G L|(j^(M'^), and the following identity holds: 

= {/ G iyi'/f(M'';L^’(L!)) : f{x + y,u) = f{x,T.yUj)ffx,yM- (3-3) 

This justifies that in the sequel we simply use the notation Df = V/. 

3.1.2. Ergodic Weyl decomposition. Ergodicity of the measurable action r of (M'^,-|-) on the probability 
space (n,P, P) is crucial in the sequel. Let 1 < p < oo. In analogy with the classical Weyl subspaces 

we define the subspaces of potential and solenoidal random fields with respect to the 
differential calculus associated with the group action in the following way: for p' = p/{p — 1), 

Lpot(^^) = {/ e LP(0;M'') : E[f • (V x p)] = 0, Vp G M'^)}, (3.4) 

L^„/L!) = {/ G LP(fl;M'') : E[f ■ Vg] = 0, V 5 G W'’P'(L!)}. 

Reinterpreting these definitions in physical space, we easily obtain the following reformulations in 
terms of stationary extensions: 

Lpot(^) — {f ^ L^(fl;M'^) : for almost all cj, a; i-A /{txOj) G L[(j^(M'^;M'^) is potential}, (3-5) 

Lg^l(n) = {/ G L^(fl;M'^) : for almost all cj, a; i-A /{txOj) G L[(j^(M'^; M'^) is solenoidal}, 

where a function h G L[(j^(M'^) is said to be potential (resp. solenoidal) if V x h = 0 (resp. V • / = 0) 
in the distributional sense. Constant functions belong to both subspaces, and we further define 

= {/ G : nn = 0} and F^JQ) = {/ G : E[/] = 0}. (3.6) 

The spaces LpQ^(n), L^^/n), and P^/fl) are all closed in L^(fl;M'^), and the following Banach 

direct sum decomposition holds 


LP{n-R'^) = FP^ffn)(BFPjn)(BR'^ 


( 3 . 7 ) 
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as well as the following density results: 

= adhLP(n;R.){V5 : g S (3-8) 

F^jn) = adhLP(0;M.){V xg: ge 

Since we did not find a suitable reference for these results (besides the Hilbert setting), we provide a 
proof of (3.7) and (3.8) in Appendix A.3. 


3.2. F-convergence of truncated energies. Since the approximations of V all satisfy standard 
polynomial growth conditions, we can appeal to the classical stochastic homogenization result of [19] 
(which could be reproved via a direct adaptation of the (periodic) arguments of [36, Theorem 1.3]). 
More precisely, there exists a subset Hi C Hq of maximal probability such that, for all w G Hi, all 

k, and all A G the following limit exists (as a consequence of the Ackoglu-Krengel subadditive 

_^ 

ergodic theorem) and defines the homogenized integrand V : 

H^(A) = lim inf -Z V’"{y,A+ V(p{y),uj)dy, (3.9) 


where Qr := [—)'^. By dominated convergence, this convergence also holds when taking the 
expectation of the infimum. In addition, for any bounded Lipschitz domain O C and for all 
cu G Hi and all k, the functionals F-converge, as e 0, to the functional defined 

by 

4{u,u;-,0) := [ V>^{y/e,Vu{y),u;)dy, and j\u-0):=[v\Vu{y))dy. 

Jo Jo 

Since k e->■ is increasing, k eA is increasing as well, and for all A G we may define 

17(A) ;= limfc_).oo H (A). In particular. 


V{A) = supH''(A) 

k 


sup lim inf 



V>^iy,A + VcPiy),uj)dy. 


(3.10) 


It remains to pass to the limit fe t oo in the F-convergence result. The key is to prove the commutation 

of homogenization and truncation, which we do in Subsection 3.4 below. 

_^ 

Alternative formulas for V are obtained in Lemma 3.7. Since V is convex and everywhere finite, 
it is continuous on and we may directly deduce from the definition 17(A) := sup^. 17 (A): 


Lemma 3.1. The map V : R™^'^ —)• [0,oo] is convex and lower semicontinuous. 


□ 


3.3. F-liminf inequality. In view of the definition of 17, the F-liminf inequality is an elementary 
consequence of the monotone convergence theorem: 

Proposition 3.2 (F-liminf inequality). For all a; G Hi, all hounded domains O C R'^, and all 
sequences {u^)^ C II7^’^’(0; R™) with ^ u in H7^’^’(0; R™'), we have 

lim inf O) > J{u; O). 


Proof. Let O, {us)e, and u be as in the statement. Then, for all w G Hi and all A: G N, using the 
F-liminf result for towards (see Section 3.2), and that 17 > 17^, 

Ymiuii Ji;{us,iw]0) > liminf Jg (u£,a;; O) > J^{u]0), 

eiO e4.0 

so that, by monotone convergence, 

lim inf O) > lim J^(tt; O) = J(u; O). □ 

£4-0 fctoo 


From this F-liminf result, we deduce the locality of recovery sequences, if they exist. 
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Corollary 3.3 (Locality of recovery sequences). If for some uj € rii, some bounded domain O C 
and some function u € M™'), there exists a sequence {ue)e C with ^ u 

in and J^{u£,L 0 ]O) —>• J{u;0), then we also have j£{u£,io;0') —?■ J{u]0') for any 

subdomain O' C O. Hence, by an extension result, the L-limsup inequality on a bounded Lipschitz 
domain O implies the L-limsup inequality on any Lipschitz subdomain O' d O. □ 

Proof. Choose a subdomain O' C O, and define O" := 0\ O'. We then have by assumption 

J{U]0) = lim J{Ue]0) = hm [j^{Ue,i 0 ;O') + Je{Ue,UJ]0")) 

> lim inf O') + lim inf J^{ue,U}; O"). 

etO ei-O 

Now by Proposition 3.2 we have liminf£4,o OO ^ J{u;0') and liminfg^o^ 

J(tt; O"). The conclusion then follows from the identity J(u; O') + J{u] O") = J{u] O). □ 


3.4. Commutation of truncation and homogenization. The crucial ingredient to prove the 
commutation of truncation and homogenization is the reformulation of the asymptotic homogenization 
formula in the probability space. For that purpose, we first introduce the following proxy for V: 

P(A):= inf E[C(0, A+ /,-)]. (3.11) 

Likewise, for all A: € N, we set 

P"(A):= inf £[^"^(0, A+ /,.)]. (3.12) 


In this case, due to the growth condition (3.2), we may prove (see Lemma 3.7 below) that 


limE 

£ 4.0 


inf -/ V^{y,A + V(i>{y),-)dy 
0GWo’P(O/£;R”i) Jo/e 


= P^{A). 


However, for V itself, this equality has no chance to be true if A G domP \ domM since the left-hand 
side could be inhnite (because of the Dirichlet boundary condition) while the right-hand side is not 
— see Example 2.3. We thus rather use a “relaxed version” of the Dirichlet boundary conditions and 
set for all e > 0 


P^{A,u}]0) ■.= mf -/ V{y,A + A/(t){y),u)dy. (3.13) 

,^ewLP(0/e;R^) jQ/e 

fo/e^<P=0 

As opposed to the case of Dirichlet boundary conditions, there is no natural subadditive property in 
this definition (two test-functions on disjoint domains cannot be glued together). This difficulty will 
be overcome by using a more sophisticated gluing argument that relies quantitatively on the following 
sublinearity property of the correctors. 

Lemma 3.4 (Sublinearity of correctors). For all A G there exists a corrector field G 

Mes(D;VL/jf(M'^;]R"^)) such that V(^a(0, •) G F^„t(D)™, and 

P(A) =E[D(0, A+ V(^a(0, •),•)] • 

In addition, </Ja is sublinear at infinity in the sense that, for almost all ui € Ll, 

eipA{-/£,uj) ^ 0 (3.14) 

weakly in VL^’P(0;M”^) for all bounded domains O C M'^. □ 

Remark 3.5. Although the space {Vg : g G ly^’P(D)} is dense in PpQ^(D) (see (3.8)), the inh- 
mum (3.11) dehning P{A) cannot be replaced in general by an infimum over this smaller dense sub¬ 
space because of a possible lack of strong continuity of the functional, see however Proposition 5.1. □ 
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Proof. Let A S be fixed. By convexity and by the lower bound (3.1) onV ,x'~^ ^[1/(0, A + x, •)] 

is lower semicontinuous and coercive on FpQ^(n)™', and therefore attains its infimum. Let g £ 
be a minimizer. The r-stationary extension {x,uj) i—>• g{TxUj) of (7 is a potential field on for almost 
every ui. Hence, there exists a map (fA £ Mes(H; M"*)) such that g{TxOj) = S/(pa{x,u)) for 

almost all x^oj (see indeed Proposition A. 10). The claim now follows from the combination of the 
following two applications of the Birkhoff-Khinchin ergodic theorem: for almost all cj, 


VipA{-/£,uj) —^ 0, (weakly) in L^(0;M”^), (3.15) 

e f <^A(y/e,w)dy 0. (3.16) 

Jo 

Indeed, by Poincare’s inequality and (3.15), the sequence y i-A eif\{y/£,uj) — £ (/ 9 A(z/e, a;)d^; con¬ 
verges weakly to zero in for almost every uj. Combined with (3.16), this implies (3.14). 

To conclude, we turn to the proofs of (3.15) and (3.16). The weak convergence (3.15) to zero is a di¬ 
rect consequence of the Birkhoff-Khinchin ergodic theorem in the form V(/?a(‘/£j w) —^ E[V(/?a( 0, •)] = 
0 in LP(0;M”*). It remains to prove (3.16). We may assume wlog ¥?a( 0, •) = 0 almost surely, so that 


e 



ipA{y/£,uj)dy 


£f (^a(-,w) 
Jo/e 

= 

Jo/ 

/ X ■'V(pA{tx,uj)dtdx 

£ Jo 



rl/e 

r 



< 

i 

'T X ■ V(pAitx,uj)dx 

Jo 

dt. 


(3.17) 


For almost all a;, the function fjcoit) := ^qX ■ VipA{tx,uj)dx is continuous on (0,oo). By (3.15), 
—>• 0 as 11 00 for almost all oj. By joint measurability and (local) integrability of Vi^^Aj and by 
stationarity, 0 is a Lebesgue point of V(/9A('ja;) for almost all uj, and hence limsup^^Q |'0a;(i)| < 00 for 
almost all uj. The result (3.16) then follows from (3.17). □ 


For all A £ domP, let tpA be defined as in Lemma 3.4, and let Ha C Hi be a subset of maximal 
probability such that (3.14) holds on Ha for all bounded Lipschitz domains. Restricting Ha further, 
the Birkhoff-Khinchin ergodic ensures that, for all uj £ Ha, we have for all bounded subsets O C M'’* 
and all t £ Q, 


and 



\/iptA{-,Uj) 



(3.18) 



tA ViptA{y,^),<^)dy 



E[K(0,tA -h V99tA(0, •))] 


P{tA). 


(3.19) 


We now turn to the proof that lime O) = P{A) for all A for almost all a; £ H. The following 

inequality is the most subtle part. 


Lemma 3.6. For all A £ intdomP and all hounded domains O C there exists a sequence 
'4’A,o,e £ Mes(H; iy^’P(0/e; M™)) such that, for all uj £ Ha, V'0A,o,e(', <^) = 0, 

weakly in 1 F^’^( 0 ; M”*) as e 0 , and 


P(A) > limsup 


elO Jo/e 


/ V{y, 


A -|- V'ijjA,o,e{y-,^)-,^)dy > hmsupP'^(A,a;; O). 


etO 


(3.20) 

□ 


Proof. Let A £ intdomP be fixed, and let uj £ Ha- For all t £ [0,1) H Q and e > 0, set 


A 


UJ 

0,£,t 


= —t 



and ilJA,o,e,t{x, uj) := tipA/t{x, uj) + Aq^^ ^ ■ x. 
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By definition, = 0, and by Lemma 3.4, ^ 0 in as e 0. 


Hence, 




(A, ui;0) < 'f V{y,A + VipA^^sAv ^=: O). 

JOje 


By convexity 

Pl{A,u-,0) = / V{y,AptSIipAit{v,^)P ^o,ep^Ay 

Jo/e 


< 


tj'o/ + V(/?A/t(y,w),a;)dy + (1H ^y, j, dy. 


Since 0 G intdomM, there exists d > 0 such that adhH^ C intdomM. As t is rational and u G Oaj 
we have ^ 0 as e 4- 0 by the Birkhoff-Khinchin ergodic theorem in the form of (3.18). Hence 

there exists ^a o t ^ ^ such that, for all 0 < e < £)!( q j, we have 


_ A 

l-t 0 ,e,t 


<s, 


and therefore. 


This implies that 


O/e 


^ [y^ ) dy < sup M(A') < oo. 

V / |A'|<5 


hmsuplimsupP/(A,a;; O) < limsuplimsup-r V{y,A/t + 'S/<fAit{y,ijj),ijj)dy. 
itl,t6Q £4.0 Hl,iGQ £4.0 Jo/e 

By the Birkhoff-Khinchin ergodic theorem in the form of (3.19) and the continuity of P in the interior 
of its domain (as a consequence of convexity), this yields 

limsup limsup PA^i ^'i O) < hmsupE[K(0, A/t + V(/?A/i(0, •))')] = limsup P{A/t) = P{A). 

£4-0 H14GQ H14GQ 

We have thus proved: 

limsup limsup ( (pf{A,uj;0) - P{A)) + ||eV’A,o,£,t(Ve>‘^)llLP(0;M-) ) = 0- 

Itl.ieQ £i0 ' J 

By Attouch’s diagonalization lemma (see [8, Corollary 1.16]), this implies the existence of a sequence 
(AA,o,e)e with V'a,o,£ G {O/e-,MP) such that VV'a,o,£ = 0, limsup^ Pf (A,a;; O) < P(A), and 

ei/)A,o,£('/^i‘^) —^ 0 in L^(0;M™), for all to G Oa- By the choice of A, P{A) < oo, so that the lower 
bound (3.1) on V implies that the sequence (V'0 a,o,£('/^i‘^))£ bounded in L^(0;M”*). We thus 
conclude that eil^A,o,e {'^ weakly in LL^’^(0;M™), as claimed. □ 

In the case of standard growth conditions (thus e.g. for the K^’s), the corresponding inequality 

(3.20) in Lemma 3.6 is indeed an equality. The following lemma gives equivalent definitions for the 

_ 

V ’s, which will be crucial in the sequel. 


Lemma 3.7. Let O be a bounded Lipschitz domain ofW^. For all oo G fli, all k, and all A G 
the following quantities are well-defined: 


vliA) 

F2(A,to) 

^3'(A) 


lim inf -f V’^{y,A+ V(j){y),oo)dy, 

efO 0ew4’^(O/£;IR^) Jo/e 
lim inf 

£4.0 ,#)6Wl'P(0/e;K"*) 

So/e 

_ inf E[K^(0,A+ /,.)], 



(3.21) 

(3.22) 


(3.23) 
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and we have 


(/‘(A) = 1/;(A) = ^‘(A) = ^^(A). 


□ 


lim inf 
\0\ efO , 


This result is standard (see for instance [31, Chapter 15]) and we display its proof for completeness. 
Note that the formulas (3.22) and (3.23) for V will be shown to be equivalent to V, whereas formula 
(3.21) is in general larger than V. 

Proof. Let O C be a bounded Lipschitz domain, and k £ N. By the dehnition of T-convergence 
for on LLq’^(O; M™') and the convergence of infima with Dirichlet boundary conditions, for all A 
and w € we have 

F^A)=^ inf [v\A + Vf^) 

/ V\y/e, A + ■)dy = v\{A). (3.24) 

)Jo 

Likewise, the L-convergence result holds on {u G W^’^{0) : JqVu = 0} so that the identity 

v\K)=vI{K) 

also follows from the convergence of infima. Since Lemma 3.6 (applied to instead of V) yields 

1 ^2(■^) — it remains to prove that V^{A) < l^i(A) for all A. 

Let O' C be a bounded domain. By the coercivity and the lower semicontinuity of the integral 
functional (which follow from the growth condition (3.2) and the convexity of V^), there exists a 
minimizer f G L°°(n; VLq’^(O'; M™')) (where measurability follows from Proposition A.15) such that, 
for almost all u, 

/ C^(?/,A +VC(y,a;),w)dy = inf -f V’"{y,A + V(j){y),uj)dy. 

JO' (/.GVEo Jo' 


Set 


I7TT / C{x + z,TzU!)dz = f C{x + z,TzUj)dz. 
\ iR-^ J-x+O' 


Clearly, ^ is well-defined and stationary, belongs to M”^), and 


V^(x,a;) = 


Hence 


TTk, 


—X+O' 


= E 


VC(® + z, TzUj)dz. 


(^0,A + £vCiz,T,-)dz, 


C3(A) <E[y'=(0,A + ve(0,-),-) 

and by convexity of 

F5(A)<e[/ vHo,A + vaz,T,-),-)dz 


U O' 


By stationarity and the Fubini theorem, we may conclude 


vl{A)<i E[H"(z,A + VC(z,-),-)]d^ = IE inf / V\y,A + Vcf{y),-)dy 
JO' L(/>GWp’^(0';R"*) Jo' 

With O' := O/e, the claim H 3 (A) < ^^(A) follows by the dominated convergence theorem and the 
growth condition from above (3.2). □ 

The following result proves the equivalence between formulas (i), (ii) and (hi) in Theorem 2.2. 
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Proposition 3.8 (Commutation of limits). For all bounded Lipschitz domains O C and all 
A G we have for almost all oo 

V{A) = P{A) = limlimP^(tA, w; O). 

tti eJ-O 

By convexity, for all A ^ cldomP, this takes the form P(A) = P(A) = lim^j^oO)- D 

Remark 3.9. Although not stated explicitly, this result proves the commutation of truncation and 
homogenization. By monotone convergence (cf. the proof of V = P below) we have for all e > 0 and 
almost every uj, 

inf / P(y,A +V(/>(?/),a;)(iy = sup inf / V^{y, A + V(t){y),uj)dy 

0ewi.P(O/e;iR’") Jo/e k </>eWp’^’(0/e;iR’") Jo/e 

/o/^v^=o 

so that Proposition 3.8, combined with (3.22) in Lemma 3.7, yields the desired commutation result 


limsup inf -r V^{y,A + V(p{y),uj)dy 
eto k <Aewi.P(o/e;K™') Jo/e 

fo/e^'f>=0 

= P(A) = supP^(A) = suplim inf -f V’^{y,A + V4>{y),uj)dy. 

k k efO ^ewi,P(0/s-,R^) Jo/e 

fo/e 

□ 


Proof of Proposition 3.8. We split the proof into two steps. 

Step 1. Proof of P = P. 

Let A G domP. By (3.23) in Lemma 3.7, for all k, 

v\A)= inf E[P"(0,A + /)]. 

By convexity, / i—)• E[P^(0, A + /)] is lower semicontinuous on PpQj.(n)™', and by coercivity, the 
infimum is attained. Hence there exists gk G such that 

F'^(A) =E[H'=(0,A + 5fc)]. 

By the uniform growth condition from below (3.2), {gk)k is bounded in L^(n; 

^2-P+^E[\gkn - ^\A\P < 1e[|A + gk^ < E[V\0, A + gk)] = v\a) < V{A). 

Let g G PpQ^(n)™ be a cluster point of (gk) for the weak convergence of L^(f2; M™'^'^). We have along 
the subsequence 

F(A) =supF^(A) = limE[H^(0,A + 5fc)]. 

k 

Since k i—)• is increasing and / i—)■ E[H^(0, A + /)] is lower semicontinuous for the weak convergence 

of L^(f2; this yields for all i 

V{A) >liminfE[H^(0,A + 5 fc)] > E[i/^(0, A + c/)]. 

fctoo 

We then conclude by monotone convergence that 

F(A) >E[H(0,A + 5 )] > inf E[H(0,A + /)] = P(A). 

For A ^ domP, the above inequality is trivial so that V{A) > P{A) for all A G For the 

converse inequality, note that for all A, 

P(A)>sup inf E[i/*^(0,A + /)] = supF''(A) =F(A). 
k k 
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Hence, V = P, as claimed. 

Step 2 . Proof of limt|i limej^o P^{tA, uj;0) = V (A). 

Since for A G domH and t G [0,1), tA G intdomH, Lemma 3.6 and Step 1 yield for almost all uj 

V{tA) = P{tA) > limsup P^(tA, cj; O). (3.25) 

e4-0 

By (3.22) in Lemma 3.7, for all A G and almost all oj, 

liminf P^(A,a;; O) > suplim inf H^(y, A + V(/>(y), a;)dy = sup H^(A) = H(A). (3.26) 

£4,0 k efO </,evrLP(o/e;K"*) Jo/e k 

fo/e 

Combined with (3.25), this yields lim^ P''(tA, ca; O) = V{tA) for almost all a;, for all A G domH and 
t G [0,1). By convexity and lower semicontinuity of V, this implies for all A G domH 

limlimP'^(tA,a;; O) = limH(tA) = H(A), (3.27) 

ttl efO ttl 

and (3.26) ensures that this equality also holds for A ^ domH. By convexity and by (3.27), the 
function A i—)• lim^ P^(A, to; O) is continuous outside 9domH, so that the limit t f 1 can be omitted 
for A ^ ddomV. □ 

3.5. Proof of Theorem 2.2: L-convergence with Neumann boundary data. It only remains 
to prove the P-limsup inequality. 

Proposition 3.10 (L-limsup inequality with Neumann boundary data). Assume p > d. There 
exists a subset C Hi of maximal probability with the following property: for all to G H', all bounded 
Lipschitz domains O C and all u G IL^’^(0; M™), there exists a sequence {ue)e C VL^’^(0;M"*) 
such that Ue ^ u in and J^{ue,uj;0) ^ J{u]0). □ 

Proof. We split the proof into three steps. We first treat the case of affine functions, then the case 
of continuous piecewise affine functions, and finally the general case. The novelty of our approach is 
the careful gluing argument needed to pass from affine to piecewise affine functions. 

Step 1. Recovery sequence for affine functions. 

In this step, we consider the case when n = A • x is an affine function. More precisely, we prove the 
existence of a subset H' C Hi of maximal probability with the following property: given a bounded 
Lipschitz domain O C M'^, for all a; G H' and all A G intdomP, there exists a sequence {uf^^)e C 
M™') with ^ A-x weakly in IP^’^(0; M™') such that, for all Lipschitz subdomains O' C O, 
we have Je{uf^^,u}', O') —)• J(A • x; O'). By Corollary 3.3, it suffices to prove this for O' = O. 

By Lemma 3.4 and Proposition 3.8, there exists a sequence (pA £ Mes(H; iyjQ^(M'^; K™')) such that, 
for all UJ G Ha, we have e(/9A(-/e, ta) —^ 0 weakly in IP^’^(0; M™) and, by the Birkhoff-Khinchin ergodic 
theorem in the form of (3.19), 

V (A) = P(A) = lim / V(y,A + VpAiv, uj),uj)dy. 

efO Jo/e 

In particular, by a change of variables, this yields 

J{A ■ x;0) = |0|P(A) = lim J^{A ■ x + e(pA{-/e, ca), w; O). 

£4,0 

The function Ue’'^{x) := A • x + eipA{x/e,uj) thus satisfies u^’‘^ A • x in IP^’^(0;M™) and 
Je{ut'^,uj] O) —)• J(A • x; O) as e 4 , 0, for all u G Ha. 

We then define H' C Hi as the (countable) intersection of all Ha’s with A G n intdomP, 

which is still of maximal probability. Let A G int domP and a; G H' be fixed. Choose a sequence 
(Ari)n C n int domP such that A„ —)• A. For all n, we have already constructed a sequence 
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{'^e,n)e C such that ^ A„ • x in hF^’^(0;M™') and Ji.{u^,^,ui]0) J(A„ • x;0). 

Since by convexity, V is continuous on int domi^, we have 

limsuplimsup {\Je{u‘^^^,ui-,0) - J(A • x;0)\ + - A • x||lp( 0 ;K-)) 

nfoo etO 

= limsup (|J(A„ ■X]0) - J(A • x;0)| + ||A„ • x - A • x||lp( 0 ;K-)) 

ntoo 

< limsup (|0||F(A„) - y(A)| + Co|A„ - A|) = 0. 

ntoo 

By the Attouch diagonalization lemma (see [8, Corollary 1.16]), this implies the existence of a sequence 
{v‘^)e such that O) J(A • x; O) and ^ A • x in L^(0;M™') for all u £ O'. By the p-th 

order lower bound for V, we conclude that converges weakly to A • x in VC^’^(0; M™). 

Step 2. Recovery sequence for continuous piecewise affine functions. 

Let tx £ O', O C be a bounded Lipschitz domain, and u be a continuous piecewise affine 
function on O such that Vu £ int domR pointwise. We shall prove that there exists a sequence 
{u‘^)e C with ^ u weakly in M"*), such that Je{u‘^,uj;0) —)• J{u;0). For 

that purpose, the major issue consists in gluing the recovery sequences for the different affine parts 
together, which requires a particularly careful treatment. 

Consider the open partition O = ]+)f=i O; associated with u (note that the 0;’s have piecewise 
flat boundary outside dO), and define Ci + Ai ■ x := u\oi, with A/ £ intdomR, for all 1 < I < k. 
Let Ai := (Uf=i'90;) \ be the interior boundary of the partition of O, and for all r > 0 set 
Air ■= {Ai + Br) n O = {x £ O : dist(x. Ad) < r}, the r-neighborhood of this interior boundary. By 
Proposition A. 17, for all 0 < k < 1 and r > 0, there exists a continuous piecewise affine function 
on O with the following properties: 

(i) Vxre,r = Vx pointwise on O \ Air, and 

limsup sup ||xK,r — i.i||L°°(o) = 0; (3.28) 

rto 0<K<1 

(ii) A/u^^r £ conv({A; : 1 < Z < k}) (s intdomC pointwise (where conv(-) denotes the convex hull); 

(hi) denoting by O := r open partition associated with u^^r, and denoting ^+A(,. ^-x := 

XK,r|ol, ^ ^ then, for any i,j with 90^,r bl dOi^r / 0, we have |AJ.,, — Ai,^r\ < k. 

We shall approximate u with these refined continuous piecewise affine functions u^^r having smoother 
variations; in the sequel, we shall successively take the limits «; 0 and r j, 0. 

Since tx £ fl' and O C are hxed in the argument, we drop them from our notation. Fix 
K,r >0. By Step 1, for all 1 < f < nK,r there exists a sequence (x* f, ,,)^ C LFjq^(M'^; M™) with 
^e,K,r ^K,r + ^K,r ' ^ ™ VFj^’^(M'^; R™') and such that, for all Lipschitz subdomains O' C O, we have 
j£(x* ,,,,, cx; O') —)• d(A(.,, • x; O'). For all p > 0 and all 1 < f < n^^r, define the sets 

0'^r,ri :={3: G O : dist(x, < r]} = O n {O^^^ + R^), 

Olrr,rj :={a^ G 0^,r ^ dist(x, 90^^^) > rj} 

Let then Yl'i=i XK,r,r] = 1 be a partition of unity on O, where, for all 1 < f < riK^r, the smooth cut-off 
function X^nr-q has values in [0,1], equals 1 on O'Arq and vanishes outside 0)Aj.^, and satisfies the 
bound \ A'xK,r,q\ — pointwise for some constant O' > 0. We now set 

Ue,K,,r,q •= UK,r + ~ {^K,r A A^ ,, • x)) XK,r,q- 

i=l 
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By the Sobolev compact embedding for p > d, we have ^ ^ cj. ^ ^ • x in L°°(0; M"*) as e 0, 

and hence limsup^ limsup^ ||tte,K,r,r; — 'WK,r||L°°( 0 ) = 0; so that (3.28) yields 

limlimsup limsup limsuplimsup \\tUe,K,r,ri — ^^||l°°(o) = 0- (3.29) 

rtO fttO »?4-0 e4-0 

Let us now evaluate the integral functional O) at tUe^K,,r,r) for t G [0,1). Since 

t'^Ue,K.,r,ri = ^XK.,r,r)^'^e,K,r 


2=1 


t - 

+ ^ • x))Vx\,r,^ + (Vu^,r “ K,r)xl,r,v) > 

2=1 

and (1 — t) + r -77 = 1) we have by convexity and non-negativity of V 

< (1 — -|~ t / XK,r,r]{y)^ {y/^'}'^'^e,K,r{y)^ 

/d*+ 

2=1 '~^K.,r,r} 


— (1 't)Ee,K.,r,rj,t '^e('^e,K,r5 


(3.30) 


2=1 


where the error term reads 


Ee,K,r,ri,t = j ^ («A,r-(y) “ + ^Jc,r ' y))VXK,r,r,(y) 


+ (Vu«,,.(y) - ,.,,(?/) j,wjdy. 

For all i, set A^K,r ,?7 ■— {j ■ J 7 ^ L Ol^r,ri H We then rewrite the argument of the energy 

density in the error term as 


Se,K,r,ri{y) ■ — 


'i'K,r 

Y1 {('^l,KAy) - (A,r + K,r • y)WXK,rAy'^ + (Vw^,r(2/) “ K,r)AK,r,n 


< 


2 = 1 
c 


\K,K,r - (4,r + K,r ' a;)||L°°(0) + SUp SUp |A^_^ - 


^ tl 


Since by definition limsup^^g ~ Ai,r\ — ^ fo'^ L we have 

limsuplimsuplimsup5£^K^r,r7(y) = 0 

fttO ?7t0 etO 

for all r,p > 0. By assumption, there exists (i > 0 such that adhB^ C intdomM. Hence, for all 
r,t > 0 there exists > 0 such that for all 0 < At < /t^y there exists ry^.r > 0 such that for all 
0 < rj < there exists En^r,r],t > 0 with the following property: for all 0 < e < £K,r,ri,t^ we have 


t 


1 - t 


:S, 


e,K.,r,ri 


< 6. 


L°°(0) 


This yields the bound 


E£,K,r,v,t < |0| sup M(A') < oo, 
|A'|<5 


lim lim sup lim sup lim sup lim sup (1 — = 0 , 

HI rtO fctO »?t0 ei-O 


and proves 
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SO that (3.30) turns into 

lim sup lim sup lim sup lim sup lim sup Je{tUe^K,r,ri-, w; O) 
r4,o /^4-0 774,0 £4'0 


< lim sup lim sup lim sup lim sup u; 

r4,0 ^4-0 774-0 ^4-0 


r,r}J 


For all i, we have by construction 




e4-0 


SO that, by definition of 


= \Ol, \V{A^ ). 

774,0 £4-0 


Hence, summing over f, 1 < f yields 


^K,r ^K,r 

) = '^\Ol \V{A^ ). 

^ 7?4,o £4,0 ^ 

2=1 2=1 


(3.31) 


On the one hand, VuK,^r = Vu holds on O \ Air- On the other hand, for all i,K,r, S AT := 
conv({A; : 1 < / < k}), which is a compact subset of intdomH. Using in addition the non-negativity 
of the energy density, one may then turn the above equality into 

, r 

^ ^ lim lim J^{u^ ^ J,, ca, ^ A — J{Uf^^ri O) — J (l^i O \ Air) A J{Un^ri -Air) 

^^ 774.0 £4-0 
i=l 

< J{u] O) + lAdj-l sup V. 

K 


Combined with (3.31), this yields 


limsuplimsuplimsuplimsuplimsup j£(tUe^K,r,»7)<^; O) < J{u;0). (3.32) 

tfi T’l-O kIO 774,0 £4,0 

We are now in position to conclude. By coercivity of V, the sequence V(fU£^K,r,77) is bounded in 
L^(0; M™'^'^). Combined with (3.29) (convergence in L°°(0;M™)), this shows that any weakly con¬ 
verging subsequence of {tu£^K,r,ri)e,ri,K,r,t in 1U^’P(0;M™') converges to u. Hence the F-liminf inequality 
of Proposition 3.2 yields 


liminf lim inf lim inf lim inf lim inf K r 775 O) > J(u:0). 

tti 7-4.0 k 4 ,o ?? 4 .o £4.0 ' ’ - -' 

These last two inequalities combine to 

lim sup lim sup lim sup lim sup lim sup (| j£(tn£_K,7-,775 O) — J{u\ 0)| -|- ||tUe^K,r,77 — '^^IIlp(o ;R’")) — 0) 
t'f‘1 7^4,0 7^4,0 774,0 £4-0 

and we conclude as before by the Attouch diagonalization lemma. 


Step 3. Recovery sequence for general functions. 

We claim that, for all ui G 17', all bounded Lipschitz domains O C and all u G kF^’^(0;M™), 
there is a sequence {ue)e C VF^’P(0;M"*) with ^ u in VF^’^(0;M™') and A{ue,uJ]0) J{u]0). 
By the locality of recovery seqnences (cf. Corollary 3.3), we may consider that O is a ball of to 
which we may apply the approximation result of Proposition A.16. By the P-liminf inequality of 
Proposition 3.2, we can further assume that u G VU^’^(0; ffi™') satisfies 

J{u;0)= [ V{Vu{y))dy < 00 , 

Jo 

so that Vtt G domU almost everywhere. Let u be such a function and let cj G 17' be fixed. 

Since O is bounded, Lipschitz and strongly star-shaped, V is convex, and 0 G intdomU, Proposi¬ 
tion A.16(ii) shows that there exists a sequence {un)n of continuous piecewise affine functions with 


STOCHASTIC HOMOGENIZATION OF NONCONVEX UNBOUNDED INTEGRAL FUNCTIONALS 


25 


Vun G intdomy pointwise such that Un ^ u (strongly) in and J{un',0) —> J{u]0) as 

n t oo. By Step 2, for all n, there exists a sequence {u£^n)e C such that Ue^n Un in 

W^'P{U;W^) and Je{ue,ni^]0) —)• J{un',0), as e I 0. In particular, 

lim lim {\Je{Ue,n,^;0) - J{u]0)\ + \\Ue,n “ 'w||lp(0;E-)) 

nfoo etO 

= lim (|J(ti„;0) - J{u;0)\ + \\un - ii||lp(0;K"*)) = 0. 

nfoo 

We then conclude as before by the Attouch diagonalization argument. □ 


3.6. Proof of Corollary 2.4(i): lifting Dirichlet boundary data. We split the proof into two 
steps. We first consider the case when J{au; O) < oo for some a > 1, and then turn to the case when 
in addition M{Wu) < oo or M{aS/u) < oo for some a > 1. 

Step 1. Case when J{au; O) < oo, for some a > 1. 

As u G ri + Vlg^(0;M™) and J{au;0) < oo, Proposition A.16(ii)(a) yields the existence of a 
sequence (ufc)fc <G « + with —)■ u in iy^’P(0;M™') and J{vk]0) —)• J{v,0). For all 

r > 0, set Oj. := {x ^ O : dist(3:,(90) > 2r}, ;= {x G O : dist(x,(90) > r}, and choose smooth 

cut-off functions XriXr with the following properties: the functions take values in [0,1], Xr equals 
1 on 0} and 0 on \ O^, Xr equals 1 on and 0 on \ O, and |Vx);|, \VXr\ ^ C"/r for some 
constant C'. For all a; G Q.' , Proposition 3.10 provides sequences {u^)e and ia 

such that ^ u and Xr^k + (1 “ xl)u in and , 0 J-,O') —)• J{u;0') 

and —>• J{x]-Vk + (1 — Xr)u]0')^ for any subdomain O' C O. We then set : = 

Xr'^erk + (^ “ Xr)'^e ■ Given t G [0,1), Using the decomposition 


t^W^,r,k = tXr'^Ve,r,k + ^(1 “ Xr)V< + (1 - i)“ O 


and convexity, we obtain 

O) < (1 — O) + Je{u^,0J] O \ O^), 

where the error term reads 


(3.33) 






For all y G O \ O^, since x):(y) = 0, we have 

Kr,fe(2/) - <(2/)l < IKr,fc - {xW + (1 - Xr)^^)llL-(O) + IK “ ^^IIl°°(0)- 

By assumption, there is some 5 > 0 with adhi?^ C intdomM. Hence, for all fixed r,k,t, there exists 
£r,k,t > 0 such that for all 0 < e < £r,k,t we have 


Y3^Vx^(<r,fc -<) 


<< 5 , 


L°°( 0 ) 


and therefore 


limsupE'^^^ j < |0| sup M(A') < oo. 
elO ’ ’ ’ |A '|<<5 


Inequality (3.33) then turns into 


lim sup lim sup lim sup lim sup 
ffl fctoo rtO etO 

< lim sup lim sup J{xlvk + (1 ~ xl)^', O) + lim sup J(ri; O \ O^). 

/cfoo r4-0 r4-0 
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The second term in the right-hand side vanishes since J{u;0) < oo, and it only remains to study the 
first term. By definition, for fixed k, we have Vk S so that for all r > 0 small enough, 

xl^k + — xl)u = Vk pointwise on O. This implies 


and thus 


limsuplimsup + (1 — Xr)^; O) = limsup J(ufc; O) = J(u; O), 

/cfoo r4-0 k^oo 


limsuplimsuplimsup limsup ^w; O) < J{v;0). 

tfl /cfcxD r4,0 e:4'0 


Combined with the T-liminf inequality of Proposition 3.2 and a diagonalization argument, this proves 
the hrst part of the statement. 

Step 2. Cases when M{Vu) < oo or M{aVu) < oo for some a > 1. 

If u is chosen in such a way that Jq M{'Vu{y))dy < oo, then we can repeat the argument of Step 1 
with u‘^ := u, and bound the last term in the right-hand side of (3.33) by 

limsup lim sup j£(tt, cj; O \ O^) < lim sup / M{Vu{y))dy = 0. 
r4,0 eto r4,0 J 0\0^ 

We conclude with the case when u satishes M{a'Vu{y))dy < oo for some a > 1. Let Vk,xljXr be 
chosen as in Step 1, and let uj £ Q' he hxed. For any t G [0,1), Proposition 3.10 shows the existence 
of a sequence {v^rkt)s ™ such that Xr'^k + (1 “ xl-)u/t in 

and —)■ J{xlvk + (1 - xl)u/t]0'), for any subdomain O' C O. Set ■= 

Xr^erk t “ Xr)^/^- before, we obtain by convexity 

Je{tW^,r,k,U^) < Je{v'^,r,kp'^'^0) + Je{u/t,UJ]0\Ol) 

+ ( 1 - 


{1-t) [ M 
Jo 


1 -1 


^Xr{y)K,r,k,t{y) - uiy)/t) dy, 


and the conclusion then follows, using the convexity once more in the following form: for t > l/ce, 

limsuplimsup O \ O^) < limsup / M{Vu{y)/t)dy 

rlO elO rlO J 0\0di 


o\op 

<—limsup f M{aVu{y))dy+ l\msup\O\O'^\M{0) = 0 
rtO io\Op riO 


This completes the proof. 


□ 


Remark 3.11. As can be seen in the proof, the assumption that J{au;0) < oo can be relaxed to 
J{au; O') < oo for some open neighborhood O' C O of dO in O. □ 

3.7. Proof of Corollary 2.4(ii): soft buffer zone for Dirichlet boundary data. We split the 
proof into two steps. For all s > 0 and O C we use the notation Og ■= {x £ O : dist(x, dO) > s}. 

Step 1. F-liminf inequality. 

Let uj £ 0,', let O C be a bounded Lipschitz domain, let u £ LF^’P(0;M”*) with J{u;0) < oo, 
and let {ue)e C LF^’^’(0;M™) be a sequence with ^ u in VF^’^’(0;M™). By the F-liminf inequality 
for Jg in Proposition 3.2, 

liminf J^(u£,a;; O) > liminf j£(u£,a;; O^) > J(rt; O^) = / V{Vu{y))dy, 

£4.0 £4-0 Jo^ 

that is, using that V{'Vu{y))dy —>• 0 as r/1 0, 

liminf liminf O) > J{u\0). 

r;4,0 £4-0 


Step 2. F-limsup inequality. 
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Let oj € Q.’, let O C be a bounded Lipschitz domain, and let u G R™) with J{u] O) < oo. 

By Proposition 3.10, there exists a sequence {we)e C LP^’^(0;R™') such that tCg —^ 0 in VP^’P(0;R”*) 
and Je{u + rCe, w; O) —?■ J(rt; O). Given r] > 0, choose a cut-off function Xrj with values in [0,1], such 
that Xri equals 1 on O,^ and 0 outside O, and that satisfies |Vx^| < C/rj for some constant C' > 0. 
Set Vi;^r] ■= X-qWe G VLq^’^( 0; R"*). For all t G [0,1), we have 

tVu + tVVe,n = + We) + t{l - Xr7)Vtt (1 - 


so that by convexity and the definition of 

J^{tu + tVe,n,w-,0) < {1 - t)Ee^r,,t + Je{u + We,u}-,0) + [ {I - Xri{y))V^^’'^{y,Vu{y), u;)dy 

Jo 

< {1-t)Ee^r,,t + Je{u + We,uj-,0) + [ \Vu{y)\Pdy, 

Jo\On 

where the error is defined by 

Es,r,,t ■= 


(3.34) 


< |O^|M(0) + 


IO\Or, 


l-t 


We{v)yXr,{y) 


dy. 


(3.35) 


By the Rellich-Kondrachov theorem, ^ 0 (strongly) in L^(0), so that limsup^ < |O|M(0) 
for all t, y. Passing to the limit in inequality (3.34) thus yields 


lim sup lim sup lim sup J^(ttt-|-tUe, w; O, r/) < lim sup J£(ri-|-uig, w; O) = J{u;0). 
tfl r;4.0 £4.0 £4.0 


We then conclude by the same diagonalization argument as before and Step 1. This proves the first 
part of the statement. 

Now consider the case when u satisfies J{au;0) < oo for some a > 1. Then, for all t G [0,1), 
Proposition 3.10 provides a sequence {we,t)e E VF^’^(0;R"*) such that Wey —^ 0 in W^’^{0]'SE') and 
Je{u/t + We,t-,w-, O) — >• J{u/t] O) as 6 0. Define Ve^t,r] ■= XvWe,t, where Xri is the same cut-off function 

as above. We then have now 

Vu + tVVe,t,ri = tXr^'S/{u/t + We,t) + t{l - Xr))^u/t "h (1 - i) 


SO that by convexity and definition of , 

J^{u + tVe,r„w,0) < (1 - t)E'^ + Je{u/t + We,t,w;0) + t^~P [ \Vu{y)\Pdy, 

JO\Or, 

where the error is defined by 


El„, := 


-'e,ri,t 


:=|O,|M(0)+ [ 

JO\Or, 


t 

1 - t 


We,t{y)'^Xr,{y) 


p 

dy. 


(3.36) 


By the Rellich-Kondrachov theorem, We^t 0 (strongly) in L^(0) for all t, so that limsup^ ^ = 

|O^|M(0) for all t,rj. Passing to the limit in inequality (3.36) then yields 


limsup limsuplimsup (u -|- tVs,t,ri,w‘, O) < limsuplimsup Je{u/t + w^^t, w]0) = limsup J(M/t; O). 
tfl »?4.o £4.0 ’ ’ tti £4.0 ’ tti 


Since u satisfies J{au; O) < oo, we deduce by convexity that the map t ^ J{u/t\0) is continuous on 
(1/a, 1]. This implies limsup(.|-]^ J{u/t]0) = J{u]0), and the conclusion follows. □ 


Remark 3.12. As can be seen in the proof, the assumption that J(au]0) < oo can be relaxed to 
J{au; O') < oo for some open neighborhood O' C O of dO in O. □ 
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4. Proof of the results for nonconvex integrands 

In this section, we study the case when W is nonconvex but admits a two-sided estimate by 
a convex function (which may depend on the space variable) and prove Theorem 2.6. Let W be 
a (nonconvex) r-stationary normal random integrand, which is further assumed to be ru-usc (in 
the sense of Definition 2.5, with respect to some r-stationary integrable random field a). Up to 
a translation, for simplicity of notation, we can restrict to the following stronger version of (2.10) 
and (2.4): for almost all cj, y, we have, for all A, 

< U(y, A,cn) < W(y,A,u;) < C(V(y,A,u;) + 1), (4.1) 

for some C > 0 and d < p < oo, and for some convex r-stationary normal random integrand V. Also 
assume that 0 belongs to the interior of the domain of the convex function M ;= supess^^^ U(?/, ■,u}). 
We can then apply Theorem 2.2 to V, yielding an homogenized energy density V with the following 
property: defining 

Je{u,uj-,0)= [ V{y/e,Vu{y),u;)dy, J{u;0) = [ V{Vu{y))dy, 

Jo Jo 

for almost all w, the integral functionals J^{-,co;0) T-converge to J{-;0) on IU^’P(0;M™'), for any 
bounded Lipschitz domain O C Let Dq C D be a subset of maximal probability on which all 
these assumptions and properties (of V, W) are simultaneously pointwise satisfied. 


4.1. Definition of the homogenized energy density. We need to define in this section a candidate 
for the homogenized energy density W . As before, the standard homogenization formula with Dirichlet 
boundary conditions does not hold because of the generality of the growth conditions considered here. 
Instead, we use the corrector for the convex problem as a boundary condition for the nonconvex 
problem, which is indeed admissible because of the two-sided growth condition (4.1). 

More precisely, for all A G Lemma 3.4 yields a function G Mes(D; (M*^; M”^)) such 

that Vv;a( 0,-) G and 

F(A) =]E[U(0, A+ VyPA(0, •),•)]■ 

Now, for any t G [0,1), consider the function defined by 

p\{0,u]):= inf / W{y,tA + tVipK{y,u:)+ A/v{y),u:)dy. 

(0-,w^) Jo 

As this quantity is stationary and subadditive, the Ackoglu-Krengel ergodic theorem implies: 


Lemma 4.1 (Definition of the homogenized energy density). Let t G [0,1) be fixed. Then, there 
exists a function Wt '■ domU —)• [0, oo) such that, for all A G domU, the following holds for almost all 
w G Dq." for all bounded Lipschitz domains O C 


Wt(A) = lim 


p\{0/e,uj) 

\0/e\ 


(4.2) 


where convergence also holds for expectations. Now define W{A) := liminfi|i liminfA'-j-A 
any A G domU, and further set W{A) = oo for all A ^ domU. Then, W satisfies U < IT < C{l + V) 
on the whole of and for all A G the following holds for almost all uj and all bounded 

Lipschitz domain O C M'^.' 

W (A) = lim inf lim inf lim ^, (4-3) 

tti A'^A £4-0 \0/e\ 

where the lim inf as t f 1 can further be restricted to t ^ Q. Finally, in the particular case when W 
is convex, then W coincides with the various definitions for the homogenized integrand as given by 
Theorem 2.2. □ 
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Proof. We split the proof into three steps. 

Step 1. Definition of Wt{A.) and proof of (4.2). 

First consider the case when A E domi/. Let t E [0,1) be fixed. The upper bound in (4.1) 
then implies E[/i)^(0,-)] < CjOKl + tV{A) + (1 — t)M(O)) < oo. As the function is obviously 
stationary and subadditive, and as •) is measurable by Hypothesis 2.1, the Ackoglu-Krengel 

subadditive ergodic theorem (see e.g. [33, Section 6.2]) can be applied and asserts the existence of 
some WtiA) E [0,oo) such that, for almost all ui, we have 

W,(A) = lim 

I In I 

for any regular sequence {In)n C X := {[a, 6) : a,b E such that lim„.[-oo/„ = (in the usual 
sense of [33, Section 6.2]), and moreover this convergence also holds for expectations. In particular, 
we easily see that the same result must hold for the choice In = nQ^, where Qq is any cube aligned 
with the axes. Further note that, for all bounded Lipschitz subsets O' <Z O C. we can estimate, 
as /e]W^) C Wg’^(O/e;M"^), 

Lo) < e'^fAiO'fe, u) + f W{y, tA + tVipA{y, u}),uj)dy 

J(o\o>)/e 

< e'^y\{0'/e,u]) + C\0 \ O'j ( 1 + / V{y,A + VipA{y,u}),uj)dy] , 

V J{0\0')/e j 

where the last expression in brackets converges to 1 + H(A) < oo as e 0. Now based on this 
estimate, an easy approximation argument (see e.g. [30, Step 4 of the proof of Theorem 3.1]) allows 
us to conclude as follows: for almost all uo (for all w E Da, for some subset Da C D of maximal 
probability, say), we have for all bounded Lipschitz domains O C and all t E (0,1) 

(«) 

Step 2. Definition of W and proof of the bounds H(A) < W{A) < C(1 + 14(A)) for A E domH. 

Let A E domlL be fixed, and let O C be some bounded Lipschitz domain. We dehne W{A) = 
liminfthminfA'->.A 114j(A'). The bounds 14(A) < W{A) < C(1 + 14(A)) directly follow from the 
two-sided estimate (4.1) together with the following equality, for almost all cj, 

F(A)=Fo(A,a;), (4.5) 

where we have defined 

14o(A, w) := liminf liminf lim inf -r 14(y/e, tA' -|-tV(/?A'(y/e)W)-|-Vu(y),w)dy. 
ttl elO ^gvi/o’^(0;lR^) Jo 

Let us give the argument for (4.5). On the one hand, we can estimate 

14o(A, w) > lim inf lim inf lim inf -f V{y/e,tA' + tA' {oj)+ Vv{y),oj)dy, 
ttl A'^A eiO vew^’P(o-,R"^) Jq 
fo Vd =0 

where we have set A(.(a;) := f^VipA'i'/SjUj). For almost all lo, since A(,(a;) —>■ 0, we can write, for 
any k > 0, 

14o(A,a;) > inf liminfliminf inf -r V{y/e,tA' + Vv{y),uj)dy, 

A':|A'-A|<k tfl efO i,evri>P(0;K’") Jq 

fo Vd =0 

so that formula (2.8) yields 14g(A,a;) > infA':|A'-A|<K 14(A'). Passing to the limit k 0, the lower 
semicontinuity of 14 directly gives 14o(^)W) ^ 14(A). On the other hand, the convexity of 14, the 
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Birkhoff-Khinchin ergodic theorem and the definition of give for all t G [0,1] and all A.' G 

lim inf -r V{y/s,tA' + t'VipA'iy/e,u}) +'Vv{y),uj)dy 
etO t;eVl/o-P(0;K"‘) Jo 

< lim '/ V{y/e,A' + V(p\r{y/e,u!),ui)dy + (1 — t)M(O) = V{A') + (1 — f)M(O). 
etO Jo 

Passing to the limit A' —)• A and t J 1, and using the lower semicontinuity of V in the form of 
V{A) = liminfA'-j-A ^(^0) gives ^ ^(^)- The desired identity (4.5) is proven. 

Step 3. Case when A ^ domP. 

For A ^ domP, arguing as in Step 2 above, we can estimate, using the pointwise bound P < W, 

lim inf lim inf lim inf inf -r lP(?//e, fA + fV(/9A(y/e, w) + Vu(?/), > P(A) = oo, 

ttl A'^.A eto ,)eWp’’’(0;R"‘) Jo 

so that (4.3) trivially holds with fP(A) := oo. Moreover, the bounds P < W < C{1 + V) holds as 
well. □ 


Although the definition of the homogenized energy density IP (A) may a priori depend on the choice 
of a corrector ip a, it would follow a posteriori from the F-convergence result that the value of TP (A) is 
independent of that choice. As this independence will actually be useful in the proof of the F-limsup 
inequality (see the proof of Lemma 4.4(c)), we display a direct proof. 


Lemma 4.2 (Independence upon the choice of a corrector). Assume p > d, and let t G [0,1) and 
A G domP be fixed. For almost all lo, given a bounded domain O C if {u^)^ C 1P^’^(0;M™') 
satisfies ||rte||L°°(o) 0 limsup^, P(y/e, A + Vu£(y), t<;)dy < Ca\D\ for all subdomains D cO 

and some constant Ca > 0, then, for all Lipschitz subdomains D <Z O, 

IPt(A)=lim inf -f W{y/e,tA + tVus{y)+Vv{y),u})dy, (4.6) 

efO JD 

where the limit is well-defined. □ 


Proof. Let t G [0,1) and A G domP be fixed. Let a; G be fixed such that (4.2) holds on all bounded 
Lipschitz domains and such that moreover, for all bounded domains D C M'^, 

||ev?A(-/e,^)||L°°(D) ^ 0, / V{y/e,A + VipA{y/e,uj),uj)dy ^V{A), 

Jd 

which follows almost surely from Lemma 3.4, the Sobolev embedding and the Birkhoff-Khinchin 
ergodic theorem. Let {us)s be as in the statement of the lemma. Also denote vf := e<pA{'/£,^) G 
IPjq’^(]R'^;M™'). By the choice of w, the sequence {vf)e satisfies the same properties as on any 
bounded domain, with Ca replaced by = P(A), and moreover, for all bounded Lipschitz domains 
D C 

IPt(A) = lim inf -j' W{y/£,tA-\-tVvf {y)-\-Vv{y),u})dy. (4.7) 

Let D C O he some fixed Lipschitz subdomain. On the one hand, define 

IPf(A,a;; Z?) = limsup inf -r W{y/£,tA-\-tVus{y) -\-Vv{y),oj)dy. (4.8) 

eto vewfiP{D-,R^) Jd 


Given t? > 0, set := {x G D : dist(x, dD) > 77 } and consider the difference 

inf [ W{y/£,tA + tVueiy)+ Vv{y),uj)dy 
veWfiP{D-,K.^)JD 



W{y/e,tA-htVvf{y) + Vw{y),uj)dy. 


(4.9) 


inf 
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Choose a smooth cut-off function x-q such that Xq equals 1 on and vanishes outside D, with 
|Vx? 7 l < C'/r] for some constant C' > 0, and define := Xq^e + (1 ~ Xq)ue- Restricting the first 
infimum in (4.9) to those v’s that are equal to t{vc — u^) on dDrj, we obtain 


< 


inf 


W{vie, tA -h tVw‘^^^{y) + Vv{y),uj)dy. 


JD\Dr, 

Hence, choosing u = 0, using the upper bound W < C{1 + V) and decomposing 

= ^Xq'^V^ + t{l - Xq)'^Ue + (1 “ t) y^Vx, 7« “ ^e), 

we obtain by convexity 




, A -h Vueiy),uj)dy 


where the error reads 


R 


\ JD\Dr, 

+ / y{y/^,-^ + 'yv‘^(.y),<^)dy + ], 

JD\Dr, J 

\t,q ^ ^ ^y- 


Since and go to 0 in L°°(Z); MR), we can prove that, for any t £ (0,1), 

limsup limsup Ag j < limsupC'|L) \ A>r)|(l + C'a + C'a) = 0- 
r; 4,0 etO f)t 0 

In view of equalities (4.7) and (4.10), this implies W[{A,u}; D) < Wt{A). 

On the other hand, define 


W^ (A, w; D) = lim inf inf 

AO ,;eVEo’^(D;R"*) Jd 


W{yje, tA + tVueiy) + 'yv{y),ui)dy, 

and repeat the same argument as above with := {x : dist(x, D) < t/} and 

f W{y/e,tA + tVv‘^{y) + Vv{y),uj)dy 
I J D‘n 

inf f W{y/e,tA + tVueiy) +Vw{y),uj)dy, 
Vq'^{D;RR jD 


(4.10) 


A"^, := 

e,t,q 


inf 


«;€Wn 

which then yields W'l{A,uj;D) > Wt{A). This shows that W'l{A,uj',D) = w\{A,uj-,D) = Wt{A), and 


the result is proven. 


□ 


Let rii C ffo be a subset of maximal probability such that (4.2) holds for all a; G fli, t G Q O [0,1) 
and A G n domR, such that (4.3) holds for all cu G ffi and A G and such that we further 

have, for all cj G ffi, A G and all bounded domains O C M'^, 

V{A) = lim'f V{y,A + V(pA{y,u}),ui)dy. 

AO Jo/6 

4.2. T-liminf inequality by blow-up. In this section, we prove the T-liminf inequality for Theo¬ 
rem 2.6 by adapting the blow-up method introduced by [25] (see also [7, Section 4.1] and [14]). In 
the present context, a subtle use of the corrector for the convex problem is needed. 

Proposition 4.3 (T-liminf inequality). For any u G ffi, any bounded Lipschitz domain O C 
and any sequence (rie)e C fT^’^(0;M"*) with ^ u in fT^’^’(0; M™'), we have 


liminf/e(rie, w; O) > l{u]0). 
eJO 


□ 
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Proof. For all r > 0 and x € define Qr{x) = x + rQ and Sr,K.{x) = Qr{x) \ QrK{x) for all k > 0. 
For all e > 0 and all A, ui, define Xe A — ^Pa {-For all cj G fii and A G the sequence 

(Xg a)^ satisfies x^ a —^ 0 in W^'P{0]W^) and J^{Xe A + A • x, tu; O') —)• J(A • x; O') = |0'|y(A) as 
e 4, 0, for any subdomain O' C O. 

From now on, let a; G be fixed, let O C be some bounded Lipschitz subset, and let {u^)^ C 
VF^’^(0; M”*) be some fixed sequence with ^ rt in bF^’P(0;M™). We need to prove 

liminfw; O) >/(tt; O). (4-11) 

etO 

It does not restrict generality to assume liminfe/^(ue, w; O) = lim^a;; O) < oo and also 
supgw; O) < oo. Hence, Vu^{x) G domVF(x/e, •, w) = domH(x/e, •, w) for almost all x. Fur¬ 
thermore, the F-convergence result for V yields J{u',0) < liminfe j£(u£,a;; O) < I^{u^, 0 J',O) < 

oo, so that \/u{x) G domH for almost all x. 

Step 1. Localization by blow up: we prove that it suffices to show for almost all xq that 

liminf liminf lim-r W{y/e,tVue{y),io)dy >W{Vu{xo)). (4-12) 

itl etO JQ^(^xo) 

For all e > 0, consider the positive Radon measure on O defined by dp^(x) = W{xje, \/ufix),oj)dx. 
As sup£P£(adhO) < oo by hypothesis, the Prokhorov theorem asserts the convergence ^ p np 
to extraction of a subsequence, for some positive Radon measure p on adhO. (The extraction will 
remain implicit in our notation in the sequel.) By Lebesgue’s decomposition theorem, we can consider 
the absolutely continuous part of the positive measure p, and the Radon-Nikodym theorem allows to 
define the density / G L^(C/) of the latter. As O is open, we then have by the portmanteau theorem 
(see e.g. [10, Theorem 2.1]) 


liminf/^(ue;, w; O) = liminf/ 0 £(O) > p{0) > / f{x)dx. 


'O 

Hence, in order to prove (4.11), it suffices to show that f{x) > Wifi7u{x)) for almost all x. Since 
p(adhO) < oo, we have p{dQr{x)) = 0 for all r G (0,1) \ Dx, where Dx is at most countable, so that, 
for almost all x, Lebesgue’s differentiation theorem and the portmanteau theorem successively give 

/(^)= ii.„ 

riO riO elO 

r^Dx r^Dx 

Hence, it suffices to show that, for almost all xq, 

liminflirn'T W{y/e,Vue{y),oj)dy > LF(Vu(xo)). 

HO elO Jq^(^xo) 

Using the ru-usc assumption on W, we easily deduce the following inequality: 

limsupliminf lim-r VF(y/e,tVue(y),a;)(iy < liminf lim-r W{y/e,'S/u^{y),(jj)dy. (4.13) 

t^l r-to eto Jq^(xo) HO eiO JQ^(xo) 

Indeed, as Vufiy) G domVF(y/e, -, 0 ;) for almost all y, we can write 

/ W{y/£,tVue{y),uj)dy < {1 +A^r{t)) J- W{y/e,Vue{y),uj)dy 

J (5r(^o) Qr(^o) 

+ ^wii)r a{y/£,u})dy, 

d Qri^o) 

and thus, by r-stationarity of a, the Birkhoff-Khinchin ergodic theorem yields 

liminflirn'T W{y/£,t'S/Ue{y),oj)dy 
HO el-O Jq^{xo) 


< (1 + A^(t)) liminf lim 


/ 


HO elO Jq^{xo) 


W{y/£,Vus{y),u})dy + A^(f)IE[o(0, •)], 
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SO that inequality (4.13) now directly follows from the ru-usc assumption with respect to a (meaning 
indeed that limsup^^^ A^(t) < 0). Using (4.13), we finally conclude that it is sufficient to show (4.12) 
for almost all xq. 

Step 2. Proof of (4.12) by truncation. 

The idea is to truncate at the boundary, in order to make appear in the left-hand side of (4.12) 
precisely W{tVu{xQ)), which will then allow us to conclude. 

Let t,K ^ (0)1) bs fixed. Since p > d, the Sobolev embedding yields tie —)• tt in L°°(0;M™'). 
Moreover, combining the Lebesgue differentiation theorem for Vu and the Sobolev embedding for 
p > d, we deduce that, for all xq ^ M (for some null set M C |AA| = 0), 

lim -||u - u(xo) - Vu(xo) • {x - Xo)||L°°(QUa;o)) = 0- (4.14) 

r4,o r ^ ^ " 


Enlarging the null set AA, we can also assume that \i'u{xo) G domU for any xq ^ Af. From now 
on, let xo G O \ AA be fixed and write for simplicity A := Vu(xo)- Since V is convex and lower 
semicontinuous, we have U(A) = liminfy\^/_j,^ U(A'), and a diagonalization argument then 

allows us to choose a sequence (Aj.)^ C such that A^ —)• A and U(Ar) —)• U(A) as r 0, and 

simultaneously 


l™-|k - u(xo) -Ar-(x- xo)||l°°(q,(xo)) = 0 - 


(4.15) 


Let c/>r^K be a smooth cut-off function with values in [0,1], such that (^r,K equals 1 on QrK(xo), vanishes 
outside Qr(xo), and satisfies ||V())r,K||L°° < r(i-K) ' then set 

^e,r,K •— 4 ^ r , K'^s T (1 4 ^ r , K ){ d ^( xo ) T Aj. • (x Xq) T X£:,Ar.(^))' 

Since v^^r,K, coincides with on QrK{xo) and 0 < VU <C{1 + V), we have 


i 


Qr{xo) 


lU (y/e, tVve,r,K{y)^^)dy 


< 


Qr(xo) 


W{y/e,tVu^{y),u})dy + ^ [ V{y/e,tVve,r,Kiy)^^)dy + C{1 - (4.16) 

^ JSr,K{xo) 


Defining '■= '^(pr,K{x) <S> (u^ix) — u(xo) — Ar ■ (x — xq) — a^) decomposing 

t 






t'^Ve,r,K = -V t{l - 4>r,K){Ar + Vx^^aJ + “ ^)y_ ^ 

we obtain by convexity of V 

V{y/e,tVve,r,K{y)^^) ^ t(j)r,KV{y/e,Vus{y),uj) +t{l - (j)r,^)V{y/e, Ar + Vxe,A^{y),uj) 

+ (1 - tW (^y/^: 

< W{y/e,Vueiy),uj) + V{y/£,Ar + Vxe,Ariy)^^) 

+ (1 - tW (y/^, 

Combined with 


£,?-,K||L°°(5r,«Uo)) < (^Ike - '“IIl°°(0) + Ik - '“(^^o) - Ar - {x - Xq) ||l°°(Q,.(xo)) 


\ u;, — k) \ 

+ IIXe,AAlL°°( 0 ) ^) 

the convergences ^ u and Xe a^- 0 ia L°°(0;M™) as e Y 0, and (4.15) yield 

limsuplimsup „(xo)) = 0- 

r4,0 £4-0 







34 


M. DUERINCKX AND A. GLORIA 


By assumption, we can find <5 > 0 with adhi?^ C intdomM. Hence, for all t,S (0,1), there exists 
> 0 such that, for all 0 < r < r^,*, there exists some > 0 such that for all 0 < e < £r,K,t 




< 6 . 


L°°(S,,„(xo)) 


This implies 


Sr,n (^0 ) 


V vie., 


1 -1 


^£,r,K(l/),w ) dy < |5',.,K(a:o)| sup M(A') = r'^(l - k'^)| sup M(A'), 




|A'|<5 


|A'|<5 


where the supremum is finite, by virtue of the convexity of M and our choice of 5 > 0. Combined 
with inequality (4.17) and the definition of the correctors Ar lim^ H(Aj.) = H(A) < oo), this 

yields 


< 


lim inf lim inf liminf lim inf / V{y/e,t\7vs r K{y),^)dy 

Ktl ifl r-4.0 eio J Sr,^{xo) 

lim inf lim inf lim inf / W{y/e,\7us{y),uj)dy. 

Ktl rto eto JSrA^o) 


This turns inequality (4.16) into 


lim inf lim inf lim inf lim inf -r 
Kfl HI rtO etO Jq^{xo) 


W{y/e, tVve,r,K{y)^^)dy 


< lim inf lim inf lim inf 


W{y/£,tVue{y),u})dy 


HI HO HO Jq^(xo) 

C f 

+ lim inf lim inf lim inf—7 / W{y/e,'Vui;(y),uj)dy. 
HI HO HO r'=‘ J Sr,^{xo) 


(4.18) 


Since we have chosen cj S Hi, t € Q Cl (0,1) and A,. S (4.2) holds and reads 

ITt(Ar) = liminf inf -r lT(y/e,fAr + fV(/?A,.(-/e,w) + Vu(y),a;)(iy. 

HO JQrixo) 

Hence, since Ve,r,K. - u{xq) - K ■ {x - xq) G xt,Ar + WQ’^{Qr{xo);W^) with Xe,Ar = ^¥^A,(-/e,w), 
(4.18) yields 


lT(Vu(xo)) = VH(A) < liminf liminf VH 4 (Ar) 

HI f'lo 


< lim inf lim inf lim inf -r 

HI y’lO HO j Qr^xo) 


W{y/e, tVve,r,K{y),u})dy 


< lim inf lim inf lim inf -r W{y/£,t'S/u£{y),uj)dy 
HI y’lO etO J Qv^xq) 


(4.19) 


+ lim sup lim sup lim sup 


C 


L 


W{y/e,Vus{y),uj)dy. 


HI HO HO ^ JSr,Kixo) 

It remains to get rid of the second term of the right-hand side of (4.19). By the portmanteau theorem. 


lim sup ■ 


1 


eiO X JSr,^(xo) 


L 


W{y/£,Vue{y),io)dy = lim sup 

ei-O 


Pe{Sr,K.{xo)') 


P£(adhS'r,K(xo)) ^ /9(adhS'r,K(xo)) 
< hmsup-j-. 

HO ^ 


Since the singular part of p must be supported in a closed subset of adh O of measure 0, we deduce, 
for almost all xq G O/Af, the existence of some tq > 0 sufficiently small such that adhQr(xo) has no 
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intersection with that support for all 0 < r < tq. Hence, for all 0 < r < ro, 
lim sup I 

£- 1-0 ^ J Sr,K.{xo) 

p{a.dh.Sr,K{xQ)) 


< 


W{y/e,Vue{y),u})dy 
\ f f{y)dy = i fiy)dy - i 

^ «/ Sr. K.( XC)^ J OrixC)^ J (Z 


Qrn (^0 ) 


fiy)dy, 


' Sr^nixo) 

where, for almost all xq, the right-hand side converges to (1 — k'^)/(xo) as r -1, 0 by Lebesgue’s 
differentiation theorem. Hence, for almost all xq, (4.19) turns into 


iy(Vu(xo)) < lim sup lim inf lim inf-r W{y/e,tVue{y),u:)dy, 
t-fl r-l-O £-1-0 Jq^{xq) 

and the desired result (4.12) is proven. 


□ 


4.3. F-limsup inequality with Neumann boundary data. In this section, we prove the F-limsup 
inequality, first considering the affine case, and then deducing the general case by approximation. For 
this approximation argument to hold, we would however need to know a priori that the homogenized 
energy W satisfies good regularity properties (i.e. lower semicontinuity on and continuity on 

intdomH). Since this is not clear at all a priori, our strategy (inspired by [7]) consists in introducing 
some relaxations of W that enjoy the required properties, and then in deducing a posteriori from 
F-convergence (or a weaker form of it) the equality of W with its relaxations (so that W itself has 
all the desired properties). Motivated by the work of Fonseca [24] (see also [7]), we thus consider the 
following relaxation of W : 

ZW{K) := inf ^-j' VF(A -|- \/(j){y))dy : cj) continuous piecewise affine on O and 4>\do = 0 

where the definition does clearly not depend on the chosen underlying (nonempty) bounded Lipschitz 
domain O C Also write ZW for the lower semicontinuous envelope of ZW (defined by ZW{K) : = 
liminfA'^-A-2^kF(A^) for all A). Now define the integral functionals corresponding to all these relaxed 
integrands: for any bounded domain O C and u € VF^’^(0; M™), 

ZI{u;0)-.= [ ZW{Vu{y))dy, ZI{u-,0) := [ ZW{Vu{y))dy. 

Jo Jo 

The following result gives some properties of these relaxations, which will be crucial in the sequel: 



Lemma 4.4 (Properties of relaxations). Assume p > d. Then the following holds: 

(a) ZW (and thus also ZW) is continuous on intdom^^VF. 

(b) V <ZW <W <C{l + V). 

(c) W and ZW are ru-usc. 

(d) For any t £ (0,1), we have fadh dom^^VF C intdom^^VF, and the following representation result 
holds: 

ZW{A), 
limj-j-i ZW (tA), 
oo, 

where, in particular, the limit in the second line does exist. 

(e) Let A £ O C be a bounded Lipsehitz domain. Then, there exists a sequence 

{4>k)k C VFq ^(0;M”^) of piecewise affine functions such that —)• 0 in L°°(0;M™') and 

lim / TF(A + V(fkiy))dy = ZW{A). 
ktooJo 


ZW(A) = liminf ZlF(tA) = < 


if A G int domZW; 
if A G d domZVF; 
otherwise; 


□ 
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Proof. Continuity of ZW on mtdom ZW is a result due to [24], which yields part (a) (even without 
any ru-usc assumption on VC). The inequalities stated in part (b) directly follow from the definitions 
of V, W and ZW. Part (e) is standard (see [7, Proposition 3.17] for details). It remains to prove 
properties (c) and (d). 

Step 1. Proof of (c). 

We first prove that Wt is ru-usc, for any fixed t £ [0,1). Let s > 0, A G domVC = doml4, and 
let O C be some bounded Lipschitz domain. For almost all cu, note that by convexity, for all 
subdomains D <Z O, 


lim sup 
etO 


D 


V{yfe, sk + sVipAiy/e, uj),uj)dy 


< lim 
ef-O J £) 


V{y/e,A + V(/?A(y/e,a;),a;)dj/ + (1 - s)M(O) = 17(A) -|- (1 - s)M(O) < oo. 

I 

Hence for almost all uj we can apply equality (4.6) at sA with u,. = se(/?A(-/e, w), which yields 

f W (y/e, stA -V stVipA{y/e, u) + Vv{y),uj)dy. 

I Jo 


VFi(sA) = lim inf 

elO t,eW'o’P(0;I 

Given cu G 11 such that this convergence and (4.2) both hold, and choosing a sequence (h 
Wq^’^( 0; such that 


£ ye 


C 


VFt(A) = lim 

efO Jq 


/ W{y/e, 


tA + tVipA{y/£,uj) +Vvf{y),uj)dy, 


we deduce 


VFt(sA) - VFt(A) < lim 

£t0 Jq 




, s{tA + tVifAiy/e, io) + Vvf (y)), w) 

- W{y/£,tA + tVipA{y/£,u}) + Vv‘f (y), uj)^ dy 


< A“ 


{a{y/£,uj) + W{y/e,tA + tV(pAiy/£,uj) +Vvf{y),ij))dy 

eio Jo 


= A^(s)(E[a(0,-)] + VFt(A)), (4.20) 

using the Birkhoff-Khinchin ergodic theorem for the stationary field a. As Wt and the field a are 
nonnegative, as a := E[a(0, •)] is finite and as limsupj,.|.;^ A^(s) < 0 by assumption, we deduce that 
Wt is also ru-usc. Now rewriting inequality (4.20) in the form 

VCi(sA) < «A^(s) + (1 + (- 1 ) V A^(s))Wt(A), 

and taking the suitable liminf, we directly deduce VF(sA) — VF(A) < (—1) V A^(s)(a -|- VF(A)) for 
all A G doml7 and s G [0,1), proving that W is itself ru-usc with A^ = (—1) V A^(s). 

We now show that ZW is also ru-usc. Take s > 0 and A G doml7. By definition, there exists a 
sequence of piecewise affine functions {(j)k)k C VFq’^(O) such that 

f 


ZW(A) = lim 

fctoo Jo 


w{A + VMy))dy. 


As A G doml7, the left-hand side is finite, and we can thus assume A J- 'V(pk G domVF almost 
everywhere. Hence the ru-usc property satisfied by VF gives 

ZW{sA) - ZW{A) < lim / (VF(s(A + Vcfkiy))) - W{A + V4>k{y)))dy 
k^ooJo 


< AUs) hm / (a + W{A + V<fk{y))dy = AUs){a + ZW{A)) 


k^OO J Q 


Step 2. Proof of (d). 
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Since domZW = domi/ is a convex set containing 0, it is clear that, for all t S [0,1), ta.<dh.dovD.ZW 
is contained in int domZW. We first show that the limit limt-i-i ZW(tA) exists for all A £ adhdoml/. 
Given some fixed A £ adhdoml/, choose two sequences Sn t 1 tn t 1 with tnjsn t 1 such that 

lim ZW (snA) = lim inf ZW (tA) and lim ZW(tnA) = lim sup ZW (tA). 

n^oo nfcxD 

As SnA,tnA £ doml/ for all n, and as ZW is ru-usc, we have 

lim ZW{tnA) < limsup(a + ZW(s,iA))A^(t„/s„) + lim ZW{snA) 

ntoo n\oo >^100 

< limZW(snA) < YiuiZW{tnA), 

nfoo nfco 

which thus proves the existence of the limit limj-i-i ZW{tA) for all A £ adhdomy. 

We now prove the claimed representation result. First, if A £ int doml^, then liminft_>.i ZW(tA) = 
ZW{A) = ZW{A) follows from part (a). Second, if A ^ adhdoml/, then ZW{tA) = oo for any t 
sufficiently close to 1, and thus liminf^^i ZW{tA) = oo = ZW{A). Now it only remains to consider 
A £ ddova.ZW. Then ZW{tA) = oo whenever t > 1, so that we simply have ZW{tA) = 

liminft'i'i ZW{tA) = limi'i^i ZW{tA), since we have already proven the existence of this limit. Hence, 
it suffices to prove that ZW{A) = liminf(.|-i 2^W(tA). By definition of the lower semicontinuous 
envelope ZW of ZW, this equality would follow if we could show that, for any sequence A„ —)• A, we 
have 

liminfZVF(A„) > liminf ZW(fA). (4.21) 

nfoo ffl 

ft is of course sufficient to assume liminf^ 2^VF(A„) = lim„2^W(A„) < oo and sup„2^W(A„) < oo. 
Hence, A„ £ domH for all n, and thus, for all t £ [0,1), tA £ intdomH, so that, using part (a) as 
well as the ru-usc property satisfied by ZW, we have 

ZW{tA) = lim ZW{tAn) < lim ZW{An) + A^(f) lim (a + ZW(A„)). 

n'foo n'^oo nfoo 

This yields 

lim sup 2TT(f A) < liminf ZVF(A„), 
ttl 

and proves (4.21). □ 

Combining the T-liminf inequality for towards I with a T-limsup argument, we prove the 
following a priori surprising equality of W with its relaxations. 

Lemma 4.5 (Regularity of the homogenized energy density). Assume p > d. Then VF(A) = 
ZW{A) = ZW{A) for all A £ In particular, W is lower semicontinuous on and is 

continuous on mt domV. □ 

Proof. We split the proof into four steps. 

Step 1. Recovery sequence for /(A • x; O). 

Let A £ int domR and let t £ [0,1). In this step, for almost all ui, for all bounded Lipschitz domain 
O C M'^, we prove the existence of sequences f 1, A^ —>• A and {w^)^ C Hq’^( 0;R”^) such that 
0, rce 0 in fT^’P(0;M”") and Ie{teAe ■ x + £te<PA^{-/e,uj) +We,u]\0) |0|W(A) = 

I{A-x;0). _ 

By definition of W and a diagonalization argument, for almost all u and all bounded Lipschitz 
domains O C M'^, it suffices to prove the existence of a sequence {vs)^ C Wq’^( 0;M™) such that 
Uj —^ 0 in W^’P{0] W^) and leitA ■ x + et(p\(-/£, u) + v^, uj; O) —)• jOj Wt(A) as e 0. 

Let O be some fixed bounded Lipschitz domain. Given e > 0, consider the cubes of the form 
k{z + Q), z £ that are contained in O/e, and denote by Zj £ Z'^, j = 1,... the centers 

of these cubes (the enumeration of which can be chosen independent of £,k). Since O is Lipschitz, 
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we have ^ \0\ as e 0, for all k. For all j,uj, we can choose a sequence with 

G WQ’^{k{zj + such that 


i. 


k{zj+Q) 


W{y, tA + tVipAiy, uj) + Vvj:‘^{y),uj)dy 


1 

< T + 


inf 


{k{zj+Q)-,'M.'^) Jk{zj+Q) 


/ 

Jk( 


W{y, tA + tVipA{y, uj) + Vv{y),uj)dy. 


For all s,k,LJ, we then consider the function := ''^k‘^^k{zj+Q) S Wq^{0/£;W^), and we 

define := € VFq’^( 0; M”^). Up to a diagonalization argument, it suffices to show that, 

for almost all oj (independent of the choice of O, as it is clear in the proof below). 


A,P( 


limsuplimsup (\l^{tA ■ x + et(/?A(-/e, uj) + w‘^f., uj; O) — |0|lUt(A)| + ||rcj 

t-l-rvA I n V ’ 


k^oo £ 4-0 
First we argue that, for almost all uj, 


limsuplimsup/e(tA • X + O) — | 0 |VFt(A). 

kfoo e4-0 


= 0. (4.22) 

(4.23) 


Indeed, by definition of 

Ie{tA - X + et^pAi■/£,uJ) +w‘^^i,,uj;0) 

Ne,k 


< 


{ek)^N^ k + {ek)^ V] , inf 

< ^ ^TTrl,p 


/ 

)Jk( 




4(0/U\U,'il'' k(zj+Q) 

Since W < C{1 + V), the last term of the right-hand side goes to 0 as e 0 for almost all oj by 
construction of the cubes k{zj -|- Q) and definition of (pA- The Birkhoff-Khinchin ergodic theorem 
(which we apply to a measurable map by Hypothesis 2.1) then gives, for almost all oo, 


i;GVEo’^(fc(2j+Q);R™') Jk{zj+Q) 

W{y, tA + tVpAiy, uj),uj)dy. 


W{y, tA + tVpAiy, oj) + Vv{y),uj)dy 


limsup Ie{tA ■ X + etpAi'/^, oj) -|- w‘^uj; O) 
e4.0 


< ^ + |0|E 


inf W{y,tA + tV(pAiy,-)+ '^viy),-)dy 

v&WQ'^{kQ-,R^) JkQ 


(4.24) 


Lemma 4.1 then yields the desired result (4.23) as A; f oo. On the other hand, by dehnition (4.2) of 
Wt, for all k and almost all oj, we have 

liminf Ie{tA ■ x + etpAi-/£,oj) -|- uj^^k{'iOj),uj; O) 

eiO 

>|0|liminf inf W{y,tA + tVpA{y, 0 j) + '^v{y),uj;O) = \O\Wt{A). (4.25) 

efo Jo/s 

We now show that w‘^f. 0 in L^(0;M”*) as e ), 0, for almost all oj. Combining inequality (4.24) 

with the bound W < C(1 -|- V), the p-th order lower bound for W and the convexity of V, we indeed 
have 

limsup ||tA -I- tVipA{-/£,oj) + '^'>JJs,k\\LP(o) - + ^1^1(1 -I- U(A) -I- (1 - t)M(O)) < oo. 

For almost all oj, the weak 17 convergence of the sequence (V<y9A(-/e, oj))e to 0 implies the boundedness 
of this sequence in L^(0;so that (Vrt;^^)^ is also bounded in L^(0; M™’^'^), for any fixed k. 
By Poincare’s inequality on cubes of side length ke, this implies 

ll^£,fcllLP(0) < Ck{uj)£, 

for some (random) constant Ck{uj). Combined with (4.23) and (4.25), this proves (4.22). 
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Step 2. Recovery sequence for ZI{A ■ x‘,0). 

Let A S int domR and let O C be a bounded Lipschitz domain. In this step, for almost all 
UJ, we prove the existence of a sequence {us)e C such that Ue —^ 0 in and 

/^(A ■ X + Us, to; O) —)• ZI{A ■ x]0). 

Lemma 4.4(e) gives a sequence {4>k)k C VLq’^( 0; M"*) of piecewise affine functions such that (j)k ^ 0 
in LP(0;M™') (and even in L°°(0;M"*)), and 

I{<Pk + A-x-,0)= [ W{A + V(l>k{y))dy ^ \0\ZWiA). 

Jo 

Denote by the partition of O associated with the piecewise afhne function (pk- For all k and 

1 < z < nfc, considering A^ := A + 'S/pklpi on Step 1 above gives, for almost all a;, a sequence 

ts 'I 1, a sequence A^ ^ A^ and a sequence (u* j^)s C LFq^’^(P^; M™') such that ^ (-/e, a;), u* ^ ^ 0 

in W^’P(Pl;R'^) and 4(pA*_^• x + ^(-/e, w)+ u*w; P|) ^ I{Ai-x;Pl) = /(4 + A-x;P^). As 

the u* ^’s satisfy Dirichlet boundary conditions, they can be directly glued together, while for the 
we need to repeat the more complicated gluing argument of Step 2 of the proof of Proposition 3.10, 
with p > d. Although the functional Is is not convex here, as in the proof of Proposition 4.3, the idea 
is to use the bound W < C {1 + V) at all points where the cut-off functions are different from 1 or 
0, then use the convexity of V and estimate the corresponding error terms as before. We leave the 
details to the reader. 

Step 3. Recovery sequence for ZI{A • x; O). 

Let A £ domR and let O C be a bounded Lipschitz domain. In this step, for almost all w, 
we prove the existence of a sequence {Us)£ C VF^’^(0;M”^) such that ^ 0 in VF^’P(0;M"*) and 
Is{A ■ X + Us,uj; O) —>• ZI{A ■ x; O). 

By Lemma 4.4(d), ZW and ZW coincide on intdoniR, and hence the result on intdomR already 
follows from Step 2. Let now A £ d(iom.ZW. Lemma 4.4(d) then asserts ZW{A) = ZW{tA). 
By convexity of domR, for all t £ [0,1), we have tA £ intdomR, and hence, for almost all u, 
Step 2 above gives a sequence {us^t)s C W^’P{0]W^) such that Us^t —^ 0 in VF^’^(0;M™') and Is{tA ■ 
X + Us,t,i^',0) —)• ZI{tA ■ x;0) = \0\ZW{tA). The conclusion then follows from a diagonalization 
argument. 

Step 4. Conclusion. 

Let A £ domR, let O C be a bounded Lipschitz domain, and let {Us)s be the sequence given 
by Step 3 above. As ^ 0 in iy^’^(0;M™), the P-liminf inequality (see Proposition 4.3) gives, for 
almost all cj, 

\0\ZW{A) = hm4(A • X + <, w; O) > /(A • x; O) = |0|Ty(A). 

e4,0 

This being true for any A £ domR, we conclude that ZW = W everywhere. □ 

With Lemma 4.5 at hands, we may prove the P-limsup inequality. 

Proposition 4.6 (P-limsup inequality). Assume p > d. There exists a subset H' C Di of maxi¬ 
mal probability with the following property: for all oj £ , all strongly star-shaped (in the sense of 

Proposition A.16) bounded Lipschitz domains O C and all u £ LF^’P(0; M™'), there exist a se¬ 
quence {us)s C VF^’^(0;M™') and a sequence {vs)s C LIq’^(0;M™') such that Us ^ u and Vs ^ 0 in 
LF^’^(0;M™'), and such that Is{us + Vs,uj;0) ^ l{u]0) and Js{us,u};0) ^ J{u]0) as e | 0. □ 

Recall that the P-liminf inequality implies the locality of recovery sequences (see the proof of 
Corollary 3.3). Hence, the P-convergence result on a Lipschitz domain D for Neumann boundary 
conditions follows from the P-limsup on a ball B D D and the P-liminf inequality on B\D. For the 
adaptation of Corollary 2.4, the approach is similar and we leave the details to the reader. 


Proof. Step 1. Recovery sequence for affine functions. 
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In this step, we consider the case when u = A ■ x. is an affine function. More precisely, we 
prove the existence of a subset Q' C fli of maximal probability with the following property: given 
a bounded Lipschitz domain O C for all to G O' and all A G intdomld^, there exist sequences 
{ue)e C and {ve)e C bFo’^(0;R”") such that -- A • a; in M'") and 0 

in iy^’P(0; R™), and such that + 'Ce,to; O) —)• /(A • x; O) and O') —)■ J(A • x; O') for all 

bounded domains O' C R*^. 

Let A G int domIL. For almost all to G Oi, and all bounded domains O' C R'^, we have by convexity, 
the Birkhoff-Khinchin ergodic theorem, definition of and continuity of V at A: 

limsup lim JeitA' ■ x + stip\i{-ls, to), to; O') < limsup lim Je^A' ■ x + to), to; O') 

ttl,A'—^-A sfO A'—^-A stO 

= \0'\ lim F(A') = |0'|F(A) = J(A • x; O'). 

A'-^A 

Combined with the F-liminf inequality for J^{-,uj‘,0') towards J[-]0') (for to G Oi), this yields 

lim Wm Ji.{tA' ■ X + /£,uj), 0 j;O') = J{A ■ x-.,0'). (4.26) 

tti,A'->-A eto 

By definition of VF, we may choose sequences A„ —)• A and fn t 1 such that Wt^{An) —)• W{A). For 
this choice, (4.26) yields for almost all to and all bounded domain O' C R'^ 

lim limj£(f„A„ • x + etn<y^A„ (-/e!O') = J(A-x;0'). 

nfoo £4.0 

For all n and almost all to, set ^ := t^An ■ x + etn'^Kn ('^) ■ By Step 1 of the proof of Lemma 4.5, 
for any bounded Lipschitz domains O C R'^, there exists a sequence (x^ „)£ C VFq’^(O; R"*) such that 
- 0 in bFi’P(0;R-) and /,«„ + to; O) ^ \0\WtSAn). 

By a diagonalization argument, we then conclude that for almost all to and all bounded Lipschitz 
domains O C R'^ there exist sequences {ui.)e C IFjof (R‘^;R'”) and {ve)e C 1F)’^(0;R'") such that 
Xg —^ A • X and Xg —^ 0 in and such that /g(xg + Xg,to;0) —)• /(A • x;0) and Js{ue,i 0 ]O') —)• 

J(A • x; O') for all bounded domains O' C R'^. 

Now define fl' C fli as a subset of maximal probability such that this result holds for all A G 
Qmxd Pi intdomC and all to G 14'. Arguing as in the end of Step 1 of the proof of Proposition 3.10, 
and using the continuity of both W and V in the interior of the domain (see Lemma 4.5), the 
conclusion follows. 

Step 2. Recovery sequence for continuous piecewise affine functions. 

We now show that, for any tx G 14', any bounded Lipschitz domain O C R'^, and any continuous 
piecewise affine function x on O with Vx G intdomR pointwise, there exist a sequence (xg)g C 
4F^’P(0;R"*) and a sequence (xg)g C IFq^’^( 0; R™) such that Xg —^ x and Xg —^ 0 in iy^’^(0; R™), and 
such that /g(xg + Xg, ex; O) —)• /(x; O) and Je{us, ex; O) —)• J(A • x; O). This follows from an immediate 
adaptation of Step 2 of the proof of Proposition 3.10. Again, the functional Ig is not convex, but we 
may use the bound W < 0(1 + V) at all points where the cut-off functions are different from 1 or 
0, and use the convexity of V to estimate the corresponding error terms. We leave the details to the 
reader. 

Step 3. Recovery sequence for general functions. 

We show that, for all cx G 14', all strongly star-shaped bounded Lipschitz domains O C R*^, and all 
X G VF^’P(0; R"*), there exist a sequence (xg)g C IF^’^(0;R™') and a sequence (xg)g C IFq^’^( 0;R™) 
such that Xg —^ X and Xg —^ 0 in kF^’^(0;R™), and such that /g(xg -|- Xg,fx;0) —)• /(x;0) and 
Jg(xg,tx;0) —)• J{u;0). Let O C R'^ be some fixed strongly star-shaped bounded Lipschitz domain. 
By the P-liminf inequality of Proposition 4.3, we can restrict attention to those x G IF^’^(0;R™') 
that satisfy 

I{u-,0) = [ W(Vu(y))dy < oo, 

Jo 

so that Vx G domR almost everywhere. Let x be such a function and let cx G 14' be fixed. 
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Let t € (0,1). Since O is Lipschitz and strongly star-shaped, and since W is lower semicontinuous 
on continuous on intdomlL, ru-usc, and satisfies V <W < C(1 -|- V) (see indeed Lemmas 4.4 

and 4.5), the nonconvex approximation result of Proposition A.16(ii)(c) yields a sequence {un)n of 
continuous piecewise affine functions such that Un ^ u (strongly) in VL^’^’(0; M™'), I{un', O) —)• I{u; O) 
and J{un',0) —)• J{u]0) as n f oo, and such that G intdomf/ pointwise. Now Step 2 above 
gives, for any n, sequences (rt£,n)e C LP^’^(0; M™) and {ve,n)e C Wq{0] M™) such that Un and 

Ve,n —^ 0 in M™'), and such that Ifr{ue,n + Ve,n-, w; O) —)• I{un', O) and J^{u£^n, w; O) —)• J{un] O) 

as e 4 , 0. The result then follows from a diagonalization argument. □ 


5. Proof of the improved results 

5.1. Subcritical case 1 < p < d. In this section, we prove Corollary 2.9. We shall use truncations in 
place of the Sobolev compact embedding. For such truncation arguments to work, we need to restrict 
to the scalar case and to assume that the domain is fixed, i.e. domlL(?/, ■,u}) = domM for almost all 
y,uj. 


Proof of Corollary 2.9. In the proof of Theorem 2.2 and Corollary 2.4, the Sobolev compact embed¬ 
ding into bounded functions is used both in Step 2 of the proof of Proposition 3.10 and in Step 1 
of the proof of Corollary 2.4(i) (see Section 3.6). We only display the argument for Proposition 3.10 
(the argument for Corollary 2.4(i) is similar). 

We use the notation of Step 2 of the proof of Proposition 3.10. For all s > 0, define the truncation 
map Ts : M —>■ M as follows: 


Ts{x) = sign(a:)|xl A s = < 


s, if X > s; 

X, if —s < X < s; 

—s, if X < — s; 


and for all s > 0 consider the following s-truncation of u^^K,r,p- 

'^e,K,r,rj,s ■— T S IT ’^(Oj 1 ^)- 

Since \tUe,K,r,r],s — rt| < s -|- |fUK,r — u\, we may replace (3.29) by 

lim lim sup lim sup lim sup lim sup lim sup ||trte,K,r,»?,s ~ 'w||l°°(o) =0. 
r4,o ^40 s40 sf-O 


(5.1) 


(5.2) 


e,K,r,r],s — '^K,r')^'^e,K,r,ri (1 (l^£,K,r,r; 'WK,r))^I^K,r; deduce by COnvexity, 


Since Vu, 
noting that takes values in [0,1], 


Je{f''^e,K,,r,ri,si^iO^ ^ {V/^jt^'^e,K,r,rii^')dy (h'3) 

Jo 

+ / (1 - ri(«£,K,r,r,(y) - ?^K,r(y)))f^(y/e,tVM«,,,.,a;)dy 

Jo 

< Je{tUe,K,r,ri,i0;O) + \{y G O : \Ue,K.,r,rj{y) - Ui^,r{y)\ > s}| max M{tKi) 

l<l<k 

^ ^ 1 O') S ^||ri£^K,r,77 llLf’(O) ^)-^(^) T ^mm^ iVL(A;)^ . 

Since by definition (and by the Rellich-Kondrachov theorem) we have rte,K,r ,»7 (strongly) in 

L^(0) as e 4, 0, since the Afs all belong to domP, and since by assumption domM = domP, we 
deduce, combining this with (3.32), that 

lim lim sup lim sup lim sup lim sup lim sup Je{tUe,K,r,r),s-,^'-, O) < J{u; O). 
pi r4,0 k 4,0 774-0 sfO e4-0 


The rest of the proof is unchanged. 


□ 
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5.2. Minimal soft buffer zone for Dirichlet boundary data. In this section, we prove Corol¬ 
lary 2.10. In view of the error term (3.35) in the proof of Corollary 2.4, it seems that the speed of 
convergence of r/ to 0 with respect to £ must depend quantitatively on the speed of convergence of 
to 0 in In the case when the target function is affine (A • x, say), then ■= £^k{-/£, •) 

is the rescaling of the corrector and its convergence to zero is strictly related to the sublinearity of 
(^A at infinity, cf. Lemma 3.4. Even in the linear scalar case when V{y,A) = A • A{y)A for some 
matrix random field A, this sublinear growth can be arbitrary, and we expect that for all 7 < 1 , there 
exists a field A such that E[|(/?A(a^)P]^ ~ \xy as lx| S> 1 (see recent results in [28] and the example 
of Gaussian fields with non-integrable covariance). Yet, if instead of using the corrector (^a itself, 
which is in general not stationary and well-behaved, we use a proxy that is stationary, then the size 
of the buffer zone can be (optimally) reduced, at least for affine target functions, as the following 
proposition shows. 

Proposition 5.1. If for all A G yj^ have 

V(A)= inf E[V{0,A + Df),-)], (5.4) 

then the conclusion of Corollary 2.10 holds (and we can further replace Oe by any sequence % 0 

satisfying liminf^ ij^fe > 0). □ 

Identity (5.4) is essentially a regularity statement on quasi-minimizers of / 1 —)• E[l/(0, A -|- /, •)] on 
Lpot(n)™. By Poincare’s inequality, periodic gradients with mean-value zero are gradients of periodic 
functions, and hence in that case the space coincide with {Dcf : f G LPr’^(12)™'}, so that (5.4) 

is trivially satisfied. This already proves Corollary 2.10 under the additional assumption (1). 

On the other hand, the following result shows that (5.4) is also satisfied in the scalar case m = 1 
if the domain of V is fixed, in which case truncations are available. This proves Corollary 2.10 under 
the additional assumption ( 2 ). 

Lemma 5.2. If m = 1 and i/domE(y, -, 00 ) = domM is open for almost all y,oj, then (5.4) holds 
true for all A. □ 

Proof of Proposition 5.1. For all A G by assumption (5.4), for all 5 > 0, there exists a sta¬ 

tionary random field S LL^’^(f2; M™) such that E[Y(0, A -|- V(/?a,( 5 ( 0 , •), •)] < V{A) -|- 6. Set 
^A,(5,a; jg^ gy stationarity, for almost all a;, the Birkhoff-Khinchin ergodic theorem 

asserts that, for any bounded domain O C 

hm [ = |0|E[|99A,d"], Hm [ = |0|E[|V<^A,5r]. (5.5) 

etO Jq etO Jq 

Let O C be some fixed bounded domain. For 77 > 0, set Oy := {x £ O : dist(x, dO) > 77 }. For any 
sequence % i 0, (5.5) yields 

lim f \u^’^’‘^/£\P = 0 = lim / (5.6) 

eiO Jo\Or,^ elO JO\Or,e 

Fix such a sequence 77 ^ f 0. As in Step 1 of the proof of Proposition 3.10, for all A G 
we obtain the following, for some subset Oa C of maximal probability: for all w G Ha and all 
> 0, there is a sequence C W^’P{0]'BF) such that ^ 0 in W^'P{0]W^) as e | 0, 

limsup^ Je{A ■ X + ,uj] O) < J(A ■ x-,0) + 5, and such that (5.6) is satisfied. 

Let A and u G Oa be fixed, and let be as above. For all e > 0, choose a smooth cut-off 

function Xs with values in [0, 1], equal to 1 on = {x G O : d{x,dO) > rje}^ vanishing outside O, 
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and with jV^el < Cji]^ for some constant C'. Defining := € VDq^’^( 0; M”*), we obtain 

Jo\Or,e 

< J,(A • x + O) + 3P-^|A|?’|0 \ O^J 

Jo\Or„ 

and hence, if the sequence % | 0 is further chosen such that liminfg rj^/e > 0, 

limsup Je'"(A • X + ,uj] O) < limsup Je(A • x + ,oj] O) < J(A ■ x;0) + 6. 
etO e-fO 

Therefore, lim sup^ lim sup^ Jg ® (A • x + ,uj;0) < J(A • x;0). Combined with the T-liminf 


inequality of Proposition 3.2 and a diagonalization argument, this proves the result. 
Proof of Lemma 5.2. We split the proof into two steps. 

Step 1. Preliminary: we claim that it suffices to prove that, for all A G domD, 

limsup inf lElDlO,tA + Dd) •)] < 1/(A). 
tn <Aeiyi.r(L!) 


□ 


(5.7) 


Define V (A) := mf^g^i,p(Q) IE[D(0, A + Dcf, •)]. By definition V (A) > V{A) for all A, and hence 
property (5.7) together with the lower semicontinuity of V directly yields D(A) = lim^'i-i V (tA) for 
all A (and in particular the limit exists). Since v' is obviously convex, it is continuous on the interior 
of its domain. Since the domain is assumed to be open, this yields D(A) = V'{A) for all A. 

Step 2. Proof of (5.7). 

Let A G domf/ be fixed. Lemma 3.4 gives a measurable corrector u := ipA ^ Mes(D; VLjq^(M'^)) 
such that Vu G and D(A) = E[l/(0, A + Vu(0, •), •)]. For all i? > r > 0, choose a smooth 

cut-off function XR,r taking values in [0,1], equal to 1 on vanishing outside Qn and satisfying 

|Vx_R,r| < ‘Ifr. Also recall the definition (5.1) of the truncation Tg. We then set 

Ur{x,Uj) = u{x,uj) - f u{-,U!), Vr^{x,0j) = XR,rix)TsUR{x,Uj), 

Jqr 


and 


WR^r{x,uj) =j v% {x + y,TyU})dy=f v% {x + y,TyU))dy. 

\Qr\ JRd J-x+Qr 


Clearly, is well-defined, stationary, and belongs to LF^’^(D), with 

^w%r{x,u}) = 4 - Vvji,.{x + y,TyUj)dy. 

J-x+Qr 

Let t G [0,1). By Jensen’s inequality, 


Kl,^,it) := E[D(0,tA + tVn;|j,,(0, •),•)] =E 

< E 

and hence, by stationarity and the Fubini theorem, 


v(o,tA + t'f Vv%{y,Ty-)dy, 

V JQr 


V (0, tA + tVvR^^iy, Ty-), ■)dy 


-JQr 


KfiAt) < E 


-JQr 


V{y,tA + tVvR^^{y,-)r)dy 
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Decomposing 

tA + uj) =tA + txR,r{y)'^TsUR{y, w) + (1 - t)Y^VxR,r{y)'^sUR{y, to) 

= tXR,r{y)T^{uR{y,uj)){A + Vu{y,uj)) + t{l - XR,r{y))T^{uR{y,u}))A 
+ ^(1 - T'siuRiy, w)))A + (1 - i) uj), 

with Tg taking values in [0,1], we may then bound by convexity 


KkAt) < 


/ y{yA + ^u{yr)-,-)dy 
-JQr 


+ (1 - t)E%.^{t) 


+ M(A)/ {l-XR,r) + M{A)E 
J Qr 


/ {l-T'iuR{y,-))) 

-JQr 


(5.8) 


where the error term reads 

^R,r(i)=E 

By stationarity of Vu, note that 


M 


-J Qr 


l-t 


'yXR,r{y)TsUR{y,uj) dy 


E 


V{y,A + Vu{y,-),-)dy 


Q Qr 

For the error term, note that 


= E[D(0, A+ Vn(0, •),•)] =nA). 


1 - t 


'yXR,ri-)'ysUR{-,U}) 


< 


2t s 


L°°(0) 


1 — t r 


(5.9) 


(5.10) 


It remains to treat the last two terms of (5.8). Noting that fQ^(i — XR,r) = R < dr/R 


we obtain 


/ (l-XR,r)+E {I-Tg{uR{y,-))) < —+E ^\uii{y)\>sdy 
J Qr \-d Qr \ ^ VJ Qr 


dr , _ 

<^+; p 


R Qr 

1 

R 


/ P ^ u{Ry,-) - 'f u{Rz,-)d. 
Jq [R Jq 


> 


R 


(5.11) 


dy. 


Lemma 3.4 (together with the Rellich-Kondrachov theorem) gives ^\u{R-,uj) — j-Qu{Rz,uj)dz\ —)• 0 
(strongly) in L^(Q) as R '[ oo, for almost all uj. Hence up to an extraction in R (implicit in the 
sequel) we deduce that, for almost all y ^ Q, ^\u{Ry,-) — j-Qu{Rz, ■)dz\ —)• 0 almost surely. Since 
almost sure convergence implies convergence in probability, we deduce by dominated convergence, for 
all e > 0, 


lim 

^toojQ 


R 


u{Ry, •) “ 7 u{Rz, ■)dz 


Q 


> e 


dy = 0. 


A diagonalization argument then gives a sequence e_R i 0 such that 

1 


lim / P 


R 


u{Ry, •) “ 7 u{R^^ ')d' 

Jq 


> £r 


dy = 0. 


(5.12) 


Choose s = Sr := Rer and r = vr := R^/Fr. By assumption, there exists some 5 > 0 with 
adhi ?5 C intdomM. By (5.10), for all t S [0,1), there is some Rt > 0 such that for all R > Rt 


l-t 


,rR {■)TsrUr{-,uj] 


2t 


L°°(0) 




(5.13) 


Combining this with (5.8), (5.9), (5.11), (5.12), and noting that M(A) < oo follows from the choice 
A G domH, we obtain 

lim sup lim sup (t) < V (A), 
m Rtoo ’ ^ 
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and the result (5.7) follows. 


□ 


5.3. Approximation by periodization in law. This last section is devoted to the proof of Corol¬ 
lary 2.13. For that purpose, we first prove the following version of the r-limsup inequality of Propo¬ 
sition 3.10 in expectation. As Lemma 3.4 only gives almost sure (and not L^(r2)) control of the 
sublinearity of the corrector, we need to use truncations to prove this result (and therefore restrict to 
the scalar case with fixed domain). 

Proposition 5.3 (F-limsup inequality in expectation with Neumann boundary data). Let V, Jg, J, M 
be as in Theorem 2.2 for some p > 1. Also assume that we have m = 1 and doml7(y, •, a;) = 
domM for almost all y,io. Then, for all bounded Lipschitz domains O C and all u £ 
there exists a sequence {u^)^ C Mes(n; VF^’^(O)) such that E[j£(n£, •; O)] —>• J{u;0) and ^ u in 
L“(L1;L°°(0)). □ 

Proof. We split the proof into three steps. 


Step 1. Preliminary. 

We claim that, for all bounded domains O C and all A £ int domP, there exists a se¬ 
quence {uA,e)e C Mes(fl; VF^’^(O)) with UA,e{',^) —^ A • x weakly in for almost all uj and 

E[j£(uA,e, •) O)] ■ x]0') for all subdomains O' C O. This is indeed a simple reformulation of 

Lemma 3.4 with the notation u\^cix,oj) := e^p/£, iw). 

Step 2. Recovery sequence for continuous piecewise affine functions. 

Let O C be a bounded Lipschitz domain and let u be a continuous piecewise affine function on O 
such that Vn £ intdomP pointwise. We prove the existence of a sequence (ue)^ C Mes(r2; IT^’^(O)) 
with Us -i- u in L°°(n; L°°(0)) and E[Js{ue, •; O)] —)■ J(rt; O). 

Let us modify u as in Step 2 of the proof of Proposition 3.10 to make its variations smoother (using 
Proposition A.17), and use the same notation. By Step 1, for all 1 < i < nK,r) there exists a sequence 
('^£,K,r)e Mes(n; iPjof (E'^)) with n® ,, ,.(•, w) —^ c).-|- A).• x in lTjQ’|f(E'^) for almost all lo and such 
that, for all Lipschitz subdomains O' C O, we have 'E[Js{ul ^ •; O')] —)• T(A). ^ • x; O'). Consider the 

same partition of unity ™ Step 2 of the proof of Proposition 3.10, and also recall the 

definition (5.1) of the truncation Tg. We now set, for s > 0, 




:= U 


K,r 


+ 


E 

2=1 


+ A),• x)) 


£ Mes(L!;W^’P(0)). 


On the one hand, since \tUs,K,r,ri ~ u\ < s -|- |tuK,r ~ u\ pointwise, we deduce 

limlimsup limsup limsuplimsup limsup \\tUe,K,r,ri — 'w||l°°(iI;L°°(o)) = 0- (5-14) 

r4.0 f;10 s4-0 sfO 

On the other hand, since 


tVu, 


£,K,r,r],s 


2 = 1 


'^tXn,r,v'^sK,K,r - «,r + K,r ' 

+ ^XK.r.yi^ ~ ^s('^e,K,r ~ ('^K,r + ^K,,r ' 4" (^ ~ t)Se,K,r,T],s,t, 


2=1 


where we have set 


Se,K,r,ri,s,t ^ i^K,r + ' ^))^XK,r,r] 4" ^ XK,r,r](^''^K,r ^K,r)i 


2=1 


2=1 
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we deduce by convexity 

, r 

¥.[JSue,K,r,r,,s, S O)] < (1 “ •; 0^+ 


(5.15) 


2 = 1 


+ ^M(A*^,.)E[|{y € : K ,,,.(y,w) - (c),,. + ,, • y)\ > s}l] 


where the error reads 


2=1 


Ee,K,r,r],s,t — E 


'O 


^{yISe,K,r,r],s,t{y^ ')) ')dy 


For all i, set A^K,r -,?7 •= {j • i / b Oi^r,ri n OK^r,r; / ■ We estimate S'e,K,r,r),s ,4 for all y,uj as follows: 

t Cn t 

\Se,K,r,7^,s,tiy^^)\ ^ -1-^S + Tj-7 SUp SUp |A^ - A* |, 

and hence, since limsup^j^g suPjeAf* \-^K,r — AJ. ,.1 < k for all i, 

limsuplimsuplimsuplimsup ||5e,K,r,»?,s,t|lL°°(0;L°°(0)) =0- (5.16) 

mJ'O rfl-O s 4,0 

The last term of the right-hand side of (5.15) is estimated by 

E[|{y G : \ul^^^^{y,uj) - ■ y)\ > s}|] < [ P[K^^^^(y, w) - ■ y)\ > s]dy. 

Jo 

For almost all w, the Rellich-Kondrachov theorem shows that u* ^ ^(-j w) ^ cj. ^ -|- A). ^ • x (strongly) 
in L^(0), and hence, up to an extraction in e (implicit in the sequel), this convergence holds almost 
surely, almost everywhere on O. The dominated convergence theorem then yields 

limsupE[|{y G 0^+,, ; K_^,.(y,a;) - + ■ y)\ > s}|] = 0. (5.17) 

eiO 

Finally we note that by construction, 

EiK,,r f^K,r 

lim^E[J,(n^„.; 0 *+^)] = ^ J(A),,, • x; 0 *+,) = | 0 *+,|F(A),,,). 


2=1 


2 = 1 


Hence, arguing as in Step 2 of the proof of Proposition 3.10, combining this with (5.15), (5.16) 
and (5.17), and recalling that domH = domM, we obtain 

limlimsuplimsuplimsuplimsuplimsupE[Je(tUe^K,r,?7,s) •; O)] < d{u]0), (5.18) 

rtO K^O f;4-0 stO £4-0 

so that the conclusion follows as in Step 2 of the proof of Proposition 3.10 (note that by Fatou’s 
lemma the F-liminf inequality also holds in expectation). 

Step 3. Conclusion. The existence of recovery sequences for general functions can now be deduced as 
in Step 3 of the proof of Proposition 3.10, using the result of Step 2 above. □ 


We now turn to the proof of Corollary 2.13 itself. In case (1), we use Proposition 5.3 to control 
the energy close to the boundary of the cube, where the periodization in law may have modihed 
the integrand (which however is the same in law locally and can therefore be handled by taking the 
expectation). In case (2), we use the particular geometric structure to explicitly estimate the energy 
density near the boundary and prove its uniform integrability. 

Proof of Corollary 2.13. Let (I7 ^)_r>o be an admissible periodization in law for V in the sense of 
Definition 2.11, and for all e > 0 and all A G denote by J|’®’^(-, A, •) the following random 

integral functional on the unit cube Q = [— 5 , 5 )“^: 

JP®"('u,A,w) := [ V^^%y/e,A + Vu{y),u))dy, 

Jq 


uG Wp'’?(Q;M"*). 
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We split the proof into three main steps. 

Step 1. F-liminf inequality. The results of this step hold without the additional assumptions (1) 
and (2). Let n" C be a subset of maximal probability such that the almost sure property of 
Definition 2.11(ii) holds pointwise on D". We claim that for all u £ Wpe?(Q;M™), all Ug —^ u in 
fLper(Q) and all ^ A, we have liminfgj^o > J(u + A • x; Q) for all a; G D". 

For all 9 G (0,1) and lo G D", we thus have for all e > 0 small enough 

= Jeiue + ■ x,uj-,Qe) + [ /e, A^ + Vueiv), u!)dy. 

jQ\Qe 

By the F-liminf inequality for (see Proposition 3.2) and by the non-negativity of 

this turns into 

liminf Ag,w) > J{u + A - x^Qg), 

e\.0 

and the result follows by the arbitrariness of 0 < 1 and the monotone convergence theorem. 

Step 2. P-limsup inequality under assumption (1). We assume m = 1 and domV{y,-,oj) = domM 
for almost all y, 0 J. We prove the existence of a subset Vt'" C D" of maximal probability such that, 
for all u G bFper(Q), all A G and all ui G D'", there exists a sequence {us)s C Mes(D; Wpe?(Q)) 

such that Ug —>■ 0 in L°°(D; L°°(Q)) and 

limliminf J^^'^{tUe{-,uj) + tu, tA,uj) = limlimsup {tUe{--,oj) + tu, tA,Lo) = J{u + A - x;Q). 
tn £4.0 tn e^o 

Moreover, the limits f f 1 can be dropped if we have J{au + aA ■ x]Q) < oo for some a > 1. We split 
the proof of this step into three parts. 


Step 2.1. Local recovery sequence for affine functions. 
For all i? > 0, define the event 


{w G D : W(-, •, L^ 

for some fixed 6q > 0. Definition 2.11(ii) ensures that P[D/j] f 1 as i? f oo. Let A G let O C Q 

be a domain of diameter smaller than 6o/2, and let R > 0. By the diameter condition, there exists 
xq £ Q such that O C xq + 9oQ. We claim that there is a sequence {ue)^ C Mes(D; IF^’^(O)) such 
that u, ^ 0 in L°°(D; L“(0)) and Jo A + Vu^, •)] ^ P[f^R] J(A • x; O). 

By definition, for all w G Dr we have •, w)|q l^ < 1/.^- 

^uQ/e ^oQ/e 

Therefore, applying Proposition 5.3 on —xq + O restricted to Dr yields the following: there exists a 
sequence {uc)e C Mes(D; xq + O)) with ttg —)• 0 in L°°(D; L°°(—xq + O)) and 


limE 


[ DV-(-/e,A + Vn„-) 

= limE 


£4-0 

L j 

-xo+O J 

£4-0 

J—xo+O 


= P[Dr] J(A • x; -xo + 0)= P[Dr] J(A • x; O). 


Now define := u^{—xq + £ Mes(D; xq + O)) for all e. By construction —)• 0 in 

L°°(D; L°°(0)), and the stationarity property of Definition 2.11(i) further gives 


limE 

£4-0 

= limE 
£4-0 

= limE 
£4-0 

= limE 
£4-0 


^ Jo 

lOfl V^/‘^{yle, A + Vue{-xo + y, ■),Tl^JjOdy 

lOfl / P^/^((-xo + y)/e, A + Vue{-xo + y, •), Jdy 

Ififl / V^/%y/e,A + Vue{y,-),-)dy 

J-xn+O 


= P[Dr]J(A-x;0). 


Step 2.2. Global recovery sequence for (piecewise) affine functions. 
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Let A G and let u be a continuous piecewise affine Q-periodic function on with A + Vtt G 
intdomf/. We claim that there exists a sequence (ue)e C Mes(n; Wpe?(Q)) such that ^ u m 
L°°(Q)), and, for almost all cu, 

limliminf w) + tu,tA,uj) 

tfi efo 

= limlimsup w) + tu, tA, w) = J(ri + A • x; Q). (5.19) 

tn £4.0 


Consider the partition of Q associated with u, and let Q = 1+jf^^ Oi be a refined partition such that 
all the Ois have diameter at most 9q/2. For all I, define ci+ Ai ■ x := Vtt|o; with A + A; G int domf/, 
and choose xi € Q such that Oi G xi + OqQ- Given R > 0, define 

k 


Xe,R 


n 

1=1 


xi/e 




Step 2.1 then gives a sequence {u\)^ C Mes(fl; fy^’^(0;)) such that —)• 0 in L°°(n; L°°(0;)) and 


lim sup E 
etO 


Xe,R [ V^^^-/e,A + Ai+Vui,-) <limElve„ [ V^/%-/e,A + Ai + Vui, ■) 

Joi J =10 L ^Joi 


= ¥[i}n]J{u + A-x-,Oi). 


We can then repeat the argument of Step 2 of the proof of Proposition 3.10, and glue the recovery 
sequences for the (small) affine parts, using this time Proposition A. 18 instead of Proposition A. 17 
(applied to the function u). This allows to deduce the following: there exists a sequence {us)e C 
Mes(fl; bFpe?(0)) such that Ug —)• 0 in L°°(n; L°°(Q)) and 

limsuplimsupE[x£,Rj£'^'^(ttf£ + tir, tA, •)] < J(ri + A • x; Q). (5.20) 

4tl efo 

Now the P-liminf inequality ofStep 1 yields liminftliminfe J|’®’^( lie (•,t<;) +A, w) > J{u + A-x]Q) 
for almost all uj. Assume by contradiction that this inequality is strict for all uj in some set of positive 
probability (even for some subsequence in £,t). In this case we would find £o,6 > 0 and to ^ (0,1) 
such that > 0, where we have defined the event 

nfo,io := S {tue{-,uj) + tu,tA,uj) > (5 + J{u + A • x; Q), for all e < eq, t > to}, 

By definition we have 


lim inf liminf E[x£ + tit,tA, •)] 

ifl efo ’ 

> J(ii + A • x; Q) liminf E[x£ r] + liminf liminf E[x£ + tu,tA, •) — J{u + A • x; Q))] 

e4-0 ’ efo ’ 

> J(it + A • x; Q) liminf E[x£ r] + (51iminf EIIqs Xe r] 

eto ’ efo eo.‘o ’ 

+ liminf liminf E[lQ\oi Xe RiJf’^^itu^ + tit,tA, •) — J{u + A • x; Q))]. 

Ifl efO ' ’ 


Now a union bound argument yields 

E[x£,r] > 1 - kF[n \ L!fi], E[1^S Xe,R] > P[<,io] " \ 

^0 ’^0 

and hence, by Fatou’s lemma and by the F-liminf inequality of Step 1, we obtain 


liminf liminf E[x£ RJf^^{tUe + tu, tA, •)] 

tfi eto ’ 

> J(ii + A • x; Q)(l - km \ ftjj]) + \ ^r])- 

However, taking the limit R '[ oo, and recalling P[H/j] '[ 1, this contradicts (5.20). The result (5.19) 
is then proven. 


Step 2.3. Recovery sequence for general target functions. 
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For all A £ intdomy, let C be a subset of maximal probability such that equalities (5.19) 
hold for all u £ for the choice u = 0. We then define il.'" C as the (countable) intersection 
of all Q\s with A £ n int doml/. Repeating on this basis the gluing argument of Step 2 of the 
proof of Proposition 3.10, for all A £ and for all continuous piecewise affine Q-periodic functions 
u on with A + Vu £ n intdomR, we may prove that (5.19) holds for all ui £ 11'". 

Let now u £ WpS{Q) and let A £ By Step 1, we may assume J{u + A- x; Q) < oo, so that by 

convexity A £ domR. By convexity and lower semicontinuity of V, lim^ii J (tu + tA - x',Q) = J{u + A- 
x',Q). For all t £ (0,1), since A £ domR, we have 0 £ intdoml/(- + fA) and Proposition A.16(ii)(b) 
then gives a sequence {un,t)n of continuous piecewise affine Q-periodic functions such that A + A/un^t £ 
Qixdnint domR, Un,t tu in and J{un,t + tA-x] Q) —)• J{tu + tA-x] Q). We can then apply 

the result for the (rational) continuous piecewise affine approximations on H'", and the conclusion 
then follows from a diagonalization argument. 

In the case when J(au + aA ■ x]Q) < oo for some a > 1, the limits f f 1 can be dropped. Similarly 
as in the proof of Corollary 2.4 (see Sections 3.6 and 3.7), it is enough to replace u,A by u/t,A/t for 
f £ [1/a, 1) and to adapt Steps 2.2 and 2.3 above accordingly. We omit the details. 


Step 3. F-limsup inequality under assumption (2). We assume p > d, and V{y,A,co) <C{1 + [Aj^) 
for all uj and all y ^ E{co), for some random stationary set satisfying almost 

surely, for all n, for some constant C > 0, 


25 " := 


1 

J- 


inf dist{BR^{q^),BR^{q^))> 

m,m^n 


1 

c’ 


Rn<C. 


(5.21) 


We prove the existence of a subset Q,'" C Al" of maximal probability such that, for all u £ WpS{Q), 
all A £ and all oj £ Al'", there exists a sequence ^ 0 in LFpe?(Q;M™') such that 


limliminf + tu, tA,uj) = limlimsup J^^'^{tu‘^ + tu, tA,Lo) = J{u + A ■ x;Q). 

m £ 4.0 tn £ 4.0 


Moreover, the limits f f 1 can be dropped if we have J(au + aA ■ x;Q) < oo for some a > 1. Finally, 
we have for all A £ and almost all w 

’F(A) = lim inf f {y/e, A+ Vu{y),uj)dy, (5.22) 

elO tiGVUphflQ;®"*) Jq 

where the convergence also holds in expectation. We split the proof of this step into three parts. 

Step 3.1. Recovery sequence for affine functions. We prove that for all A £ for almost all oj 

there is a sequence u" —^ 0 in LFpe?(Q;M™) such that Jg®’'(A, u", w) —)■ J(A • x;Q). 

Let A £ be fixed. Lemma 3.4 and the Birkhoff-Khinchin ergodic theorem give a random field 

qjR £ Mes(n;lF[Q/’(M'^;M'")) such that, for almost all oj, eip\{-/e,Lj) — ^ 0 weakly in 1F[q^(M'^; M'") and 
also J£(A • X + S(p\{-je, oj),oj; O) —)• J(A • x; O), for all bounded domains O C Let oj £ fl" be fixed 
such that both properties hold, and set u" := A ■ x + £yj\{-/£,oj). 

For all e > 0, consider the Q/e-periodic random set E^ = U^i .^1?“„ (?£ n) associated with the 
periodization in law at scale 1/e. By stationarity, these periodized inclusions also satisfy (5.21), 
that is 25"^ > l/C and < C almost surely. For all £,n, define := BR<^^{q‘^,^) and : = 

ER'^^{i+iA5^„){Qe,n)- ^y definition, the are all disjoint. For all e,n, choose a smooth cut-off 

function Xen that equals 1 on vanishes outside e.B"„ and satisfies | Vx"„[ < 2/(ei?"„(lA5"„)) < 

4C'/(ei?"„). We choose these cut-off functions in such a way that x"n i® Q-periodic. 

Let 6 £ (0,1). Denote by (resp. M/') the set of all n > 1 such that {Q \ 9Q) n ^ 0 

(resp. dQ n £B^^ 0)- Choose a smooth cut-off function xe that equals 1 on 6Q, vanishes outside 

Q and satisfies [Vx^j < 2/(1 — 6). Set 

Xe,e{y) ■= xeiv) - xe{y) xtniy) + Y1 A£,n(y)/_ {xe{z) - xe{y))dz. 

J c- Ou; 

n£M^ nON^g\M^ 
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This defines a smooth cut-off function 0 that equals 1 on {6 — 2C£)Q, vanishes outside Q, is constant 

on each inclusion and satisfies < C/{1 —9) for some constant C > 0. We then define 


‘‘efi 


(y) = (l-X£, 0 (y))A- 2 / + (l-X£, 0 ( 2 /)) X] A£,n(y)/. A-{z-y)d; 

neN^„ -te-Be.n 

e,u 

+ xZm(i- E xUv))<iv) + <>(«) E xUv){- 

V n£N^ ^ J sB] 


n^N't 


By construction G A • x -|- lTpe?(Q;M™'). Since at least on {0 — 2eC)Q, and since 

Vttg g = 0 on all eB'^^ with n G we have, for all e > 0 small enough, 


[ yV.(-/e,V<g,a;)< [ 

JQ 


Vi-/s,Vu‘^,oj) 

M{0)\Q\{e-2sC)Q\+C [ 

JO' 


{9-2eC)Q 

+ 


lQ\[e-2eC)Q 

< j£«,w;Q) + C(l-0 + 2eC) + C 


(1 + iVn-gD 


lQ\{e-2eC)Q 

We need to estimate the last term. A straightforward calculation yields 




(5.23) 


L 


Q\{e-2eC)Q 


IVu"*!” < C(1 - 0 + 25C)|A|>’ + ^ !„-(„) _ A. 

+ C E ((1/- A. ( 2 - 9 )* 

nGM“ 

+ C E ((1-9)-''+(£<„)-”) [_ <-/, 

J SB'£ rr, J SB^ rr. 


Q\(e-2eC)Q 
P 

dy 


|VA 


CJ |p 




and hence, by the Poincare inequality and the estimate \^^e,n\ ^ IQ \ (1 ~ 2eC')Q| < 2d£C, 

C 


I 

Jc 


Q\(B-2eC)Q 


|V<g|P <C{l-e + 2eC)|A|P + 


{i-OY 


+ C 1 + 


L 


(1 - 9)Pj jQ\[e- 2 eC)Q 
For all £0 > 0, the Birkhoff-Khinchin ergodic theorem yields, for almost all uj 

[ (|A|^ + |V(^A(-/6,a;)n<lim f 

jQ\(d-2eC)Q efO JQ' 


[ \£^A{-/£,Lor (5.24) 

Jq 

(|A|^ + |V^A(./e,a;)n. 


lim sup 

etO J Q\{9—2eC)Q 


{\k\P + \V^A{-le,u:)?) 

^ Q\{0—‘^sqC)Q 

< |Q\(0-2eoC)Q|(|Ar + E[|Vv;An) 

<C{l-e + 2eoC)(|A|P + EOVi^aII). 

Combined with (5.24) and the sublinearity of the corrector in the form of £^pa{'I£) —)• 0 in LP((5; M™) 

(see Lemma 3.4 and the Rellich-Kondrachov theorem), this yields, for almost all uj, 


lim sup 


L 


|V<g|P < C(1 - 0)|A|P + C{l-9 + 2eoC')E[|A + Vc/paH. 

etO JQ\{e-2eC)Q 

We then insert this estimate in (5.23) and to pass to the limits 9 1 and £9 i 0 to obtain for almost 

all oj 


[ V^^^{-/£,Vu‘^ 0 ,oj) < lim Je(tt£,w;Q) = y(A) = J(A • x;Q). 

Jq ’ etO 


lim sup lim sup 
efi £4-0 JQ 

The conclusion then follows from Step 1 and a diagonalization argument. 
Step 3.2. Recovery sequence for general target functions. 
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As we assume here p > d,we may use the gluing argument of Step 2 of the proof of Proposition 3.10 
to pass from affine to piecewise affine functions. The conclusion follows by approximation similarly 
as in Step 2.3 above. We omit the details. 


Step 3.3. Proof that the almost sure convergence (5.22) also holds in expectation. 

Using the same notation and cut-off functions as in Step 3.1, we define € Vlper(f5; 1^”*) as 


u, 


’(y) ■='^Xe,niy)'f _ A-iz-y)dz. 


Testing the infimum in (5.22) with u = u^, we obtain by similar calculations as in (5.23) and (5.24) 
that for all e > 0 and all oo 


inf 


[ V^/^{-/e,A + Vu,u}) 

Jq 

<c{i + \Ar)+cY,i^Rip-^ [ 

„ jQi 


f . 

f A-{z- y)dz 




dy<Cil + \An. 


The conclusion then follows by dominated convergence. 


□ 
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A. Appendix 

A.l. Normal random integrands. In this appendix, we briefly recall the precise dehnition of 
normal random integrands (as defined e.g. in [23, Section VIII.1.3]) and we prove their main properties, 
mentioned in the beginning of Section 2 and used throughout this paper. Let (fl, P) be a complete 
probability space. We denote by the (not completed) Borel cr-algebra on 

Definition A.l. A normal random integrand is a map VU : X x fl ^ [0, oo] such that 

(a) W is jointly measurable (i.e. with respect to the completion of i3(M'^) x i3(M™'^'^) x F)] 

(b) for almost all cu, there exists a map V),; : x —)• [0, oo] that is x .S(M™'^'^)-measurable 

and such that W{y, -, 0 ;) = I4;(y, •) for almost all y; 

(c) for almost all y, there exists a map Vy : x fl —)• [0, oo] that is .S(M™^'’*) x J^-measurable and 

such that W{y, -^oj) = Vy{-,uj) for almost all cu; 

(d) for almost all y, 0 J, the map W{y, -jU;) is lower semicontinuous on 

It is said to be r-stationary if it satisfies (2.1) for all A,y,z,uj. □ 

As shown e.g. in [23, Section VIII. 1.3], a simple example of normal random integrands is given by 
the so-called Caratheodory random integrands, that are maps LU : M'’* x x U —[0,oo] such that 

W{y, -jOj) is continuous on for almost all y, oj, and such that W (•, A, •) is jointly measurable on 

X n for all A. 

As already advertised in the beginning of Section 2, the reason for these technical assumptions is 
that they are particularly weak but still guarantee the following properties: 

Lemma A. 2. Let W : ^ ^mxd ^ 

fl —>■ [0, oo] be a normal random integrand. Then, 

(i) for almost all io, the map y i—)• W{y,u{y),io) is measurable for all u € Mes(M'’*, 

(ii) for almost all y, the map lo i—)• W{y, u{io), lo) is measurable for all u S Mes(fl, M™'^'’*). □ 




52 


M. DUERINCKX AND A. GLORIA 


Proof. For almost all ui, part (b) of Definition A.l gives a x i3(M™^'^)-measurable map 

on X such that W{y,-,io) = K;(y, •) for almost all y. Hence, for u G Mes(M‘^; the 

map y i—>• W{y,u{y),uj) is equal almost everywhere to the map y i—)• Vu){y,u{y)), which is necessarily 
measurable since Id xrt : X —?> x is measurable and since is Borel-measurable. This 

proves (i), and (ii) is similar. □ 

If W is r-stationary, we may write IF(y, A,a;) = W{0, A,T-yUj), which thus receives a pointwise 
meaning in the first variable, and Lemma A.2(ii) may obviously be strengthened as follows. 

Lemma A.3. Let W : X ^mxd X 

Ll [0, oo] be a r-stationary normal random integrand. For all 
y and all u G Mes(D; the map u) i—>■ W{yoj) is measurable. □ 

A.2. Stationary differential calculus in probability. In this appendix, we precisely define the 
measurable action r that is used throughout this paper to induce the stationarity, and we discuss the 
properties of the stationary derivatives defined in Section 3.1.1, and prove in particular the useful 
identity (3.3). 

A.2.1. Stationary random fields. As usual, the standard notion of stationarity of random fields (de¬ 
fined as the translation invariance of all the finite-dimensional distributions) is strictly equivalent to 
a formulation of stationarity as the invariance under some (measure-preserving) action of the group 
of translations (M*^,-b) on the probability space (see e.g. [32, Section 16.1]). This point of view is of 
great interest, since it puts us into the realm of ergodic theory. 

Because we focus on jointly measurable random fields, which is standard in stochastic homoge¬ 
nization theory (see also Remark A.6 below), a further measurability requirement is added in our 
definition of an action, as e.g. in [31, Section 7.1], 

Definition A.4. A measurable action of the group (M*^,-b) on (D, P) is a collection r := {Tx)x£m.'^ 
of measurable transformations of D such that 

(i) TxOTy = Tx+y for all x,y £ 

(ii) P[r 3 ;A] = P[A] for all x G and all A G 

(hi) the map x D —?■ D : (x, cu) i—>• is measurable. □ 

For any random variable / G Mes(D;M), we may define its r-stationary extension / : x D —M 

by f{x,uj) := /(r_a;Cu), which is a r-stationary random field on (in the sense that /(x y,uj) = 
/(x, T-yUj) for all X, y, ui) and which is by definition jointly measurable on x Q. For / G LP(D), 1 < 
p < oo, the r-stationary extension / belongs to L^(D; L[(j^(M'^)). In this way, we get a bijection between 
the random variables (resp. in L^(D)) and the r-stationary random fields (resp. in L^(D; L](j^(M'^))). 

We may also naturally consider the associated action T := {Tx)x^M.d, of (M'^,-b) on Mes(D;M), 
defined by {Txf){uj) = f(T-xUi) for all cj G D and / G Mes(D;M). Let 1 < p < oo. The following 
gives elementary properties of this action (see e.g. [31, Section 7.1]): 

Lemma A.5. The action T defined above is unitary and strongly continuous on L^(D). □ 

In the context of stochastic homogenization theory, the measurability hypotheses made above (as 
in e.g. [31, Section 7.1]) are sometimes replaced by stochastic continuity hypotheses (see e.g. [37, 
Section 2]). As the following remark shows, both are actually equivalent. 

Remark A.6 (Measurability or continuity). It should be noted that the additional measurability 
assumption (iii) in Definition A.4 above is not inoffensive at all. Indeed, a stochastic version of the 
Lusin theorem can easily be proven: a random field h on is jointly measurable if and only if, for 
almost all x G for all 5 > 0, 

lim P[|/i(x + y,uj) — h{x, cj)! > 5] = 0. 

y^O 

Hence, for a stationary random field h on joint measurability is actually equivalent to stochastic 
continuity (and even to continuity in the p-th mean, in the case when /i(0, •) G L^(D)). In the 
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same vein, the measurability property (iii) in Definition A.4 is equivalent to the strong continuity 
of the action T of (M'^,+) on L^(D), and also to the property that all r-stationary extensions are 
stochastically continuous. □ 

A.2.2. Stationary Sobolev spaces. Let 1 < p < oo, and let the stationary gradient D and the space 
be defined as in Section 3.1.1. Now we present another useful vision for derivatives of 
stationary random fields. 

Given a random variable / S L^{Q), the r-stationary extension is an element / S L^(D)) = 

L^(D; and can thus be seen as an L^(D)-valued distribution on M'^. We may then de¬ 

fine its distributional gradient V/ in the usual way. Note that by definition, for almost all uj, 
V/(-,a;) is nothing but the usual distributional gradient of G L|^^(]R'^). As usual, the Sobolev 

space L^(D)) is defined as the space of functions / G LP(D)) such that V/ G 

LP(D; and in that case V/ is called the weak gradient. The following result shows the 

link with stationary gradients and with the space in particular proving identity (3.3). 

Lemma A.7. Modulo the correspondence between random variables and r-stationary random fields, 
we have 

= {/ G : f{x + y,u) = /(x,r_,a;),Vx,y,a;}, 

and moreover V/ = Df for all f G VL^’^(D). □ 

Proof. Denote for simplicity Ep := {/ G L^(D)) : f{x-\-y,u}) = f{x,T-yUj),'\/x,y,uj}. For all 

i, the stationary derivative Dif is defined as the strong derivative of the map M —^ L^(D) : h i— T_h(,.f, 
so that 

= {/ G : f{x + y,u^) = /(x,T_,a;),Vx,y,w}. 

Hence, for all / G we have / G Ep, and V/ = Df, since weak derivatives are generalizations 

of strong derivatives. 

We now turn to the converse statement. Let / G VF^’^(D). As in [31, Section 7.2], choose a 
nonnegative even function p G C'^(]R‘^) with f p = 1 and suppp C Bi, write ps(x) = 5~^p{x/5) for 
all (5 > 0, and define a regularization G L^(D) by 

Rs[f]{^)= j Ps{y)f{y,^)dy, (A.l) 

JRd 

or equivalently, as p^ is even, 

RslfMx,^]) = j ps{y)fix + y,uj)dy = ps{y - x)f{y,uj)dy = {ps * f{-,u})){x). 

Clearly, Rs[f] —)• / in LP(D), and hence by stationarity Rs[f]{x,-) f{x,-) in L^(D) uniformly in 

X. As by definition Rs[f] G L^(D)), we have Rslf] G and the stationary gradient is 

simply DRs[f] = Rs[Vf]. This proves DRs[f] —)• V/ in L^(D;M'^), and thus DRs[f]{x, •) —)■ V/(x, •) 
in L^(D;M'^) uniformly in x. Hence Rs[f] —> / in L^(D)), so / G L^(D)), from which 

we conclude / G □ 

A.3. Ergodic Weyl decomposition. In this appendix, we discuss the various properties of the 
ergodic Weyl spaces recalled in Section 3.1.2. In particular, we prove the ergodic Weyl decomposi¬ 
tion (3.7) as well as the density result (3.8). 

Let r be a measurable action of {W^, -|-) on the probability space (0,^", P), let 1 < p < oo, and let 
LpQ^.(D), Lg^j(D), T’pq^(D) and be defined by (3.4) (or equivalently (3.5)) and (3.6). First we 

prove the ergodic Weyl decomposition (3.7): 

Proposition A.8 (Ergodic Weyl decomposition). Let r be ergodic and let 1 < p < oo. Then the 
following Banach direct sum decompositions hold: 

LP{n-, M-') = L^„,(L!) © Elfin) = Efi^fin) © LlfiQ) = Efi^fiQ) © Efijn) © M-'. 

□ 
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Proof. Let 1 < p < oo and let q := p' = p/{p — 1). The following weaker form of the result is proven 
e.g. in [31, Lemma 15.1]: 

= Lljn) = Fl,{n) © M'', (A.2) 

= F^^,in) © m '', 

in the sense of duality between L^(11;M'^) and L'^(n;]R'^). In the Hilbert case p = 2, this is already 
enough to conclude by the orthogonal decomposition theorem. For the general case, we will need to 
use more subtle (nonlinear) decomposition results in Banach spaces. 

As the Banach space L^(n; is uniformly smooth and uniformly convex, the following (nonlinear) 
direct sum decomposition (see e.g. [3, Theorem 2.13]) holds for any closed subspace M C L^(fl;M'^): ^ 

LP(fl;M'^) = M© J(M-^), 

where M"*- C L'^(H;M'^) denotes the orthogonal in the sense of duality, and where the (nonlinear) 
map J : L'^(n;M'^) —)• L^(fl;M'^) is defined by J{u) := Let us apply this result 

to the choice M = By (A.2) we have M-*- = Lg^j(n), and we may compute DJ{u) = 

{q — Du, so that JM-*- C Lg^j(n). On the other hand, given u G we may write 

u = J{v) with V := G Lg^j(fl) = M-*-. This proves equality = Lg^j(fl), and the 

result follows. □ 

The density result (3.8) is now easily obtained, using the same (nonlinear) direct sum decomposition 
in Banach spaces: 

Proposition A.9. Let r be ergodic, and let 1 < p < oo. Then, 

FP^,{n) = adhLP(n;R4){Ilx : X G 
FPJQ) = adhLP(f,;M.){ll X X : X G 

□ 

Proof. Using the same notation as in the proof of Proposition A.8, the (nonlinear) direct sum decom¬ 
position of [3, Theorem 2.13] gives 

LP(L!;M'^) = adh DW^’P{n) © J{{adh DW^’P{n))^). 

By definition (adh ZlLF^’^(ll))-*- = LgQj(n), while as in the proof of Proposition A.8 we find J(Lg^j(n)) = 
Lg^l(ll). This yields 

LP(L!;M'^) = adhT»H'°°(!d)©L^^j(L!). 

The obvious inclusion adh. DW^’^{D) C TpQ^(ll) and the direct sum decomposition of Proposition A.8 
then imply the equality adhllLU^’^(fl) = FpQ^(fl). The result for is obtained similarly. □ 

A.4. Measurability results. This appendix is concerned with various measurability properties. 

A.4.1. Measurable potentials for random fields. The following result complements the equivalent def¬ 
inition (3.5) of LpQ^(n), and shows that potentials associated with potential random fields may be 
chosen in a measurable way with respect to the alea. 

Proposition A.10. Let r be an ergodic measurable action on a complete probability space (11,F, P). 
Let 1 < p < oo. For all f G LpQ^(fl) there exists a random field (j) G Mes(H; VF[^’^(M'^)) such that 
f{-,uj) = V(/)(-,(x:) for almost all oo. □ 

1. The direct sum means here that any element of L^(P; R'*) can be written in a unique way as the sum of an element 
of M and an element of and that moreover M n J{M^) = {0}. As shown In [3], this is actually even a direct 

sum in Birkhoff-James’ sense, that is ||u||LP(n) < \\u + tu||LP(o) for all u G M, v G J{M^) and t G R. 
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Proof. Let / £ LpQ^.(r2) be fixed. For all n,k > 1, define the space 

Xn,k ■={cf>€ W^’P{Bn) : IIV(/.||lp(B„) <k, [ 0 = 0 

I JBr, 

endowed with the weak topology. By Poincare’s inequality and by the Banach-Alaoglu theorem, this 
space is metrizable and compact, hence Polish. Consider the multifunction P^^^ : n ^ Xn^k defined 

by 

rn,fc(w) := {0 G Xn^k ■■ V0|s„ = /(•,w)|b„}- 

Clearly Tn^ki^) is closed for all oj. We first prove further properties of this multifunction, and the 
conclusion will then follow by applying the Rokhlin-Kuratowski-Ryll Nardzewski theorem. 

Step 1. For all n > 1, we claim the existence of an increasing sequence of events i^n,k *G n such that 
]P[^n,A;] t 1 S'® ^ t oo for hxed n, and such that r„^fc(a;) jtz 0 for all ui £ i^n,k- 

By Definition 3.5, there is a subset D' C D of maximal probability, such that for all u € kl' the 
function is a potential field in M'^), and hence there exists 0“^ £ such that 

= V0‘^. For all n > 1, define 0^ := (fP — (lP‘{z)dz £ By definition, Jg 0^ = 0 

and V0^ = on B^- Moreover, ||V0i;^||LP(5^) < M holds for all oo £ ^n,k, where we define the 

event 

^n,k ■= {w £ D' : \\f{-,io)\\LP(Bn) — ^}- 
Integrability and stationarity of / easily imply that f 1 as /c f oo. 

Step 2. Proof that Tn^k is measurable, in the sense that £ J- for all open subset O C Xn,k, 

where we have set 

r"fc(0) := {a; € D : r„,fc(a;) n O / 0 }. 

As Xn,k is metrizable, it suffices to check £ J- for all closed subset F C Xn^k (see e.g. [5, 

Lemma 18.2]). Given a closed subset F C Xn^k^ we may write, using Poincare’s inequality and the 
weak lower semicontinuity of the norm, 

KUP) = {o; G D : 30 £ F, V0 = fi;u;)\Bj 

OO 

= P|{a; £ D : 30 £ F, ||V0 - /(•,a;)||LP(B„) < 1/i}- 
i=i 

Separability of the Polish space ^ implies that F is itself separable, and there exists a countable 
dense subset Fq C F. Hence 

OO 

= n u ^ ^ ^ - /(g^)IIl-(b„) < 1/j}, 

j=14>&Fo 

and measurability of r~^(F) then follows from measurability of /. 

Step 3. Conclusion. 

By steps 1 and 2, for all n,k > 1, the restricted multifunction : ^n,k ^ X^^k is measurable 

and has nonempty closed values. As Xn^k is a Polish space, we may apply the Rokhlin-Kuratowski- 
Ryll Nardzewski theorem (see e.g. [5, Theorem 18.13]), which gives a measurable function (j)n,k '■ 
Q.n,k -t Xn^k such that 4>n,k{^) G rn,fc(w), that is V0„,fc(-,a;) = /(•,a;)|s„ for all lo £ Qri,k- For all 
n > 0, define a measurable function cfn : Pt, ^ VF^’^(Fn) by 

OO 

0n(a;) = loi(a;)0i,„(a;) ^ (a;)0„,fc(w). 

k=2 

By definition, we have V0n(-,c<;) = /(•,w)|_b„ for all oo £ where := U£i is a subset of 
maximal probability. Denote D" : = n:r=i^^n. 

Let re > 1. By definition, V0n —V0i vanishes on Bi, hence the difference 6 n{^o) ■= 0n(‘, w) —0i(-, w) 
is constant on Bi for all oo £ D" and defines a measurable function 6 n ■ D" —)• M. Then consider the 
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measurable function ^ ^ W^’^{Bn) defined by := 4>n{x,ui) — 5n{^)- By construction, 

for all m > n > 1, we have = V’m on Bn, so the V’n’s can be glued together and yield a measurable 
function such that V'0(-,w) = f{-,oj) for all oj £ 17". □ 

A.4.2. Sufficient conditions for Hypothesis 2.1. As will be shown below, the measurability Hypoth¬ 
esis 2.1 is automatically satisfied if the integrand is quasiconvex and has the following nice approxi¬ 
mation property, introduced by [6]: 

Definition A.11. A normal random integrand PH : x x 17 —>■ [0, oo] is said to be quasiconvex 

if W{y, -, 0 ;) is quasiconvex for almost all y,uj. Given p > 1, it is further said to be p-sup-quasiconvex 
if there exists a sequence {Wk)k of quasiconvex normal random integrands such that Wk{y,l^,aj) f 
W{y,K,uj) pointwise as k f oo, and such that, for all k, for almost all y,uj and for all A, A', 

^\A\P -C< Wk{y,A,io) < Ckil + |A|^>), (A.3) 

for some constants C, Ck > 0. □ 

Note that Tartar [40] has proven the existence of quasiconvex functions that are not p-sup- 
quasiconvex properties for any p > 1. Before stating our measurability result, let us examine some 
important particular cases: 

Lemma A.12. Let IT : x ^ [0,oo] be a normal random integrand. Assume that there 

exist (7 > 0 and p > 1 such that, for almost all oj,y and for all A, 

^\A\P-C<W{y,A,u). (A.4) 

Also assume that one of the following holds: 

(1) W = V is a convex normal random integrand, and the convex function M := sup ess^^^ H(y, -, 00 ) 
has 0 in the interior of its domain; 

(2) W satisfies Hypothesis 2.1, with a convex part V such that M := supess^^^ y{y, has 0 in the 
interior of its domain. 

Then, W is a p-sup-quasiconvex normal random integrand. □ 

Proof. Case (2) directly follows from the approximation result of case (1) applied to the convex part 
V. We may thus focus on case (1). Let W = V he a convex normal random integrand. For all A: > 0, 
consider the following Yosida transform: 

Vk{y,A,uj)= inf {V{y,A',uj) +k\A-A'\P). (A.5) 

A/gRmxd 

For almost all y,uj, convexity of P4(p, -jw) easily follows from convexity of V{y, -, 00 ). For almost all 
y,oj, the lower semicontinuity of V{y,-,Lo) ensures that Vk{y,-,uj) f V{y,-,uj) pointwise as k f 00 . 
Moreover, by definition, Vk{y, A,oj) < M(0)-\-k\A\P, while the lower bound (A.4) implies P4(p, A,w) > 
^\A\P-C. 

It remains to check that the I4’s are normal random integrands. For almost all y,co, the function 
V{y, ■,uj) is convex and lower semicontinuous, hence it is continuous on its domain (not only on 
the interior). As by assumption 0 £ mi the set is a convex subset of maximal dimension, 

and hence points with rational coordinates are dense in The infimum (A.5) defining I 4 may 

thus be restricted to As a countable infimum, the required measurability properties follow. □ 

We now turn to the validity of the measurability Hypothesis 2.1 for p-sup-quasiconvex integrands. 

Proposition A.13. Let O C be a bounded domain, let (17, P) be a complete probability space, 
and let IF : x x 17 —)■ [0, 00 ] be a p-sup-quasiconvex normal random integrand for some p > 1 

(in the sense of Definition A.11). Given some fixed function f £ LP(I7; L^(0; M™^'^)), consider the 
random integral functional I : IF^’^(0;M™') x 17 ^ [0, 00 ] defined by 

I{u,u:)= / W{y,f{y,Ljj) + Vu{y),Ljj)dy. 

Jo 


(A.6) 
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Then, I is weakly lower semicontinuous on M”*). Moreover, for all weakly closed subsets 

F C or F C {u G : JqU = 0}, the function to i-A I{v,U}) is 

F-measurable. In particular, Hypothesis 2.1 is satisfied. □ 

Proof. For all k, define the approximated random functional X 17 —>■ [0, oo] by 

Ik{u,u!)= / Wk{y,f{y,uj) + Vu{y),uj)dy. 

Jo 

As the VFfc’s are nonnegative, monotone convergence yields that Ik f I pointwise. Moreover, for 
all k, and almost all w, the quasiconvexity and the upper bound (A.3) satisfied by Wk{-,-,co) imply 
the weak lower semicontinuity of Ik{-,uj) on (see [2]). As a pointwise supremum of 

weakly lower semicontinuous functions, we deduce that I{-,u}) is itself weakly lower semicontinuous 
on 

Combining the weak lower semicontinuity of Ik with the uniform coercivity assumption (cf. lower 
bound in (A. 3)) and with Poincare’s inequality, we easily conclude, for any weakly closed subset 
F C M"*) or F C {u G :f^u = 0}, 

lim inf Ik{v,uj) = inf I{v,uj). (A.7) 

fci^oo v£F v£F 

For all k, as Wk is quasiconvex hence rank-1 convex in its second variable (see [9]), the p-growth 
condition (A.3) implies the following local Lipschitz condition: for almost all y,co, for all A, A', 

\Wk{y,A,co) - Wk{y,A',u;)\ < Ck\A - A'|(l + jAr^ + |A'ri), 

for some constant Ck > 0. The Holder inequality then gives, for all u,v, for almost all w, 

\Ik{u,u;)-Ik{v,u;)\ < CfcC||V(u - u)||lp(o)(1 + ||Vu||[-'o) + \Mino) + ll/(•>^)llLHO))• 

This proves that the map Ik{-,io) is strongly continuous on kF^’^(0;M™'), for almost all io. Given a 
weakly (hence strongly) closed subset F C VFq^’^(0; M™’) or F C {u G W^’P{0; M™') ; J^u = 0}, strong 
separability of kF^’^(0; M"*) implies strong separability of F, so there exists a countable strongly dense 
subset Fq C F. Therefore, the map 

uj I—)• inf Ik{v,uj) = inf Ik{v,uj) 
veF veFo 

is F-measurable, and the conclusion follows from (A. 7). □ 

A.4.3. Measurable minimizers. We show that in the convex case (or more generally in the p-sup- 
quasiconvex case) the random functional (A.6) admits a measurable minimizer. For that purpose, we 
begin with the following useful reformulation of the Rokhlin-Kuratowski-Ryll Nardzewski theorem, 
which essentially asserts that the measurability of the infimum implies the measurability of minimizers. 

Lemma A.14. Let X be a Polish space, let (17, F, //) be a complete measure space, and let I : A x 17 —)• 
[0,oo]. Assume that 

(i) for all CO, I{-,oj) is lower semicontinuous on X; 

(a) for all CO, I{-,co) is coercive on X (i.e. the sublevel sets {u G X : I{u,co) < c} are compact for 
all c > 0); 

(Hi) for all closed subset F <Z X, the map fip : 17 —>■ [0, oo] defined by (pp{co) := min^^p I(v, co) is 
F-measurable. 

Then, there exists an F-measurable map u : 17 —>■ A such that, for all oo G LI, 

I{u{co),co) = mmI{v,co) = 

v€X 


□ 
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Proof. By coercivity and lower semicontinuity, the minima of are always attained on all closed 

subsets F, so that the function (pp is always well-defined. 

Consider the multifunction F : ^ X defined by r(a;) = {tt € X : I{u,Ijj) = (/>x(w)}. By lower 

semicontinuity, r(w) C X is nonempty and closed for all oj. Moreover, we claim that T is measurable, 
in the sense that 

r-i(O) := {w e 0 ; r(w) no^0} 

belongs to T for all open subsets O C X. As X is metrizable, it actually suffices to check T~^{F) £ F 
for all closed subsets F d X (see e.g. [5, Lemma 18.2]). By the coercivity and the lower semicontinuity 
of /, we may write 


r ^(F) = {cj £ : 3u £ F, = (/>x(w)} 


Pi < o; £ n : Bn £ F, /(n, w) < 4>x{^^) + 

n=l ^ 

OO X 

Pi i a; £ n : < </>x(w) + 


1 


n 


n=l 


1 


where the right-hand side belongs to F, by measurability of cpx and (pp- Hence, F is measurable, 
and we may thus apply the Rokhlin-Kuratowski-Ryll Nardzewski measurable selection theorem (see 
e.g. [5, Theorem 18.13]), which states the existence of a F-measurable map u : H —>■ X such that 
u{ijj) £ F(a;) for all oj. □ 

Combining this measurable selection lemma with the measurability result of Proposition A. 13, we 
obtain the following: 

Proposition A.15. Let O he some bounded domain, let (H,F, P) be a eomplete probability space, 


id ^ ^mxd 

and p > 1 such that, for almost all y,uj and for all A, 


and let VF : x x H —)■ [0, oo] be a normal random integrand. Assume that there exists C > 0 


^\Af>-C<W{y,A,oo). 


(A.8) 


Also assume that W satisfies Hypothesis 2.1 and that is weakly lower semicontinuous on 

LF^’^(0; M”*) for almost all u (in particular this is the ease if W is convex or p-sup-quasiconvex 
in the sense of Definition A. 11). Given some fixed function f £ LP(H; L^(0;M”^^'^)), consider the 
random integral functional I : VF^’P(0;M™') x H ^ [0, oo] defined by 


I{u,uj)= W{y,f{y,u;) + Vu{y),u;)dy. 

Jo 


Then, for all nonempty weakly closed subsets F C LFg’^(0; M™) or F C {u £ LF^’^(0; 
there exists a F-measurable map u : Ll ^ F such that, for almost all oj, 

I{u{uj),uj) = inf I{v,uj). 

v&F 


): foU = 0}, 


□ 


Proof. Let X denote the Banach space VFq^’^( 0; M”^) or {u £ VF^’P(0;M™') : f^u = 0}, endowed with 
the weak topology, and, for all k > 1, consider the subset Xfc := {n £ X : ||Vn||LP(o) < A:}, endowed 
with the induced weak topology. By Poincare’s inequality and by the Banach-Alaoglu theorem, X^ is 
easily seen to be metrizable and compact, hence Polish. Let F C X be a nonempty (weakly) closed 
subset. 

Let Q' C fl denote a subset of maximal probability such that is weakly lower semicontin¬ 

uous on LF^’P(0;M™') for all ui £ Ll'. Since X^ C X is (weakly) closed, the intersection X^ n F 
is also (weakly) closed, and hence Hypothesis 2.1 asserts that the map oj —>■ inf„gFnAfc <^) is 
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J^-measurable. Applying Lemma A. 14 on the compact Polish space and on 14' then yields a 
J^-measurable map Uk ^ Xk such that, for all a; G 14', 

I(uk(uj),u}) = min I(v,io). 
vGFnXk 

The lower bound (A.8) and the triangle inequality give 

[ \Xu\P <C^2P-^+ C2P-^I{u,uj) + 2P-^ [ \f{-,uj)\P. 

Jo Jo 

Define for all A: > 1, 

14fc := jw G 14' : C^2P-^ + C2P-^ inf/(u, w) + 2^-^ ^ \f{-,uj)\P < A;| , 

and note that 14^ G by Hypothesis 2.1. By definition, for all a; G 14^ we have 

I(uk(uj),uj) = min I(v,uj) = mml(v,uj). 

vGFnXk v£F 

The sequence (14^)^ is increasing, 14^ f 14" := U^i ^k- By integrability of /, 14" C 14'" and P[14'" \ 
14"] = 0, where we have defined the event 14'" := {w G 14 ; I(v,uj) < oo}. Given some fixed 

w & F, the measurable map u :Vl ^ X defined by 

CX) 

u{uj) := ^ lQ^\Q^^_^Uk{uj) 


k=2 


satisfies by definition I(u{u}),uj) = I(v,uj) for all w G 14" U (14 \ 14" 


□ 


A.5. Approximation results. In the present appendix, we prove two general approximation results 
that are crucially needed in this paper. The first one is an extension of [36, Lemma 3.6] and [23, 
Proposition 2.6 of Chapter X]. 

Proposition A.16. Let O C be a bounded Lipschitz domain, which is also strongly star-shaped, 
in the sense that there exists xq G O such that 


—xo + O C q:(— xo + O), for all a > 1. 

Let 0 : [0, oo] be a convex lower semicontinuous function with 0 G int dom0, and let 

u G such that f^O(Vu) < oo. Then, 

(i) there is a sequence {vn)n C C°°(adhO;M'") such that Vx„ G intdom0 pointwise, 

Vn ^ u in and f Q{Vvn{y))dy ^ f Q{Vu{y))dy, 

Jo Jo 


(ii) there is a sequence {wn)n of (continuous) piecewise affine functions such that Xwn G 
int dom0 pointwise. 


\,mxd 


n 


Wr, 




and 


'U 


Q{Xwn{y))dy 


'O 


Q{Xu{y))dy. 


If in addition u belongs to for some 1 < p < oo, then the sequences {vn)n ond {wn)n 

can be chosen such that Vn ^ u and Wn u in Moreover, 

(a) if u belongs to VPq’^(0;M™), then we can choose Vn G and Wn\dO = 0 (and in that 

case the assumption that O be strongly star-shaped can be relaxed); 

(b) if O = Q and u G lTper(<5; M'"), then Vn and Wn can be both chosen to be Q-periodic; 

(c) ifF: [0,oo] is a (nonconvex) ru-usc lower semicontinuous function which is continuous 

on intdom0 and satisfies 0 < H < 0 pointwise, then the sequences {vn)n o,nd {wn)n can be chosen 
in such a way that f^E(Vvn) —t /^^(Vn) and E(Vwn} —)• /^^(Vu). □ 
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Proof. We divide the proof into four steps. 


Step 1. Proof of (i). 

First, we show that we can assume Vtt € intdom© almost everywhere. Indeed, assume that (i) is 
proven for such rt’s, and let us deduce the general case. Given u € with 0(Vu) < oo, 

we have Vtt € dom© almost everywhere, and hence by convexity tVn £ dom© almost everywhere for 
all t £ [0,1). As convexity also implies 0(tVn) < t ©(Vu) + (1 — t)0(O) < oo, we can apply the 
result (i) to tu, for any t £ [0,1): this gives a sequence {vn,t)n C C'°°(adhO; M”*) such that Vn,t tu 
in and QiVVn,t{y))dy —)• Jq Q{tVu{y))dy. Weak lower semicontinuity of the integral 

functional u i-?- ©(Vu) on (which follows from convexity and lower semicontinuity of 

©) implies liminfi-i-i 0(tVu) > /^©(Vrt). As the converse inequality follows from convexity, we 
obtain 


lim sup lim sup 

ttl nfoo 




[ G{VVn,t) - f 0(Vn) 
Jo Jo 


= lim sup 
tti 



'“lliyLi(o) + 



G{tVu) 


'o 


G{Vu) 


= 0 , 


and hence a standard diagonalization argument gives a sequence {vn)n C C'°°(adh O; R™) such that 
Vun £ intdom© almost everywhere, —)■ rt in R”*), and Jq G(Vvn) —)■ Jq G(Vu}, and 

proves the general version of (i). 

Hence, from now on, we assume Vtt £ intdom© almost everywhere. Moreover, without loss of 
generality, we may also assume that O is strongly star-shaped with respect to xq = 0. Choose 
(ak)k C (l,oo) a decreasing sequence of positive numbers converging to 1, sufficiently slowly so that 
^ ^ dist(adh O, (9(afcO)) for all k, and define 

Uk : dkO —)• R™' : x i-A Uk{x) = u{x/ak). 

Take p £ C'“(R'^) such that f p = 1, p > 0 and suppp C B(0, 1), and write pk{x) = k~^p{kx), for 
all k > 1. Consider the sequence {vk)k defined by Vk = Pk* (dkUk) (which is well-defined by virtue of 
the condition on the ak’s). Note that Vk £ C'°°(adhO;R™) and Vufc = Pk* (Vtt)^, where we use the 
notation (Vtt)fc(x) = S/u{x/ak). Moreover, we then observe Vk ^ u m. kF^’^(0;R™') since Uk ^ u et 
(Vtt)fc —)• Vtt in L^(0;R”^). Now, Jensen’s inequality yields 


0 < ©(Vufc) = G{pk * {Vu)k) < Pk* (©((V'u)fc)) = Pk* (G(Vu))k. 

As the sequence (pk * (G(Vu))k)k converges to © o Vu in L^(0;R"*), it is uniformly integrable, and 
the same thus holds for the sequence (G(Vvk))k- As Vvk —)• Vu in L^(0;R™'), the convergence 
holds almost everywhere up to an extraction. Since by convexity © is continuous on intdom©, since 
Vu £ intdom© almost everywhere, and since Vvk —>■ Vu almost everywhere up to an extraction, we 
deduce G{Vvk) —>• ©(Vu) almost everywhere up to an extraction. By uniform integrability the latter 
convergence also holds in L^(0;R™'), which allows us to get rid of the extraction. In particular, 

[ G{Vvk)^ [ ©(Vu). 

Jo Jo 

Moreover, (Vu)fc £ intdom© almost everywhere, and thus Vvk = Pk* (Vu)^ £ intdom© everywhere 
for all k, since intdom© is a convex set containing 0. This proves part (i). 

Let us now assume that u £ iy^’P(0;R™') for some 1 < p < oo, and consider the sequences (ufc)^, 
{{Vu)k)k {vk)k defined above. First, Lebesgue’s dominated convergence theorem implies Uk ^ u 
and {Vu)k —t Vu in L^(0;R”*). Further, since the Jensen inequality gives 




\akPk *Uk -u\P < 



Pk{t)\akUk{x - 


t) — u{x)\^dt dx, 


and likewise for gradients, we conclude that the sequence {vk)k converges to u in iy^’^(0; R™). This 
proves part (i) in the case when u £ iy^’P(0;R™). 


Step 2. Proof of (ii). 
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For all k, since Vk G C'°°(adliO; M”*), there exists a sequence {wkj)j of piecewise affine functions 
such that Wk,j —)■ Vk in M”^) as j t oo and ||Vrcfcj||L°° < ||Vufc|| for all j, k (see e.g. [23, 

Proposition 2.1 of Chapter X]). Further, these functions Wkj can (simply remember their construction 
by triangulation) be chosen taking their values in int dom0, since we have constructed Vufc G int dom0 
everywhere for all k. Another approximation argument further allows us to choose Wkj such that 
Vrcfcj- only takes rational values. The desired result then follows from Step 1 and a diagonalization 
argument. Finally, the particular case when u belongs to VF^’^(0;M™') is obtained similarly as in 
Step 1. 

Step 3. Proof of the additional statements. 

It remains to address the particular cases (a) and (b). First assume that u belongs to Wq’^{0] M™). 
For the corresponding result, we refer to [23, Proposition 2.6 of Chapter X]. The only difference is 
that the argument in [23] requires continuity of 0. Instead, we replace u by tu for t < 1 as in Step 1, 
so that by convexity tVu G int dom0 almost everywhere, hence all the constructed quantities have 
almost all their values in intdom0, on which 0 is continuous by convexity. No further continuity 
assumption is then needed. 

Finally, if we assume O = Q with u G VFper(Q; K™), then we can consider the periodic extension 
of u on and repeat the arguments in such a way that periodicity is conserved. 


Step 4. Proof in the nonconvex case. 

Let u G be such that 0(Vu) < oo, which implies Vu G dom0 almost everywhere. 

Let t G (0,1). The approximation result given by point (i) gives a sequence {un,t)n of smooth functions 
such that Un,t —^ tu (strongly) in LF^’P(0;M™’) and Q{Vun,t) —^ Jq &{tVu) as n t oo, and such 
that Vun^t G intdom0 pointwise. Up to an extraction, we have Vun,t tVu almost everywhere, 
and thus H(VM„^i) —)• H(fVri) almost everywhere, which follows from continuity of H on the interior 
of its domain, with indeed fVu G int domH almost everywhere. Then noting that 


0 < E{VUn,t) < C(1 + 0(VM„,i)), 


and invoking both Lebesgue’s dominated convergence theorem (for S(Vu„^t)) and its converse (for 
@{Vun,t)), we deduce convergence Jq E{Vun,t) —t Jq H(fVu) as re t oo. As H is lower semi-continuous 
and also ru-usc, we compute, by Fatou’s lemma. 


/ E{Vu{y))dy< / liminf H(fVre(y))dy < liminf / E{tVu{y))dy 

Jo Jo ifl Jo 

< limsup / E{tVu{y))dy < / E{Wu{y))dy. 

ttl Jo Jo 


Hence, 


lim lim 
1 , 


/ E{Vun,t) = lim / E{tVu) = [ H(Vre) 

Jo HI Jo Jo 


and similarly 

lim lim f QiVun^t) = [ 0(Vre), 

ifl n^oojQ ’ Jq 

so that the conclusion follows from a standard diagonalization argument. The other properties are 
deduced in a similar way. □ 


For technical reasons, we need in the proof of Proposition 3.10 to further approximate piecewise 
affine functions by refined ones with smoother variations. The precise approximation result we need 
is the following. 


Proposition A. 17. Let u be an -valued continuous piecewise affine function on a bounded Lips- 
chitz domain O C Consider the open partition O = [+j|T;^ associated with u (i.e. u is affine on 
each piece O^). Define M := (ljf=i dO^) \ dO the interior boundary of this partition of O, and, for 
fixed r > 0, also define Mj- := (M J- Br) 0 O the r-neighborhood of this interior boundary. Then, for 
all K > 0, there exists a continuous piecewise affine function on O with the following properties: 
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(i) = Vu pointwise on 0 \ Mr, and limsup^|o suPo<K<i — ^||l°°(0 ) = 0/ 

(a) S conv({Vti(x) : X € O}) pointwise (where conv(-) denotes the convex hull); 

(Hi) denoting by O := 1+jJ^^’' O^K,r open partition associated with u^^r, and Aj, ,, := ViiK,r|o* 

all I, we have |AJ. ^ — AK,r| < for all i,j with 90^ ^ n dOf^r 7^ 0- D 


Proof. Let u, O and r be fixed. Without loss of generality, we can assume 0 € O. Denote rg := 
dist(0,90) and Rq := max 3 ;ggo dist(0, x), and define > 1 by (ctr — l)i2o = '^/2. Choose a 
nonnegative smooth function pr supported in i?(a^_i)ro with f pr = 1, and consider the smooth 
function Ur on O defined by Ur = Pr * [oiru{-/Ur)]. Since by definition inequality \arX — x| < 
[ur — l)i2o = holds for any x G O, we easily check Vri^ = Vm on the set O \ Mr- As Ur is 
smooth, we can consider Lr := max^-gadhO |VVnr(x)| < oo. Choose a triangulation of O 

which is a refinement of the partition {O^ \ Mr : I = 1,..., A;} U {Mr}, such that the diameter of 
each of the i® most K/{2Lr). Now construct as usual the piecewise affine approximation 

UK,r of Ur with respect to the triangulation {0^,.r)l (®oo o-S- [23, Proposition 2.1 of Chapter X]). By 
construction, Vu^^r = Vit,. = Vu on 0\Mr, since Ur is affine on each connected component of 0\Mr. 
Also note that the construction ensures that VuK,,r belongs pointwise to the set {Vrir.(x) : x G O} 
(as a consequence of the mean value theorem), which is by dehnition included in the convex hull 
conv({Vu(x) : x G O}). Moreover, if dO^r C dOf^r 7^ 0, it implies that r C has diameter- 
bounded by n/Lr- The construction of u^^r then gives 


Vu 


K.r O® 

' K,r 


- Vm 


< sup 
xeo 


sup 

K,T y&of , 


|VMr(x) — VUr{y)\ < — sup 
Ar xGO 


VVrir(x)| < K, 


which proves property (iii). Finally, the last property of (i) directly follows from the construction. □ 


To treat the case of periodic boundary conditions, we further need the following periodic version 
of the previous result. The proof is omitted, since it is a direct adaptation on the torus of the proof 
above. 


Proposition A.18. Let u he an W^-valued continuous piecewise affine Q-periodic function on 
Consider the periodic partition = l+J^^ O* associated with u (i.e. u is affine on each piece O^). 
Define M := dO^ the boundary of this partition, and, for fixed r > 0, also define Mr := 
(Af + Br) n Q the r-neighborhood of this boundary. Then, for all k > 0, there exists a continuous 
piecewise affine Q-periodic function u^^r on with the following properties: 


(i) VuK,r = Vri pointwise on \ Mr, and limsup^^g ®uPo<K<i \WK,r — i^||l°°(Q) = 0 / 

(it) ^UK,,r £ conv({Vu(x) : x G Q}) pointwise (where conv(-) denotes the convex hull); 

(iii) denoting by := 1+)^;^ ^ the Q-periodic partition associated with UK,r, and A(. ^ := Vtt 

for all I, we have |A),,, — Mn^r\ < for all i,j with dT^ .,. n dTfr 7^ 0- 


□ 
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